CSSANB 


A  CODE  FOR  SYSTEM  SURVIVABILITY 
ANALYSIS  -  NUCLEAR  EFFECTS 
THESIS 


GNE/PH/72-3 


Robert  G. 
C,  ptain 


DeRaad 

USAF 


D  D  C 

sug 


Approved  for  public  release;  distribution  unlimited. 


Unclassified 


Security  Classification  * 


DOCUMENT  CONTROL  DATA  -  R  &  D 

(Secur/fy  classification  of  title,  body  of  abstract  and  Indexing  annotation  must  bn  entered  when  the  overall  report  Is  classllled) 


I •  0R1G:nATING  ACTIVITY  (Corpora?*  author)  j2«.  REPORT  SECURITY  CLASSIFICATION 


Air  Force  Institute  of  Technology 
Wright-Patterson  Air  Force  Base,  Ohio  4S432  l2b- MOUP 


Unclassified 


3.  REPORT  TITLE 


CSSANE,  A  CODE  FOR  SYSTEM  SURVIVABILITY  ANALYSIS  -  NUCLEAR  EFFECTS 


4.  DESCRIPTIVE  NOTES  (Typo  o!  report  and  inclusive  dates) 

AFIT  Thesis 


5.  AUTHORiS)  (Pirst  name,  middle  Initial,  last  name) 

Robert  G.  DeRaad 
Captain  USAF 


ft.  REPORT  OATE 

June  1972 


la,  TOTAL  NO.  J.F  PAGES 


176,  NO.  OF  REFS 


9a.  ORIGINATOR'S  REPORT  NUMBERS) 


GNE/PH/72-3 


9b.  OTHER  REPORT  no(S)  (Any  other  numbers  that  may  be  assigned 
this  repo.*!) 


to.  DISTRIBUTION  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


‘A|)pfbVgaTWf'P51ic  release;  IAW  AFR  19l)’-2isroNSOR,i4<5'M,L,TARY  ACT,v,Ty 
1 4, -ft;  A  L'lSSQ,^^  ;•  \  1 

vKelflrA.  Williams,  1st  Lt. ,  USAF 

A&tine:  Director  of  Information 


13.  ABSTRACT  f 


\.  ■  0 


A  computer  code  has  been  written'^to  determine  aerospace  system 
survivability  to  selected  nuclear  effects.  A  special  case  of 
survivability  is  treated  in  which  system  survival  is  based  on  a 
comparison  of  the  sure-kill  vulnerability  level  and  the  computed  free 
field  nuclear  effects  levels.  The  code  consists  of  two  functional 
parts.  One  part  is  the  guiding  program  which  conducts  the  surviv¬ 
ability  study;  the  other  part  consists  of  individual  subroutines 
which  evaluate  free  field  nuclear  environment  levels.  Subroutines 
have  been  included  to-  evaluate  the  blast  and  thermal  environments.^  > 
At  a  later  date,  subroutines  will  be  added  to  evaluate  the  effects  u 
levels  for  x-ray,  gamma  ray,  neutrons,  and  EMP .  The  present  code  i s \\ 
capable  of  handling  from  one  to  ten  nuclear  bursts  and  up  to  100 
similar  vehicles  in  a  single  program  run.  Any  type  of  system  may 
be  treated  for  which  the  effects  vulnerability  levels  are  known. 

The  code  has  been  written  in  the  FORTRAN  EXTENDED  language  for  a 
CDC  6600  computer  with  a  scope  3.3  compiler.  Central  core  memory 
required  is  approximately  40,000  octal  words  and  run  times  are  on  the 
order  of  seconds.  A  sample  problem  has  been  included  to  illustrate 
the  type  of  study  that  may  be  performed  and  to  demonstrate  the  use 
of  the  program. 


m. 


FORM 
1  nov  es 


1473 


Security  Classification 


GNE/PH/72-3 


CSSANE 

A  CODE  FOR  SYSTEM  SURVIVABILITY 
ANALYSIS  -  NUCLEAR  EFFECTS 

THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 

by 

Robert  G.  DeRaad,  B.S.M.E. 

Captain  USAF 

Graduate  Nuclear  Engineering 

June  197? 

Approved  for  public  release;  distribution  unlimited. 


GNE/PH/72-3 


Preface 

The  computer  code  presented  in  this  report  is  called  a 
code  for  system  survivability  analysis.  It  is  designed  to 
serve  as  a  tool  for  engineering  studies  of  system  survival 
rather  than  operations  analysis,  though  the  latter  interest 
may  be  served  as  well. 

The  following  analogy  is  given  to  define  and  clarify 
the  term  "survivability  code"  as  used  in  this  report.  As 
the  human  body  is  vulnerable  in  varying  degrees  to  many 
diseases,  so  is  the  aerospace  system  vulnerable  in  varying 
degrees  to  each  of  the  nuclear  effects.  The  survival  of  the 
body  is  dependent  partially  on  the  intensity  of  exposure  to 
disease.  Exposure  levels  are  often  related  to  the  environment. 
Similarly,  survival  of  systems  to  the  threat  of  nuclear 
effects  depends  on  the  system  vulnerability  and  on  the  threat 
level.  This  code  evaluates  a  rather  special  case  of  surviv¬ 
ability  in  which  the  survival  determination  is  based  on  a 
sure-kill  vulnerability  level  and  the  computed  threat  level. 

I  hope  that  this  program  will  help  to  emphasize  the 
individuality  of  each  of  the  nuclear  effects  and  the  variation 
in  the  relative  importance  of  each,  depending  on  the  burst 
and  receiver  positions. 

I  wish  to  thank  Captain  George  Kuch  and  Mr.  Alfred  Sharp 
of  the  Air  Force  Weapons  Laboratory  at  Kirtiand  AFB  for  their 
assistance  in  providing  information  on  the  thermal  and  blast 
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effects.  The  guidance  and  suggestions  of  my  thesis  advisor. 
Dr.  Charles  J.  Bridgman,  are  also  gratefully  acknowledged. 
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Abstract 

A  computer  code  has  been  written  to  determine  aerospace 
system  survivability  to  selected  nuclear  effects,  A  special 
case  of  survivability  is  treated  in  which  system  survival  is 
based  on  a  comparison  of  the  sure-kill  vulnerability  level 
and  the  computed  free  field  nuclear  effects  levels.  The 
code  consists  of  two  functional  parts.  One  part  is  the 
guiding  program  which  conducts  the  survivability  study;  the 
other  part  consists  of  individual  subroutines  which  evaluate 
free  field  nuclear  environment  levels.  Subroutines  have 
baen  included  to  evaluate  the  blast  and  thermal  environments. 
At  a  later  date,  sub' cutines  will  be  added  to  evaluate  the 
effects  levels  fs-.  ra y,  gamma  ray,  neutrons,  and  EMP.  The 

present  code  :s  capable  of  handling  from  one  to  ten  nuclear 
bursts  and  up  to  100  similar  vehicles  in  a  single  program 
run.  Any  type  of  system  may  be  treated  for  which  the  effects 
vulnerability  levels  are  known.  The  code  has  been  written  in 
the  FORTRAN  EXTENDED  language  for  a  CDC  6600  computer  with  a 
&cope  3.3  compiler.  Central  core  memory  required  is  approxi¬ 
mately  40,000  octal  words  and  run  times  are  on  the  order  of 
seconds.  A  sample  problem  has  been  included  to  illustrate 
the  type  of  study  that  may  be  performed  and  to  demonstrate 
the  use  of  the  program. 
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CSSANE 

A  CODE  FOR  SYSTEM  SURVIVABILITY 
ANALYSIS  -  NUCLEAR  EFFECTS 

I .  Introduction 

The  effects  of  a  nuclear  environment  on  aerospace 
systems  is  an  important  factor  in  systems  analysis.  The 
system,  the  threat,  and  the  flight  regime  all  enter  into  the 
analysis.  Nuclear  weapons  are  of  special  interest  because 
of  the  variety  of  effects  that  pose  a  threat.  Additionally, 
the  relative  importance  of  each  threat  varies  with  the  burst 
altitude  and  the  flight  regime. 

The  objective  of  this  study  was  to  develop  a  computer 
code  to  evaluate  the  survivability  of  systems  in  a  nuclear 
environment.  There  are  existing  codes  that  include  surviv¬ 
ability  analysis  as  part  of  an  operations  analysis  investiga¬ 
tion  and  there  are  codes  that  evaluate  the  free  field  nuclear 
effects  levels  without  regard  to  system  survivability.  This 
code  conducts  a  "special  case"  survivability  analysis  for 
systems  which  are  undergoing  a  specified  nuclear  attack. 

The  code  evaluates  levels  of  selected  nuclear  effects  and 
evaluates  system  survival  based  on  a  comparison  of  those 
levels  to  sure-kill  vulnerability  levels  which  are  specified 
by  the  user.  Information  furnished  by  the  code  includes  the 
relative  level  of  the  various  effects  in  addition  to  surviv¬ 


ability. 
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The  computer  program  consists  of  a  routine  named  GUIDE 
which  accepts  data  input  and  conducts  the  study,  and  a  set  of 
subroutines  that  computes  effects  levels  and  evaluates  the 
survival  of  systems.  Initially,  subroi  tines  have  been 
included  to  evaluate  both  blast  and  thermal  environments. 

This  type  of  organization  was  chosen  so  that  the  program 
could  be  easily  expanded  to  include  subroutines  for  x-ray, 
gamma  ray,  neutrons,  and  EMP. 

Chapter  II  of  this  report  describes  the  computer  program 
in  detail.  The  description  includes  capabilities,  limita¬ 
tions,  and  options  available  to  the  user.  Chapter  III 
describes  the  treatment  of  blast  effects,  the  subroutine 
for  blast  computations,  and  the  capabilities  of  the  sub¬ 
routine.  Similar  information  concerning  thermal  effects  is 
presented  in  Chapter  IV.  Chapter  V  presents  the  results  and 
conclusions  including  recommendations.  A  statistical  model 
for  obtaining  a  normally  distributed  variable  for  targeting 
purposes  is  presented  in  Appendix  A.  Appendix  B  includes 
a  list  of  major  program  variables  along  with  definitions 
and  the  location  where  each  variable  appears  first  in  the 
code.  Appendix  C  contains  detailed  instructions  for  making 
the  required  data  input.  Appendix  D  is  a  source  listing  of 
the  complete  code  and  Appendix  E  presents  the  input  and 
output  of  a  sample  problem. 
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I I .  The  Survivability  Analysis  Code 

A  study  of  system  survivability  to  nuclear  effects  is 
a  function  of  the  system,  the  threat,  and  the  flight  regime. 

A  survivability  analysis  requires  a  knowledge  of  the  vulnera¬ 
bility  of  a  system  and  of  the  threat  level.  Vulnerability 
has  been  defined  as  the  level  of  a  specified  effect  that  will 
cause  a  finite  degradation  in  the  capability  of  a  system  to 
perform  its  mission.  Survivability  is  defined  as  the 
capability  of  a  system  to  accomplish  a  specified  mission 
while  withstanding  the  effects  of  a  nuclear  environment. 

For  the  purposes  of  this  study,  system  survivability  is  based 
on  a  comparison  of  the  computed  threat  levels  to  specified 
"sure-kill"  vulnerability  levels. 

A  description  of  the  code  that  was  developed  to  perform 
the  survivability  study  is  presented  in  this  chapter.  In 
addition  to  the  description,  the  capabilities,  limitations, 
options,  data  input,  and  the  output  of  the  code  are  discussed. 

Description  of  the  Computer  Program 

The  code  is  named  CSSANE,  "A  Code  for  System 
Survivability  Analysis  -  Nuclear  Effects".  It  has  been 
developed  to  provide  a  simple  and  easy-to-use  method  for 
determining  survivability  of  systems  in  a  well  defined  nuclear 
threat  situation.  This  determination  is  based  on  sure-kill 
vulnerability  levels  which  are  specified  by  the  user  for  each 
nuclear  effect  threat  under  consideration.  The  free  field 


nuclear  effect  environments  are  calculated  in  subroutines  and 
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comparisons  are  then  made  to  the  specified  vulnerability 

This  comparison  results  in  a  kill,  no-kill  decision. 
This  type  of  an  analysis  does  not  lend  itself  to  operations 
*  r  a  1 y  s is  but  is  rather  intended  as  an  aid  in  engineering 
studies  of  system  survivability. 

The  code  consists  of  two  functional  parts.  The  main 
program,  GUIDE,  controls  or  directs  the  survivability  study 
ail  is  one  part.  The  second  part  of  the  program  consists  of 
a  set  of  individual  subroutines  which  evaluate  each  nuclear 
threat  level. 

Data  required  by  the  code  are  entered  in  program  GUIDE. 
This  includes  the  yield  and  the  time  of  burst  for  one  or  more 
nuclear  detonations;  the  position  in  space  and  the  velocity 
vectors  for  one  or  more  similar  vehicles;  and  the  vulnera¬ 
bility  levels  of  those  vehicles  to  each  of  the  nuclear 
effects  which  are  to  be  studied. 

The  burst  locations,  the  incremental  time  to  each 
succeeding  burst,  the  time-updated  vehicle  positions  and 
velocities,  the  vulnerabilities,  and  the  slant  ranges  from 
burst  to  vehicle  position  are  maintained  in  GUIDE.  In  addi¬ 
tion,  the  subroutines  which  evaluate  the  effects  levels  are 
called  from  GUIDE. 

Perhaps  the  best  way  to  demonstrate  the  intended  use  of 
the  program  and  its  capabilities  is  to  illustrate  a  typical 
aerospace  system  -  nuclear  weapon  encounter.  A  sketch  which 
describes  the  physical  situation  and  most  of  the  input 

> 

>  parameters  that  are  required  is  included  as  Pig.  1. 


A 
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Figure  1  describes  a  situation  in  which  four  similar 
vehicles  with  common  vulnerability  numbers  are  subjected  to 
two  separate  nuclear  bursts.  Each  of  the  vehicles  has  been 
assigned  initial  x,  y,  z  coordinates  in  an  arbitrary,  sea 
level  reference  and  velocity  components  at  time  T  *  0.  The 
ground  level  and  the  maximum  altitude  of  a  uniform  haze  layer 
over  the  earth  must  be  specified.  Both  of  these  values  may 
be  zero.  A  visibility  factor  in  miles  must  be  specified  for 
the  haze  cover. 

The  nuclear  threat  is  provided  by  two  bursts,  the  first 
of  which  is  intended  to  simulate  an  indirect  attack;  the 
second,  to  simulate  a  direct  attack  on  vehicle  number  one. 

The  indirect  attack  consists  of  a  large  yield  weapon  which 
is  detonated  at  tine  T  •  0  and  at  position  x,  y,  z  as  speci¬ 
fied  by  the  user.  The  direct  attack  (burst  nunber  two)  is 
specified  by  yield  and  tine.  Its  position  is  automatically 
determined,  based  on  the  target  point  which  is  the  updated 
vehicle  position,  and  a  spherical  probable  error  which  nust 
be  specified  for  a  given  weapon  delivery  systen. 

The  progran  is  cycled  through  its  evaluation  procedures 
for  each  burst.  In  this  cycle,  the  vehicle  positions  are 
updated  to  the  tine  of  each  burst,  the  nuclear  effects  levels 
of  interest  are  computed  for  each  vehicle,  and  the  survival 
determination  is  made  on  each  vehicle.  A  record  is  kept  of 
the  surviving  vehicles  and  printouts  are  nade  for  each  burst, 
noting  the  vehicles  which  have  been  destroyed,  the  nuclear 
effect  responsible,  and  the  levels  of  threats  encountered  by 
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each  vehicle.  Once  a  vehicle  has  been  designated  as 
destroyed,  it  is  removed  from  the  program  and  no  further 
calculations  are  performed  on  it. 

Program  Options 

There  are  five  options  in  this  code  and  they  are  speci¬ 
fied  by  data  input  to  GUIDE.  The  most  important  option 
determines  which  effects  are  to  be  considered  in  the  sur¬ 
vivability  analysis.  The  option  for  evaluation  of  any  effect 
is  chosen  by  entering  a  positive  value  for  the  vulnerability 
level  of  that  effect.  The  other  four  options  allow  a  choice 
of  a  nonstandard  atmosphere,  a  choice  of  method  of  burst 
placement,  a  choice  of  the  type  of  output,  and  a  choice  of 
employing  reactive  turning  (dodging  the  burst) .  The  method 
for  choosing  these  options  is  described  in  Appendix  C, 
"Details  of  Data  Entry". 

The  ARDC  standard  atmosphere  may  be  modified  by  speci¬ 
fying  the  temperatures  at  sea  level,  burst  height,  receiver 
height,  and  the  terrain  level.  Burst  placement  may  be 
specified  by  the  user  through  data  input,  or  may  be  obtained 
from  subroutine  TARGET.  Output  may  consist  of  a  listing  of 
all  values  of  each  effect  parameter  at  each  vehicle  plus  the 
survivability  result,  or  it  may  be  reduced  to  just  the  sur¬ 
vivability  result. 

Capabilities 

The  code  is  capable  of  handling  from  one  to  ten  nuclear 
bursts  if  they  are  placed  automatically  by  subroutine  TARGET. 
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Any  number  of  bursts  may  be  entered  if  the  burst  positions 
■>  arc  specified  by  the  user.  There  is  no  restriction  on  mixing 

the  mode  of  lurst  placement  except  that  the  number  of  bursts 
placed  automatically  cannot  exceed  ten  for  any  given  program 
run.  Each  burst  may  vary  in  yield  and  may  be  entered  with 
any  time  increment  desired.  The  bursts  are  entered  in  time 
sequence  and  each  burst  is  treated  separately  even  if  two 
bursts  are  entered  at  coincident  times. 

Any  number  from  1  to  100  vehicles  can  be  entered  for  a 
single  program  run.  They  may  be  traveling  in  different 
directions,  at  different  altitudes,  and  at  different  speeds. 
As  a  result,  it  is  possible  to  evaluate  systems  that  are 
operating  in  unrelated  flight  regimes  without  rerunning  the 
‘  program.  Any  type  of  system,  whether  in  the  air  or  station¬ 
ary  on  the  ground,  may  be  treated  provided  its  vulnerability 
is  known. 

There  is  a  capability  to  simulate  the  direct  attack  by 
calling  subroutine  TARGET  which  places  the  burst  in  a  normal 
distribution  about  the  target  point.  This  type  of  placement 
allows  for  a  margin  of  error  inherent  in  the  guidance  system 
of  the  weapon  delivery  vehicle.  That  error  is  generally 
described  as  the  spherical  probable  error,  spe.  The  spe  is 
the  radius  of  a  sphere  about  the  target  point  which  should 
contain  approximately  50%  of  the  burst  positions.  The  target 
point  is  always  chosen  as  the  first  filled  slot  in  the 
vehicle  array. 

-t  % 
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A  statistical  approximation  often  used  in  weapon  simula¬ 
tion  work  is  used  to  obtain  the  normally  distributed  x,  y,  z 
set  of  coordinate  points.  These  coordinates  then  represent 
the  burst  position.  A  development  of  the  approximation  is 
given  in  Appendix  A. 

If  the  burst  position  is  entered  as  input,  it  is  possi¬ 
ble  to  change  the  velocity  vectors  of  one  or  all  of  the 
vehicles  at  the  time  of  burst  detonation.  A  change  of  the 
velocity  vectors  away  from  the  burst  point  can  roughly  simu¬ 
late  reactive  turning.  Also,  for  each  buret,  to  eliminate 
unnecessary  calculations  each  vehicle  is  tested  for  inter¬ 
cept  with  the  fireball  and  for  intercept  with  the  ground. 
Vehicles  found  in  these  situations  are  annotated  as  destroyed 
and  are  removed  from  the  program. 

Subroutines  have  been  included  to  evaluate  survival  to 
blast  and  thermal  effects.  Subroutine  BLAST1  makes  survival 
tests  based  on  dynamic  pressure  and  overpressure  levels; 
subroutine  THERM  investigates  survival  to  thermal  fluence 
levels . 

Limitations 

There  are  some  limitations  that  exist  because  of  input 
data  that  cannot  be  changed  once  it  is  entered.  The  vulnera¬ 
bilities  are  not  variable  within  a  run,  and  only  those 
effects  are  evaluated  for  which  a  vulnerability  level  is 
specified.  Consequently,  it  is  not  possible  to  vary  the 
effects  which  will  be  evaluated  during  a  given  run.  A  single 


vulnerability  tor  each  effect  also  means  that  only  one  type 
of  vehicle  can  be  treated  per  run.  In  addition,  it  is  not 
easy  to  alter  the  flight  regime  of  a  vehicle  since  its 
heading,  speed,  and  position  are  non-variable  input.  (The 
reactive  turning  capability  excepted.) 

The  flight  regime  may  restrict  the  use  of  certain 
effects  subroutines,  primarily  because  of  maximum  altitude 
limitations  on  some  of  the  calculations.  The  capabilities 
and  restrictions  of  the  effects  subroutines  are  given  in 
Chapters  III  and  IV.  Also,  the  order  in  which  the  different 
effects  are  evaluated  is  determined  by  the  program,  not  by 
the  input. 

Data 

Data  Input .  All  input  values  that  are  required  are 
entered  into  the  program  GUIDE.  There  are  a  total  of  ten 
input  statements.  The  first  six  are  executed  once  for  each 
run.  The  last  four  are  executed  once  for  each  burst  in  a 
run . 

The  following  is  a  word  description  of  the  content  of 
each  of  the  first  six  input  statements: 

1.  This  "read”  enters  the  visibility,  the  water  vapor 
pressure,  the  albedo  (ground  surface  reflectance), 
and  altitude  of  the  haze  layer  (one  input  card 
required) . 

2.  Enters  the  number  of  vehicles  and  the  number  of 
bursts  to  be  studied  (one  input  card  required). 
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3.  Enters  the  values  associated  with  choice  of  a  non¬ 
standard  atmosphere  (one  input  card  required). 

4.  Sequentially  enters  the  values  for  coordinate 
positions  and  velocities  of  each  vehicle  (one  input 
card  required  for  each  vehicle  originally  entered). 

5.  Enters  the  list  of  integers  that  are  used  to  pre¬ 
cycle  the  random  number  generator  required  fo) 
targeting  (one  card  required) . 

6.  Enters  the  vulnerability  levels  for  survivability 
determinations  (one  card  required) . 

Data  that  must  be  entered  repeatedly  for  each  burst  is 
described  below: 

!.  This  input  statement  enters  the  incremental  time  of 
burst,  weapon  yield,  burst  height  above  ground,  and 
integers  for  the  options  of  burst  placement,  data 
output,  and  change  of  velocity  vectors  (one  card 
required) . 

2.  Enters  the  coordinate  points  of  burst  (one  card 
required) . 

3.  Enters  new  values  of  velocity  vectors  (one  card 
required  for  each  vehicle  originally  entered). 

4.  Enters  the  value  of  spherical  probable  error  (one 
card  required) . 

Note:  Specific  instructions  for  data  entry  are  given  in 
Appendix  C. 

Data  Output .  Hollerith  statements  are  used  to  supple¬ 
ment  and  clarify  the  output  data.  The  data  that  is  input  is 
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listed  and  identified  again  in  the  output.  The  vehicle 
positions  and  velocities  at  the  time  of  each  burst  are 
listed.  In  addition,  for  each  burst,  the  method  of  place¬ 
ment  and  the  burst  position  are  given. 

Results  of  the  survivability  tests  are  printed  sepa¬ 
rately  for  each  burst.  For  those  vehicles  that  are  destroyed, 
the  effect,  the  level  of  the  effect,  and  the  vulnerability 
level  are  all  given  in  a  removai  statement.  Removal  state¬ 
ments  are  also  printed  for  the  fireball  and  the  ground 
intercept  cases. 

Summary 

The  computer  code  consists  of  two  parts.  The  program 
GUIDE  accepts  data  and  directs  the  survivability  study;  the 
subroutines  compute  effects  levels  and  make  survival  analyses. 
This  organization  was  chosen  so  that  additional  subroutines 
for  other  effects  could  be  added  with  a  minimum  of  change. 

Subroutines  for  evaluation  of  blast  and  thermal  effects 
are  included.  Survivability  to  dynamic  pressure  and  to 
overpressure  levels  is  tested  by  subroutine  BLAST1;  surviva¬ 
bility  to  thermal  fluence  levels  is  tested  in  subroutine 
THERM.  The  program  allows  a  choice  of  the  effects  to  be 
evaluated,  method  of  burst  placement,  and  type  of  data 
output.  There  is  also  a  capability  to  specify  a  nonstMidurd 
atmosphere  and  to  simulate  reactive  turning. 

A  full  description  of  the  data  and  the  format  for  ,.ts 
entry  are  given  in  Appendix  C.  Capabilities  and  restrictions 
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for  each  of  the  effects  subroutines  are  given  in  the 
chapters  written  for  each  subroutine.  The  subroutine  and 
the  methods  for  calculation  of  blast  effects  are  presented 
in  Chapter  III.  Chapter  IV  contains  similar  information  for 
thermal  effects. 
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III.  Computation  of  Blast  Environments 
Blast  Phenomenon 

Blast  phenomenon  can  be  computed  using  analytic  or 
empirical  methods.  There  are  analytic  hydrodynamic  codes 
that  yield  extremely  accurate  shock  parameter  values,  but 
at  the  expense  of  lengthy  computer  runs.  Empirical  methods, 
on  the  other  hand,  scale  the  data  of  a  standard  explosion  to 
obtain  shock  values  for  other  specified  bursts.  One  such 
method  for  predicting  shock  phenomenon  is  the  Sachs  scaling 
technique.  The  blast  effects  information  presented  here  is 
based  on  Sachs  scaling  and  was  taken  directly  from 
AFWL-TR-70-85  (Ref  7).  The  report  includes  a  code  named 
SABER  which  evaluates  blast  parameters.  The  BLAST1  sub¬ 
routine  of  CSSANE  borrows  the  methods  of  SABER  to  find  blast 
effects  levels  for  given  receiver  positions  and  slant  ranges. 

A  capability  for  receiver  motion  was  added  to  allow  applica¬ 
tion  to  aircraft  in  flight.  The  calculations  of  shock 
parameters  include  methods  to  account  for  an  inhomogeneous 
atmosphere,  reflection  from  the  earth  (fused  shock  region), 
and  a  blast  efficiency  factor. 

Sachs  Scaling 

The  Sachs  equations  for  scaling  have  been  developed  in 
Ref  7,  pp.  2-12.  A  summary  of  these  equations  is  given  below. 

Overpressure  (AP)  psi 
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where  the  subscript  T  applies  to  the  explosion  under  investi¬ 
gation;  the  subscript  m  applies  to  the  1  KT  model,  and 
P^  »  ambient  pressure  (psi) 

E  »  energy  of  explosion  (KT) 

Ta  *  ambient  temperature  (/R) 

C  »  ambient  speed  rf  sound  (ft/sec) 

3 

p  b  ambient  density  (siugs/ft  ) 
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Effects  of  a  Nonhomogeneous  Atmosphere 

The  Sachs  scaling  equations  are  satisfactory  for  the 
coaltitude  burst  and  receiver.  However,  the  nonhomogeneous 
atmosphere  prevents  a  direct  application  in  the  noncoaltitude 
case.  That  case  may  be  handled  by  applying  the  Ledsham-Pike 
alpha  correction  (Ref  7,  12-14)  to  overpressure  calculations 
and  a  modified  Sachs  scaling  to  the  other  shock  parameters. 

The  alpha  correction  method  was  used  to  find  overpressure 
values  because  of  its  accuracy  and  because  the  overpressure 
is  a  value  that  is  subsequently  used  to  compute  other  shock 
parameters.  The  remaining  parameters  are  found  by  the  modi¬ 
fied  Sachs  scaling  method  which  assumes  the  blast  wave  is 
propagating  from  the  burst  to  the  receiver  under  ambient 
conditions  of  the  altitude  of  the  receiver.  Test  explosions 
have  verified  a  correlation  of  measured  effects  values  to 
those  obtained  by  use  of  the  alpha  and  modified  scaling 
methods . 

In  the  alpha  correction  method,  the  overpressure  at  the 
receiver  is  obtained  by  multiplying  the  overpressure  at  the 
burst  altitude  by  the  alpha  correction  factor. 

P  \  ^ 

apr  =  apb(p^)  (7' 

This  factor  is  the  ratio  of  the  pressure  at  the  receiver 
altitude  to  pressure  at  the  burst  altitude  raised  to  the 
power  alpha.  The  exponent,  alpha,  is  derived  from  empirical 
data  and  theoretical  calculations  end  is  a  function  of  the 
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scaled  range,  burst  to  receiver  (see  Fig.  2). 

The  modified  scaling  equations  for  the  noncoaltitude 
case  are  presented  here  in  Eqs  (8)  through  (14)  (Ref  7, 
22-26)  . 


Time  of  Shock  Arrival  sec 
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where 

t  _  *>  time  of  shock  arrival  for  model  (sec) 
am 

E^  a  energy  yield  of  burst  (KT) 
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V  w  *■-  ambient  pressure  at  receiver  (psi) 

Pcl  »  ambient  pressure  at  sea  level  (psi) 

C  =  ambient  speed  of  sound  at  receiver  (ft/sec) 

M 

Cgk  s  ambient  speed  of  sound  at  sea  level  (ft/sec) 

Note:  If  the  receiver  is  in  the  fused  shock  region  (Mach 

stem)  the  energy  term  must  include  the  amplifica¬ 
tion  factor  F  (see  section  on  reflected  shock  rein- 
r 

forcement  p.  19). 

The  Rankine-Hugoniot  relations  between  density  -  particle 
velocity,  particle  velocity  -  speed  of  sound,  and  dynamic 
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Fig.  2-,  Ledshain-Pike  Alpha  Ccrrectxon  as  a  Function 
cf  Scaled  Range  (Ref  7,  14). 
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pressure  -  particle  velocity  along  with  modified  scaling 
give  the  following  (Ref  7,  24-26): 

Shock  Front  Velocity  (Cg)  ft/sec 


/  6AP  \ 
Cs  =  Cz(1  +  JfT-J 
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where  AP  =  peak  overpressure  at  receiver 
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where  Pz  =  ambient  atmospheric  density  at  receiver 


Peak  Particle  Velocity  (U^)  ft/sec 


AP  .  .  6AP 


vK(pt)1 


7P 


♦  1 


) 


-1/2 


(13) 


Dynamic  Pressure  (q)  psi 
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Region  of  Reflected  Shock  Reinforcement 

Shock  phenomenon  for  bursts  near  the  ground  is  compli¬ 
cated  by  the  reflection  of  the  shock  front  off  the  ground 
into  the  region  of  previously  shocked  air.  The  reflected 
shock,  traveling  more  rapidly  in  the  heated  air,  overtakes 
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the  basic  shock  front  and  fuses  with  it  to  form  a  reinforced 
shock  region  called  the  Mach  stem.  The  triple-point  path, 
originating  at  a  point  on  the  earth  some  radial  distance 
from  ground  zero,  describes  a  boundary  below  which  there 
exists  a  region  of  shock  reinforcement. 

A  yield  amplification  factor,  F  ,  has  been  determined 

from  experimentation  for  both  the  region  near  the  triple- 

point  path  and  for  the  region  considered  far  within  the  Mach 

stem.  For  points  near  the  path,  the  yield  multiplying  factor 

F  ,  is  given  in  Fig.  3  as  a  function  of  scaled  altitude  of 
1/3 

burst  (Kft/(KTj  '  ) .  The  multiplying  factor  for  far  geometries 
i.e.,  points  which  are  at  scaled  ranges  greater  than 
100  from  the  triple  point  path,  is  given  by 

Fr  a  2.33  -  0.025  Rg  (15) 

a  function  of  the  scaled  range  (Rg)  in  Kft/(KT)1/3  (Ref  7,  64). 
Positive  Phase  Durations 

The  times  of  positive  phase  duration  for  overpressure, 
material  velocity,  and  overdensity  have  also  been  developed 
from  scaling  methods  (Ref  7,  26-31).  Empirical  relations 
have  been  derived  for  the  positive  duration  of  overpressure 
[ (AP) ^ ] m ,  the  ratio  of  scaled  positive  duration  of  material 
velocity  to  scaled  duration  of  overpressure  [Gt/(AP)^]^,  and 
the  ratio  of  scaled  positive  duration  of  overdensity  to 
scaled  duration  of  overpressure  [(Apj^/CAP)^^.  These 
empirical  relations  in  conjunction  with  the  general  equation 
for  positive  phase  durations,  Eq  (4),  yield  the  following: 


SCALED  HEIGHT  OF  BURST  [Kft/tKTr5] 


Fig.  3.  Yield  Amplification  (Reflection)  Factor  versus 
Scaled  Burst  Height  (Ref  7,  63). 
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Positive  Duration  of  Overpressure  (sec) 


(AP)t 


wv175 


(16) 


Positive  Duration  of  Material  Velocity  (sec) 


iymt],twt]. 

C2(F2/et)1/3 


(17) 


Positive  Duration  of  Overdensity  (sec) 
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Blast  Efficiency  Factor 

Scaling  methods  arc  very  *  seful  tools  for  finding  blast 
effects  lovels  but  there  are  limitations  to  be  considered. 

The  blast  phenomenon  is  hydrodynamic  in  nature  and  therefore 
the  methods  presented  here  are  not  useful  for  very  early 
times  when  the  dominant  mode  of  energy  transport  is  radiation. 
Also,  at  higher  altitudes  there  is  not  a  clear  separation  of 
the  shock  front  from  the  radiation  phase  before  much  of  the 
energy  is  dissipated.  A  blast  efficiency  factor  (Fig.  4) 
was  developed  from  data  obtained  by  runs  of  the  SPUTTER 
hydrodynamic  code  at  AFWL,  Kirtland  AFB,  New  Mexico  (Ref  7, 
31-32).  The  blast  efficiency  factor,  when  applied  to  the 
burst  yield,  can  extend  accurate  blast  calculations  from 
60,000  to  150,000  ft. 


ST  EFFICIENCY  FACTOR 
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Subroutine  for  Blast  Effects 

Subroutine  BLAST1 .  The  subroutine  for  calculation  of 
blast  parameters  and  for  evaluation  of  system  survivability 
to  blast  effects  is  BLAST1.  This  subroutine  is  called  by 
the  program  GUIDE  and  returns  to  GUIDE  a  record  of  those 
vehicles  that  have  been  destroyed  by  blast  effects.  Shod 
parameters  at  the  point  of  vehicle  and  shock  intercept  are 
computed  by  the  methods  presented  earlier  in  this  chapter. 
Survival  of  each  vehicle  is  based  on  the  vulnerability 
criteria  of  dynamic  pressure  and/or  overpressure. 

.  Parameters  Computed .  The  results  of  the  blast  computa¬ 
tion  for  each  burst  are  printed  from  the  BLAST1  routine. 

All  vehicle  losses  are  listed  under  "Results  of  Blast  Effects 
Computations."  The  lethal  level  and  the  type  of  pressure 
responsible  for  each  vehicle  loss  is  also  printed.  The 
following  is  also  printed,  out  if  the  option  for  data  is 
specified:  , 


Height  of  receiver  at  shock  intercept 

(ft) 

Slant  range  at  shock  intercept 

(ft) 

Time  of  shock  arrival 

(sec) 

Shock  front  velocity 

(ft/sec) 

Peak  dynamic  pressure 

(psi) 

Peak  overpressure 

(psi) 

Peak  material  velocity 

(ft/sec) 

Peak  overdensity 

(slugs/ft3) 

Positive  duration  of  overdensity 

(sec) 
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Positive  duration  of  overpressure  (sec) 

Positive  duration  of  material  velocity  (sec) 

Sub-Programs  Called  by  BLAST1 ♦  Sub-programs  TRIPNT, 
ATMOS,  SETUP,  MACURE,  and  MOTION  are  called  from  the  routine 
BLAST1.  Subroutine  TRIPNT  establishes  the  fused  shock  region 
for  bursts  near  the  surface  of  the  earth.  TRIPNT  also 
determines  a  yield  amplification  factor  for  vehicles  in  the 
fused  shock  region.  The  ATMOS  subroutine  provides  needed 
atmospheric  data  from  the  ARDC  standard  atmosphere.  SETUP 
and  MACURE  are  table  lookup  routines  that  are  capable  of 
interpolation  and  extrapolation  from  tables  of  data.  The 
tables  of  data  are  stored  in  the  BLOCK  DATA  section.  The 
programs  TRIPNT,  ATMOS,  SETUP,  and  MACURE  were  a  part  of  the 
SABER  code  and  are  described  in  AFWL-TR-70-85 ,  (Ref  7). 

Sub-program  MOTION  is  called  from  BLAST1  to  find  the 
point  of  vehicle,  shock  front  intercept.  Since  the  shock 
front  moves  outward  at  speeds  near  the  speed  of  sound,  the 
change  in  vehicle  position  during  shock  front  propagation  is 
important . 

The  model  for  finding  the  intercept  point  by  vector 
methods  is  presented  below.  Initially,  the  scaling  laws  are 
used  to  establish  the  shock  front  velocity  based  on  the 
vehicle  position  at  the  burst  time.  The  shock  and  vehicle 
positions  after  a  time  increment  AT  will  establish  a  closing 
rate 
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where 

CR  *  closing  rate 

AS  =  change  in  separation  distance,  (shock  to  receiver) 
AT  *  time  increment 

A  predicted  time  increment,  ATp,  for  intercept  is  then 
given  by 

AT  =  —  (20) 

where 

SD  =  separation  distance 

The  intercept  point  and  the  shock  parameters  are  finally 
obtained  through  an  iterative  process  involving  the  predicted 
time  increment.  The  no-intercept  case  will  be  indicated  by 
increasing  separation  distances  with  each  ;.<e  increment. 

Limitations  on  Blast  Calculations .  The  data  base  from 
which  blast  parameters  are  calculated  is  the  free-air,  sea- 
level,  1-KT  test  model.  The  range  of  altitudes  which  can 
be  spanned  by  scaling  techniques  is  limited.  The  application 
of  a  blast  efficiency  factor  extends  the  range  of  accurate 
calculation  from  60,000  ft.  to  approximately  150,000  ft. 

The  program  can  calculate  results  up  to  250,000  ft.  but  no 
verification  is  available.  Calculations  made  in  the  fused 
shock  region  are  based  on  a  flat  reflecting  surface  and 
variations  from  this  case  may  have  indeterminate  effects  on 
the  actual  overpressure  values.  Also,  it  is  important  to 
realize  in  evaluating  results  that  the  scaling  techniques  do 
not  hold  for  regions  in  and  on  the  boundary  of  the  fireball. 
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This  is  true  because  the  scaling  techniques  are  not  applicable 
in  regions  where  radiation  is  the  primary  means  of  energy 
transport  (Ref  7,  31). 

Accuracy  of  Blast  Calculations.  Subroutine  BLAST1 

"  "  V  ”  ' 

computes  blast  parameters  using  the  same  techniques  that  are 
used  in  the  SABER  code.  The  calculated  results  will  be  the 
same  for  the  two  programs.  A  correlation  between  SABER 
results  and  the  SAP  hydrodynamic  code  has  been  made  (Ref  7, 
42-46).  The  results  of  that  correlation  are  presented  in 
Figs  5,  6,  and  7.  The  first  two  figures  are  for  the  coalti¬ 
tude  case  and  the  last  figure  is  non-coaltitude  for  verifica¬ 
tion  of  the  correction  for  the  nonhomogeneous  atmosphere.  An 
additional  correlation  of  SABER  results  to  measured  values 
from  actual  tests  is  presented  in  the  appendix  of  Ref  7.  The 
close  correlation  of  SABER  results  to  the  extremely  accurate 
hydrodynamic  calculations  indicates  that  the  scaling  techniques 
which  are  employed  are  adequate  to  represent  the  blast 
phenomenon. 

Summary 

The  methods  of  blast  effects  computations  were  drawn 
from  AFWL-TP-70-85  (Ref  7).  The  approach  includes  an 
accounting  for  inhomogeneous  atmosphere,  fused  shock  region, 
blast  efficiency,  and  receiver  motion.  Scaling  techniques 
and  empirical  data  ferra  the  basis  for  all  the  blast  parameter 
calculations . 
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Fig.  5.  Comparison  of  SABER  and  SAP  Results  for  300  KT 
at  65,000  feet  (Ref  7,  43). 
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Fig.  6.  Comparison  of  SABER  and  SAP  Results  for  300  KT 

at  131,240  feet. 
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Fig,  7.  Comparison  of  SABER  and  SAP  Results  for  6  MT 

at  60,000  feet  (Ref  7,  46). 
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A  subroutine  named  BLAST1  was  developed  to  make  the 
effects  calculations  and  to  perform  the  system  survivability 
analysis.  Survivability  is  based  on  the  threat  levels  and 
vulnerabilities  to  peak  overpressure  and  peak  dynamic 
pressure.  The  subroutine  is  capable  of  evaluating  blast 
effects  up  to  250,000  ft. 
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IV.  Computation  of  Thermal  Environments 
Thermal  Phenomenon 

The  investigation  of  thermal  effects  of  nuclear  weapons 
on  systems  depends  on  three  basic  factors.  They  are  the 
source  of  thermal  energy,  the  transmission  of  that  energy 
through  a  medium,  and  the  absorptivity  of  a  receiver  in  that 
medium.  The  scope  of  this  study  is  limited  to  the  considera¬ 
tion  of  the  first  two  factors. 

Analysis  of  the  thermal  environment  is  at  best  an 
approximate  effort.  This  is  so  because  of  the  variability 
from  day  to  day  and  hour  to  hour  of  atmospheric  properties 
which  affect  the  energy  transport.  Consequently,  there  is 
no  single  approach  that  can  account  for  all  the  factors  and 
geophysical  conditions.  Simplified  model  environments  have 
been  developed  from  which  results  should  be  carefully 
evaluated.  The  method  presented  in  this  report  is  such  a 
model . 

The  first  topic  treated  here  describes  the  determination 
of  the  thermal  yield.  The  remaining  paragraphs  describe  the 
methods  used  to  determine  thermal  fluence  levels  depending 
on  the  location  of  both  the  burst  and  the  receiver.  A  method 
is  presented  to  determine  thermal  levels  of  both  reflected 
and  direct  energy  for  burst  altitudes  less  than  30  kft.  with 
receivers  below  100  kft.  There  is  no  change  in  the  approach 
for  bursts  between  30  and  50  kft.  except  that  reflected 
energy  is  negligible  in  this  region.  Procedures  and  relations 
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have  been  included  which  are  applicable  to  the  case  of  a 
burst  above  50  kft.  with  the  receiver  below  50  kft.  Finally, 
the  case  for  both  burst  and  receiver  above  50  kft.  is  treated. 

Thermal  Yield 

The  thermal  yield  is  a  function  of  the  altitude  of  burst, 
and  the  total  yield  of  the  weapon.  In  general,  the  thermal 
yield  is  a  fraction  of  the  total  yield  and  this  fraction 
rises  with  increasing  altitude  and  then  drops  to  nearly  zero 
around  350  kft.  (100  km).  The  low  and  intermediate  altitude 
burst  (<  100  kft.)  is  characterized  by  two  distinct  thermal 
pulses,  the  second  of  which  contains  the  bulk  of  the  thermal 
energy.  At  higher  altitudes  the  first,  and  second  pulse 
combine  to  form,  a  single  pulse  of  very  high  intensity  and 
short  duration. 

Data  from  the  SPUTTER  AFWL  weapons  effects  code  has 
provided  the  thermal  efficiencies  for  a  range  of  burst  yields 
and  altitudes.  An  empirical  relation,  a  function  of  altitude 
alone,  has  been  obtained  from  the  collapsed  sputter  data  and 
is  plotted  in  Fig.  8.  Relations  for  thermal  efficiency  and 
thermal  yield  were  furnished  from  unpublished  information, 
courtesy  of  Mr.  A.  Sharp  of  AFWL,  Kirtland  AFB,  New  Mexico. 
The  thermal  efficiency,  T,  is  given  by 

T  o  e(-.3579  -  . 008H  +  7.136  x  10“4H2)2.30  ^ 

(-1.25  x  10'V  +  6.42  x  10‘8H4)2.30  ..... 

e  (21) 
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where 

a  =  height  of  burst  (Km) 

The  thermal  yield  (air  burst)  is  given  by 

W  h  =  TW*94  (22) 

where 

W  =  total  yield  in  KT 
The  thermal  yield  (ground  burst)  is 

1  03 

W  h  =  -35W (23) 

The  time  of  the  second  thermal  peak  (maximum  power 
output)  is  taken  as  the  time  at  which  all  the  thermal  energy 
is  released.  An  empirical  relation  for  that  time  is  given 
by  (Ref  2,  175) 

T(sec)  =  . 0491 (Wn) *  42  (24) 


where 

n  =  ratio  of  air  density  at  burst  altitude  to  that  at 
sea  level 

Atmospheric  Transmission  (Burst  <  50  kft.) 

A  model  for  transmission  of  thermal  radiation  in  the 
atmosphere  has  been  developed  in  Ref  3.  That  model  is 
presented  below  and  is  followed  exclusively  for  thermal 
calculations  except  as  noted  by  individual  references. 

Six  basic  assumptions  are: 

1.  The  fireball  of  an  air  burst  is  essentially  spherical 
and  is  treated  as  a  point  source  black  body  radiator 
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at  6000  degrees  Rankine.  The  fireball  for  the 

r 

'  surface  burst  is  treated  as  a  hemispherical 

Lambert-type  emitter  of  uniform  temperature.  The 
black  body  temperature  is  taken  as  3000  degrees 
Rankine. 

2.  The  transmission  of  visible  energy  is  evaluated 
for  a  collimated  beam  of  energy.  This  assumption 
is  made  for  the  purpose  of  describing  the  geometric 
path  of  reflected  energy  and  to  allow  application 
of  Lambert’s  Law. 

3.  The  visible  energy  is  assumed  to  pass  through  a  non¬ 
absorbing  but  scattering  media. 

4.  Attenuation  of  infrared  energy  is  assumed  to  be  the 
result  of  absorption  by  the  water  vapor  in  the 
atmosphere. 

5.  The  terrain  surface  .s  taken  as  a  Lambert  plane, 
infinite  and  flat. 

6.  The  albedo,  or  surface  , cflecticn  factor,  is  known 
and  is  the  sam:  for  all  frequencies  of  the  reflected 
energy  (see  Table  I  on  the  next  page). 

An  application  of  Lambert’s  Law  was  made  to  yield  an 
analytic  expression  for  the  ieflected  thermal  energy. 
Lambert's  Law  expresses  the  intensity  of  reflected  radiant 
energy  as  a  function  of  the  angle  from  which  it  is  viewed 
(Ref  1,  111-112).  The  surface  of  reflection  is  taken  as  a 
black  body  radiator  and  as  a  perfectly  scattering  surface. 

4  V 
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Table  I 

Values  of  Albedo  for  Earth  Surfaces 


Surface 

Reflectance  (per  cent) 

Desert 

24-28 

Fields  (Various  Types) 

3-40 

Forest  (Green) 

3-20 

Grass  (Various  Conditions) 

14-47 

Ground  (Bare) 

7-30 

Snow 

65-89 

Water  (Average) 

10 

Rough  Water  (White  Caps) 

10-31 

Shock-Frothed  Water 

40-80 

Sand  (Dry) 

18-28 

Sand  (Wet) 

9-19 

(Ref  3,  54} 

Reflected  Radiation  (Air  Burst) .  The  thermal  fluence  of 
an  air  burst  reflected  to  a  plane  horizontal  receiver  from  an 
infinite  Lambert  plane  is  given  by 

w 

Q«  ”  (^i?)TtptY(HRo'XRo’l  (25) 

where 

=  reflected  thermal  fluence  (cal/cm^) 

S  e  slant  range,  ground  zero  to  receiver  (cm) 
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T  *  atmospheric  transmission  fraction 
p  =  albedo 

^^Ro'*Ro^  =  Samma  function  for  air  burst,  (a  function 
of  scaled  height  of  the  receiver  and  the 
scaled  x-coordinate  range  of  burst  to 
receiver) 

The  gamma  function  has  been  evaluated  by  numerical  techniques 
for  a  range  of  values  of  H_  and  X_  .  Both  Hn  and  X„  are 
scaled  by  the  height  of  the  burst  above  the  ground  level. 

The  gamma  function  described  above  is  not  to  be  confused  with 
the  well  known  y  function  of  mathematics. 

Reflected  Radiation  (Ground  Burst) .  The  thermal  f luence 
of  a  ground  burst  reflected  to  a  plane  horizontal  receiver 
from  an  infinite  Lambert  plane  has  been  developed  in  Ref  4 
and  is  given  by 


w 

°R”  ( 


(26) 


where 

Y(So,4>)  =  gamma  function  for  ground  burst,  (a  function 
of  the  scaled  slant  range  and  the  angle  <}>) . 

The  gamma  function  has  been  evaluated  numerically  for  ranges 
of  the  two  variables  and  <p  (Ref  4,  288),  The  slant  range, 
S,  is  scaled  to  the’ radius  of  the  fireball  and  4>  is  the  angle 
between  the  vertical  at  the  receiver  and  the  slant  range  to 
ground  zero.  The  thermal  fluence  obtained  for  both  the  air 
and  the  surface  burst  is  that  for  normal  incidence  on  a  plane 
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horizontal  receiver.  Equations  (25)  and  (26)  are  applicable 

for  burst  heights  up  to  and  including  30  kft. 

Direct  Radiation .  The  thermal  fluence  directly  incident 

on  a  plane  horizontal  receiver  is  attenuated  both  by  radial 

divergence  and  by  the  transmissivity  of  the  intervening 

2 

atmosphere.  The  incident  thermal  fluence  (cal/cm  )  for  an 
air  burst  is  given  by 


QD  =  — ^  <V 

U  4 it  (SR)  Z  U 


where 

Tp  »  direct  thermal  energy  transmittance 
SR  =  slant  range,  burst  to  receiver  (cm) 

Wth  »  thermal  yield  (KT) 

2 

The  incident  thermal  fluence  (cal/cm  )  from  a  surface  burst 
is  given  by 


2 7T  (SR) 


i  <V 


The  normally  incident  fluence  may  be  found  by  multiplying 
equations  (27)  and  (28)  by  the  appropriate  trigonometric 
relations . 

Transmissivity  Fractions .  Thermal  fluence  for  bursts 
below  50  kft.  is  composed  primarily  of  visible  and  infrared 
energy.  Transmission  of  the  two  energy  types  is  treated 
separately  since  the  visible  energy  is  primarily  subject  to 
scattering  and  the  infrared  energy  is  primarily  subject  to 
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absorption  by  water  vapor.  Both  forms  of  radiation  are 
affected  by  the  amount  of  haze,  described  as  a  visibility 
factor,  that  is  present  in  the  environment. 

The  following  equations  give  the  transmissivity  func¬ 
tions  for  visible  light  in  air  free  of  haze  including  both 
coaltitude  and  non-coaltitude  cases  (Ref  4,  28-29). 


cos  Ip  = 


(Z  -  Z  ) 
r  s' 


Z  t  1 
s  r 


T  *  e 
r 


-.0875 
cos  ip 


-4.57  x  10'3(Z  ) 


-4.57  x  lO^CZ  )  . 

.  r  )) 


z  «  z 
s  r 


T  o  e 
r 


[-«  x  10-6(SR)e-4-57  *  10  <V] 


where 


*  transmissivity  factor  (visible) 

Z^  =  height  of  receiver 

Z  «*  height  of  source 
s 

The  transmissivity  factors  for  visible  light  in  air  with 
haze  present  are  given  by 


Z  /  Z 
s  r 


Th  ■ e 


-16.4 
V  cos  ip 


[e”4 • 57  x  10”b(Zs) 


e-4.57  x  10  (Zj.)  j 
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i  f 

V/ 


z  =  z 

s  r 


-7.5  x  10‘4(SR)  r_-4.57  x  10"5(Z_)1 

•  v  l  ®  *  J 

Th  =  e  '  (33) 


where 


=  transmissivity  factor  for  visible  (with  haze) 
V  »  visibility  at  sea  level  (miles) 


fc:* 


The  transmissivity  factor  for  infrared  energy  is 
dependent  on  the  amount  of  water  vapor  present  (given  in 
precipitable  millimeters  of  water  along  the  transmission 
path).  The  water  vapor  present  in  the  path  is  given  by 


Zs  *  Zr 


2.3  P 
_ c 

COS  l|» 


-6.1  x  10 


[10 


<V 


-6.1  x  10 


-5 


-  10 


(,r)] 


Z  »  Z 
s  r 


.  -6.1  x  10"5(Z  ) 

w  a  3.23  x  10"  (P  )(SR)[10  S  ) 


where 

w  =  precipitable  millimeters  of  water 
Pq  o  water  vapor  pressure  at  sea  level  (mmHg) 
Finally  the  total  transmissivity,  T^,  is  given  by 


(34) 


(35) 


T*  =  F  TdT,  +  F.  T  (.7  +  .  3T,  ) 
t  v  R  h  lr  w^  h' 


(36) 
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where 

Fy  *  fraction  of  thermal  fluence  in  visible  region 
F^r  =  fraction  in  infrared  region 

T  *  infrared  transmission  factor 

w 

The  fractions  of  visible  and  infrared  energies  are  functions 

of  the  black  body  temperature.  The  energy  is  split  evenly 

between  visible  and  infrared  for  the  6000  degree  air  burst. 

For  the  ground  burst  at  3000  degrees,  F  =  .1  and  F.  -  .9 

v  lr 

(Ref  3,  35). 

The  infrared  transmission  factor  is  given  as  a  function 
of  the  precipitable  water  vapor  in  the  transmission  path. 
Curves  for  both  the  air  and  the  ground  burst  are  given  in 
Fig.  9. 

Combinations  of  the  basic  equations  presented  in  this 
section  will  allow  transmission  fractions  to  be  compiled 
for  reflected  energy  and  for  situations  where  there  is  a 
haze  layer  from  ground  level  up  to  a  specified  altitude. 

The  interested  reader  is  referred  to  Ref  3,  pp.  56-60, 

The  High  Altitude  Burst 

Transmission  of  thermal  fluence  from  a  burst  above 
50  kft.  to  a  receiver  below  50  kft.  has  been  developed  by  a 
separate  approach.  The  ultraviolet  fraction  of  the  thermal 
energy  is  a  significant  portion  of  the  thermal  energy  for 
bursts  in  this  region.  The  thermal  fluence  incident  on  a 
plane  horizontal  receiver  is  given  by 
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Q  .  _ -  [F  e'1-31  *  10'4'R)  . 

(Z  -  Z  )2  uv 
v  s  r' 


,,  _  .-1.31  x  lO”5 (R) ,  .2  . 

(1  -  F  )e  ]sm  <f> 


uv 


(37) 


where 


sin  (J)  *  — — 


(Z  -  Z  ) 

'  e  y*  / 


SR 


(38) 


F  =  fraction  of  ultraviolet  energy 
uv 

and  for  Z  <  28,000  ft.  the  reduced  range,  R,  is 
r 


31746 
sin  <f> 

9987 
sin  <p 


[e 


[1 


•3. IS  x  10~5(Z  ) 


.-5 


-  .414]  + 


-4,55  x  10  ''(Z  -  28,000) 

e  S  ]  (39) 


and  for  28,000  <  Z  <  50,000 

r 


9987 
sin  4> 


(1  -  e 


4.55  x  1 0~ 5 (Z  -  28,000) 


3  (40) 


The  High  Altitude  Burst  and  Receiver 

The  thermal  transmission  factor  is  essentially  unity 

when  both  the  burst  and  receiver  are  above  50  kft.  Radial 

divergence  is  the  only  attenuation  present  in  this  region. 

2 

The  thermal  fluence  in  cal/cm  is  given  by 
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Q  -  - ^-2  (41) 

4ir(SR)  ^ 

Subroutine  for  Thermal  Effects 

Subroutine  THERM,  The  subroutine  for  calculation  of 
thermal  effects  and  survivability  to  those  effects  is  called 
THERM.  This  subroutine  is  called  by  the  program  GUIDE  and 
returns  to  GUIDE  a  record  of  those  vehicles  that  have  been 
destroyed.  Thermal  energy  levels  are  computed  at  each 
vehicle  position.  Survival  of  each  vehicle  5s  based  on  a 
comparison  of  the  vulnerability  level  and  the  free  field 
thermal  energy  level. 

Reflected  energy  from  the  surface  of  the  earth  may  be 
significant  under  certain  conditions.  A  survivability  test 
is  also  performed  based  on  the  sum  of  the  reflected  and  the 
direct  energy  normal  to  a  plane  horizontal  receiver.  This 
test  is  made  only  when  the  reflected  and  direct  energy  is 
additive,  i.e.,  for  a  burst  below  the  receiver. 

Parameters  Compute-  .  The  results  of  the  thermal 
computations  for  each  burst  are  printed  from  the  THERM 
routine.  All  vehicle  losses  are  listed  in  the  output  under, 
"Results  from  Thermal  Effects  Computations".  The  lethal 
level  is  also  included  as  printed  information.  The  following 
additional  information  is  printed  if  the  option  for  data  is 
chosen: 
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2 

Unattenuated  free  field  fluence  (cal/cm  ) 

Height  of  burst  (ft) 

Slant  range  at  time  of  peak  radiant  power  (mi) 

Time  to  peak  radiant  power  (sec) 

2 

Direct  free  field  fluence  (cal/cm  ) 

2 

Direct  normally  incident  fluence  (cal/cm  ) 

(Normal  to  a  plane  horizontti  receiver) 

2 

Reflected  fluence  (cal/cm  ) 

2 

Reflected  plus  direct  normal  fluence  (cal/cm  ) 

(Computed  only  if  burst  below  the  receiver) 

Sub  Programs  Called  by  THERM.  Subroutines  ATMOS,  TRANS, 
SETUP,  and  MACUP.E  are  called  from  the  routine  THERM.  Sub¬ 
routine  TRANS  is  called  upon  to  evaluate  the  atmospheric 
transmission  factors  for  thermal  energy  as  needed  in  the 
thermal  calculations.  The  effect  of  a  haze  layer  from  sea 
level  to  any  specified  altitude  is  included  in  the  calcula¬ 
tion  of  the  transmission  factors.  The  remaining  subroutines 
perforin  the  same  functions  as  described  in  Chapter  III. 

BLOCK  DATA  is  again  used  for  storage  of  data  tables. 

Limitations  on  Thermal  Calculations .  Methods  have  been 
presented  for  determining  transmission  of  thermal  energy  in 
all  altitude  ranges.  However,  the  determination  of  thermal 
efficiency  is  accurate  only  to  approximately  250,000  ft. 
Extrapolation  beyond  this  limit  is  possible  but  no  verifica¬ 
tion  of  accuracy  has  been  made. 

Calculation  cf  thermal  effects  for  burst  altitudes  below 
50  kft.  include  the  effect  of  visibility.  The  method  for 
finding  transmission  fractions  is  accurate  for  transmission 
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distances  up  to  and  including  the  visibility  range.  Thermal 
levels  that  are  computed  for  distances  exceeding  that  range 
are  likely  to  be  higher  than  actual  values.  Also,  the 
attenuating  atmosphere  is  considered  to  extend  throughout 
the  transmission  path  for  the  case  of  a  burst  below  50  kft. 
with  the  receiver  above  100  kft.  This  assumption  will  result 
in  some  excess  attenuation  of  the  thermal  energy. 

The  presence  of  clouds  and  stratified  layers  of  haze 
will  greatly  affect  the  thermal  energy  present  at  any  given 
point.  These  effects  have  not  been  accounted  for.  Devia¬ 
tions  from  the  assumptions  of  this  simplified  model  will 
nearly  always  be  present  and  the  results  should  be  inter¬ 
preted  accordingly. 

Evaluation  of  Thermal  Computations .  Thermal  levels 
computed  by  subroutine  THERM  have  been  compared  to  values 
computed  by  the  thermal  code  SNAPT  at  AFWL.  The  codes  are 

based  on  a  similar  theoretical  approach.  A  correlation  of 

( 

free  field  thermal  levels  was  made  for  a  range  of  burst 
altitudes  from  0  to  60,000  ft.  and  for  receivers  from  0  to 
50,000  ft.  Agreement  to  within  10%  of  the  SNAPT  values  was 
obtained  except  for  the  ground  burst  case  where  the  THERM 
values  were  a  factor  of  two  times  greater  than  those  of  SNAPT. 
A  different  approach  was  used  in  the  two  codes  for  the 
determination  of  thermal  yield  for  the  ground  burst.  The 
SNAPT  approach  is  intended  to  account  for  a  high  attenuation 
factor  near  the  earth's  surface  and  is  probably  more  correct 
in  that  region. 
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“ihe  method  for  prediction  of  thermal  fluence  levels  in 
the  lower  altitude  regions  (below  50  kft.)  has  been  corre¬ 
lated  to  actual  tests.  There  is  general  agreement  to  within 
+  10%  for  those  conditions  meeting  the  assumptions  presented 
in  this  report.  The  calculations  for  a  high  altitude  burst 
with  a  receiver  in  the  attenuating  atmosphere  (below  100  kft.) 
are  estimated  to  be  correct  to  within  a  factor  of  two. 

Thermal  transmission  for  the  high  altitude  burst  and  receiver 
case  is  not  affected  by  atmospheric  attenuation.  Radial 
divergence  is  applied  to  find  the  thermal  level  and  the 
accuracy  is  only  limited  by  the  prediction  of  thermal  yield. 

Summary 

The  methods  of  thermal  effects  computations  were  drawn 
primarily  from  Ref  3.  The  transmission  fractions  were  found 
by  separate  techniques  for  three  regions  of  altitude.  The 
regions  handled  separately  were  the  following:  burst  alti¬ 
tudes  less  than  30  kft.  with  receivers  below  100  kft.,  burst 
altitudes  above  50  kft.  with  receivers  below  50  kft.,  and 
both  burst  and  receiver  above  50  kft. 

A  subroutine  named  THERM  was  developed  to  make  the 
effects  calculations  and  to  perform  the  survivability 
analysis.  Survivability  is  based  on  the  threat  level  and 
vulnerability  for  the  free-field  thermal  fluence  in  cal/cm  . 
The  subroutine  is  capable  of  evaluating  thermal  effects  for 
bursts  up  to  approximately  250,000  ft. 
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V.  Results  and  Conclusions 

A  computer  code  has  been  developed  that  will  provide  a 
method  to  make  survivability  studies  on  systems  that  axe 
undergoing  a  nuclear  attack.  Initially,  subroutines  have 
been  included  to  evaluate  survivability  of  systems  to  the 
effects  of  blast  and  of  thermal  energy.  The  code  is  designed 
to  readily  accept  similar  subroutines  for  the  evaluation  of 
other  nuclear  burst  effects. 

The  program  was  prepared  for  use  on  the  CDC  6600 
computer  with  a  Scope  3.3  compiler  version.  The  code  is 
written  in  the  FORTRAN  EXTENDED  language.  Core  memory 
required  on  the  CDC  system  is  approximately  40,000  octal 
words  and  run  times  are  on  the  order  of  a  few  seconds. 

The  code  is  capable  of  handling  from  one  to  ten 
sequential  bursts  which  are  automatically  targeted  and, of 
evaluating  nuclear  effects  on  a  maximum  of  100  vehicles. 

The  number  of  bursts  which  may  be  entered  with  positions 
specified  by  the  user  is  not  limited.  Any  type  of  system 
may  be  studied  for  which  the  vulnerability  limits  are  known. 

Subroutines  for  the  separate  effects  each  have  separate 
limitations.  The  user  of  this  program  should  be  familiar 
with  those  limitations  to  obta>  i  the  most  accurate  results. 
The  results  from  blast  computations  have  been  correlated  to 
actual  tests  and  to  hydrodynamic  calculations.  A  close 
agreement  has  been  verified  for  all  regions  up  to  altitudes 
of  150,000  feet.  The  accuracy  of  calculations  for  thermal 
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effects  varies  depending  on  the  altitude  of  the  burst  and 
receiver  combination.  Some  correlation  to  actual  test  data 
has  been  made  and  a  close  agreement  was  found  between  this 
thermal  routine  and  the  SNAPT,  AFWL  routine.  The  values 
for  free  field  blast  and  thermal  effects  computed  by  this 
code  should  be  considered  adequate  for  systems  analysis  and 
survivability  considerations.  A  sample  problem  which  illus¬ 
trates  the  use  of  this  code  has  been  included  as  Appendix  E. 

The  presence  of  clouds,  temperature  inversions  and  other 
atmospheric  variations  can  have  a  significant  effect  on  both 
blast  and  thermal  levels.  Also  deviations  from  the  assumed 
flat  reflecting  surface  of  the  earth  are  important  for  the 
lower  altitudes.  These  factors  should  always  be  recognized 
in  conjunction  with  an  interpretation  of  results. 

Some  areas  of  this  code  that  could  be  further  investi¬ 
gated  or  improved  upon  include  the  following  items.  The 
capabilities  of  the  thermal  effects  evaluation  could  be 
improved  by  accounting  for  the  effects  of  reflected  thermal 
energy  from  cloud  layers  above  the  point  of  interest.  The 
additive  effects  of  free  field  levels  for  coincident  bursts 
has  not  been  treated  and  could  be  a  significant  factor  for 
closely  spaced  bursts.  No  provision  has  been  made  to 
evaluate  the  environment  of  special  weapons.  This  capability 
could  be  very  important  in  survivability  evaluations  and  is 
a  desirable  addition  to  the  code.  The  development  of  sub¬ 
routines  to  calculate  environments  for  x-ray,  gamma  ray, 
neutrons,  and  Electro-Magnetic  Pulse  should  also  be 
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accomplished.  This  addition  to  the  basic  survivability 
code,  along  with  the  blast  and  thermal  subroutines,  would 
then  provide  a  single  package  to  evaluate  survivability  to 
each  of  the  major  nuclear  effects. 
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Appendix  A 

A  Statistical  Approximation  for  a  Normally 
Distributed  Variable 


The  Central  Limit  Theorem  of  probability  theory  was 
used  to  obtain  a  normal,  spherical  distribution  for  the 
detonation  point  relative  to  the  chosen  target.  This  theorem 
states  (Ref  6,  251)  that  if  x^,  is  a  random  variable  taken 
from  sample  size  n,  then 


l  -  E[Ixi] 

i  =  l  1 _ l_ 

Sri  a 

x 


(42) 


where 


z 


n 


E(xi).= 
EClx.)  = 

OX  = 


value  of  a  random  variable  which  is  approximately 

normally  distributed  with  a  mean  of  zero  and  a 

variance  of  one 

expectation  value  of  x^ 

summation  over  n,  of  expectation  values 

standard  deviation  for  x. 

i 


The  coordinates  of  a  normal,  spherical  distribution  may 
be  found  in  the  following  manner: 

For  a  value  of  n  =  12  and  for  random  numbers  between  0 
and  1  we  have 


E[x.]  »  1/2 


(43) 
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Fig.  10.  Schematic  Representation  of  Random  Number 

Distribution. 


EfXxj]  =  12(1/2)  =  6 


Variance  =  cr^  =  =  1/12 


Std .  Dev.  =  ax  =  1//T2 


then 


/n  ax  =  fl2  •  1//12  =  l 


and  the  approximately  normally  distributed  random  variable  is 


given  by 


z*  * 


12 

I  x.  -  6 


This  value,  z^  ,  has  been  determined  for  a  distribution 
of  variance  of  one.  A  specified  spherical  probable  error, 
spe,  then  delineates  the  distribution.  The  value,  z^  ,  must 
be  scaled  by  the  appropriate  standard  deviation.  The 
spherical  distribution  has 
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0  =  C  =0  -0 

x  y  2 


spe  =  1.538a  (Ref  5,  301) 


therefore 


<J  «  _g 

1.538 


Finally,  for  a  given  distribution.  A,  an  approximately 
normally  distributed  point,  (x,y,z)  is  given  by 


(x,y,z)  «  aA(z(1),  zC2),  z(3))  (52) 
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Appendix  B 


Index  of  Variables  (Major  Routines) 


Note:  Variables  are  defined  in  the  first  routine  in  which 
they  appear  and  are  not  listed  again. 


Program  GUIDF: 

Variable 

Name 


BURSTX 

BURSTY 

BURSTZ 

CYCLE 

COUNT 


DISPX 

DISPY 

DISPZ 

DISPR 


HG 

HB 

HZ 

II 

IODAT 

IRC 

KTEMP 

NUM 

0PT1 

0PT2 


Definition 


First 

Appearance 
(Line) _ 


Coordinate  points  of  the  burst 
position 

Number  of  bursts  to  be  analyzed 

Record  of  burst  number  currently 
under  analysis 


XfY,Z  displacements,  burst  to 
receiver  position 

Radial  displacement,  burst  to 
receiver 

Height  of  ground  above  sea  level 

Height  of  burst  above  sea  level 

Height  of  receiver  above  sea  level 

Number  of  vehicles  input 

Input  parameter  for  data  print 
option 

Integer  number  to  precycle  the 
random  number  generator 

Input  parameter  for  non-standard 
temperature  option 

Array  of  possible  values  for  IRC 

Option  for  method  of  burst  placement 

Option  to  change  velocity  vectors 


13 


13 

125 

126 
127 


128 

8 

8 

8 

6 

6 

102 

f> 

12 

13 

13 
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GUIDE  (cont 

.) 

Name 

Definition 

Line 

PO 

Water  vapor  pressure  at  sea  level 

11 

POSX 

10 

POSY 

X,Y,Z  arrays  for  vehicle  positions 

10 

POSZ 

6 

RFB 

Radius  of  fireball 

70 

RHO 

Albedo 

11 

SPE 

Spherical  probable  error  for 
targeting 

89 

TARGX 

92 

TARGY 

X,Y,Z  coordinate  points  on  which  burst 

93 

TARGZ 

is  targeted 

94 

TEMPS 

Nonstandard  sea  level  temperature 
(degree  Rankine) 

6 

TIMEX 

Incremental  time  for  each  burst 

48 

TM(1 ) 

Temperature  at  receiver  altitude 
(degree  Rankine) 

6 

TM  (2) 

Temperature  at  ground  altitude 
(degree  Rankine) 

6 

TM  (3) 

Temperature  at  burst  altitude 
(degree  Rankine) 

6 

V BLAST I 

Vulnerability  level  for  overpressure 

6 

VEHIC 

Array  for  record  of  vehicles 

6 

VGAMA 

Vulnerability  level  for  gamma 

45 

VIS 

Visibility  at  sea  level  (miles) 

11 

VNEUT 

Vulnerability  level  for  neutrons 

45 

VQ 

Vulnerability  level  for  dynamic 
pressure 

45 

VTHERM 

Vulnerability  level  for  thermal 
energy 

45 

VXRAY 

Vulnerability  level  for  x-ray 

45 
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GUIDE :  (cont.) 


Name 

Definition 

Line 

VX 

VY 

X,Y,Z  velocity  component  vectors  for 

VZ 

each  vehicle 

10 

w 

Yield  in  KT 

8 

ZH 

Height  of  haze  layer  above  sea 
level  (ft.) 

11 

Subroutine 

TARGET: 

DUM 

Dummy  variable  for  argument  of  random 
number  function 

9 

R 

Radial  displacement  of  burst  from 
target  position 

29 

RX1 

26 

RY1 

X,Y,Z  coordinate  displacement  of 

27 

RZ1 

burst  from  target  position 

28 

SIGMA 

Standard  deviation  for  a  normal 
distribution 

7 

SUMX 

3 

SUMY 

Computational  sum: 

4 

SUMZ 

5 

X 

Y 

Z 

X,Y,Z  coordinates  of  the  target  point 

1 

Subroutine 

BLAST1 : 

ACF 

Table  of  altitudes  for  blast 
efficiency 

21 

ALFA 

Ledsham-Pike  a  correction  factor 

132 

ALP1E 

Minimum  angle  at  which  mach  reflection 
will  occur  (degrees) 

13 

ALP1ER 

ALP1E  (radians) 

13 

CF 

Table  of  altitudes  corresponding  to 
table  of  efficiency  factors,  CFF 
(feet) 

21 
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Subroutine 

BI.AST1:  fcont.) 

Name 

Definition 

Line 

H(l) 

H(2) 

HZ,  KG,  HB  respectively 

27 

H(3) 

HT 

Height  above  ground  of  point  on 
triple-point  path  (kft) 

24 

INCOMP 

64 

INTC 

Variables  controlling  program  flow 

90 

KCASE 

72 

KER 

Variable  indicating  error  in  ATMOS 

100 

LER 

Dummy  argument  in  call  to  MACURE 

77 

NIT 

Variable  for  program  flow 

59 

PCD 

P(2) 

PZ,  PB,  PG  respectively  (set  in 

P(3) 

equivalence  statement) 

27 

PB 

Pressure  at  burst  altitude  (psi) 

13 

PBPZ 

Ratio  PB/PZ 

13 

PBR 

Ratio  PB/PSL 

13 

PBRW 

(PBR/W)**.  3333 

13 

PDMV 

Positive  duration  of  material 
velocity  (sec) 

169 

PDOD 

Positive  duration  of  overdensity 
(sec) 

171 

PDOOP 

Positive  duration  of  overpressure 
(sec) 

167 

PG 

Pressure  of  ground  elevation  (psi) 

13 

PGR 

Ratio  PG/PSL 

13 

PHIR 

Angle  used  in  triple  point 
consideration 

13 

POD 

3 

Positive  overdensity  (slugs/ft  ) 

170 

?R(i) 

Ratio  P(1)/PSL)  (PZR) 

27 
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Subroutine 

BLAST1:  (cont.) 

Name 

Definition 

Line 

PR(2) 

Ratio  P (2) /PSL  (PGR) 

27 

PR  (3) 

Ratio  P (3)/PSL  (PBR) 

27 

PSL 

Atmospheric  pressure  at  sea  level 
(psi) 

58 

PZ 

Atmospheric  pressure  at  receiver 
(psi) 

13 

PZR 

Ratio  PZ/PSL 

13 

R 

Radius  of  the  fireball  (kft) 

13 

RA 

Slant  range  from  burst  to  point  on 
triple-point  path  (kft) 

13 

RBAR 

Scaled  range  (1  KT,  sea  level  burst) 
(kft) 

124 

RH0(1) 

RHO (2) 

RHO (3) 

Density  at  HZ,  HG,  HB  respectively 
(slugs/ft3) 

100 

RHOB 

RHOG 

RHOZ 

RHO (3) ,  RHO ( 2) ,  RHO(l)  respect?  v 

28 

SDELP 

Scaled  peak  overpressure  o  1  KT, 
sea  level  burst  (psi-) 

126 

SEPD 

Separation  distance  'hock  front 
to  vehicle  positio.: 

140 

SEPDT 

Separation  distance  for  time 
increment  DT 

140 

SFR 

Shock  front  range  from  burst  position 

92 

SFV 

Shock  front  velocity  (ft/sec) 

140 

SR 

Slant  range  from  burst  to  receiver 

56 

SRE 

Slant  range  from  burst  to  origin 
of  the  triple-point  path  on  ground 
(kft) 

13 
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Subroutine 

Name 

ss(i) 

SS  (2) 

SS  (3) 

SSB 

SSG 

SSZ 

SSZR 


ST 


STO 


TAB v  ,LL) 

TAB1I 

TABID 


TAB2I 

TAB2D 


TAB3I 

TAB3D 

TAB4I 

TAB4D 

TAB5I 

TABSD 


TAB6I 

TAB6D 

TDPZ 


TDRZ 


BLAST1 :  (cont.) 

Definition 


Ambient  speed  of  sound  at  receiver, 
at  ground,  and  at  burst  altitudes 
respectively  (ft/sec) 


SS(3),  SS(2),  SS(1)  respectively 


Ratio  of  speed  of  sound  at  receiver 
to  that  at  sea  level 

Horizontal  range  from  ground  zero 
to  point  on  triple-point  path  (kft) 

Horizontal  range  from  ground  zero 
to  the  beginning  of  triple-point 
path  (kft) 

Tabular  data  stored  in  BLOCK  DATA 

Values  of  scaled  range,  RBAR, 
(Independent  Variable) 

Values  of  scaled  peak  overpressure, 
SDELP  (Dependent  Variable) 

Values  of  scaled  range,  RBAR 
Values  of  Ledsham-Pike  alpha 
correction 

Used  to  invert  order  of  entry  of 
values  in  TAB1 

Values  of  angle  ALPHIR 
Values  of  angle  PHIR 

Values  of  scaled  height  of  burst 
Values  of  yield  amplification  factor 
for  Mach  stem  region 

Values  of  scaled  range,  RBAR 

Values  of  scaled  time  of  shock  arrival 

Scaled  positive  duration  of  over¬ 
pressure 

Ratio  of  scaled  positive  duration  of 
overdensity  to  the  scaled  duration 
of  overpressure 


Line 

27 

28 

114 

69 

67 

18 


154 

166 
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Subroutine  BLAST1  (cont.) 

Name  Definition  Line 

TEMP  Dummy  argument  for  ATMOS  call  100 

TMV  Ratio  of  scaled  positive  duration 

of  material  velocity  to  scaled  dur¬ 
ation  of  overpressure  161 

TR  Dummy  variable  in  ATMOS  call  100 

TSA  Time  of  shock  arrival  (sec)  139 

TSACAP  Time  of  shock  arrival  (scaled,  1  KT 

sea  level  burst)  138 

TT  Dummy  argument  in  MOTION  call  140 

HL(ALL-)  Separation  blocks  for  ranges  of 

the  curves  of  positive  phase 

duration  versus  RBAR  21 

VC  Dummy  variable  in  ATMOS  call  100 

WOR  Yield,  W  49 

XITER  25 

XKKX  Variables  controlling  program  flow  24 


Subroutine  THERM: 


AAAI  Table  of  independent  variables 

normalized  slant  range  and  angle, 

PFE  10 

AAA2  Table  of  ground  burst  gamma 

functions  (dependent  variable)  10 

BB1  Table  of  variables  -  HR0  and  XRO 

(independent)  10 

BB2  Table  of  gamma  functions,  air  burst 

(dependent)  10 

DELX  X  coordinate  range,  burst  to 

receiver  61 

DELZ  Z  coordinate  range,  burst  to 

receiver  62 
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Subroutine  THERM:  (cont.) 

Name  Definition  Line 

FUV  Fraction  of  thermal  yield  in 

ultraviolet  120 

GAM  Gamma  function  for  air  or  ground 

burst  80 

HOB  Height  of  burst  above  sea  level  20 

HRO  Scaled  height  of  the  receiver  85 

PFE  Angle  between  the  vertical  at  the 

receiver  and  the  slant  line  to 
ground  zero  78 

PI  Constant  of  multiplication  38 

PO  Water  vapor  pressure  at  sea  level 

(mmHg)  9 

PX  58 

PY  Vehicle  coordinate  positions  at  59 

PZ  time  of  the  thermal  maximum  60 

QDN  Direct  thermal  fluence  normal  to 

plane  horizontal  receiver  (cal/cm2)  96 

QDFF  Direct  thermal  fluence,  free  field 

(cal/cm2)  93 

QDFFUN  Direct  thermal  fluence,  free  field, 

unattenuated  94 

QR  Reflected  thermal  fluence  (cal/cm2)  90 

QMHR  Sum  of  reflected  and  direct  normal 

fluence  on  a  plane  receiver 

(cal/, cm2)  99 

R  Reduced  range  (high  altitude  burst)  123 

SRMJ  Slant  range  (miles)  69 

SRCM  Slant  range  (cm)  70 

SRN  Normalized  slant  range  77 

SS  Dummy  argument  in  ATMOS  call  21 
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Subroutine  THERM:  (cont.) 

Name  Definition  Line 

STHET  Angle  used  in  high  altitude  burst 

computations  121 

TBW1  Table  of  precipitable  water  vapor 

(independent  variable)  10 

TBT1  Table  for  fractional  transmission 

of  infrared  (ground  burst)  10 

TBT2  Table  for  fractional  transmission 

of  infrared  (air  buret)  10 

TBV1  Table  of  height  of  burst 

(independent  variable)  10 

TBV2  Table  of  fraction  cf  ultraviolet  energy 

(dependent  variable)  10 

TD  Direct  transmission  fraction  92 

TEFF  Thermal  yield  efficiency  factor  33 

TMAXZ  Time  of  2nd  thermal  max.  (sec)  30 

TMIN  Time  of  thermal  minimum  (sec)  29 

TR  Reflected  transmission  fraction  87 

TYPE  Variable  controlling  program  flow  86 

WMT  Yield  in  megatons  18 

WTH  Thermal  yield  (KT)  37 

XR0  Scaled  X  range  82 

Subroutine  TRANS : 

CHK  Variable  for  program  flow  1 

CTHETA  Angle  used  in  transmission  formulas  16 

FIR  Fraction  of  infrared  energy  69 

FV  Fraction  of  visible  energy  68 

LER  Error  record  75 
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Subroutine  TRANS:  (cont.) 


Name 

Definition 

Line 

TH 

Transmission  factor  with  ha*e 

19 

TR 

Transmission,  clear  air 

20 

TT 

Total  transmission  factor 

79 

TW 

Transmission  factor,  infrared 

75 

xxw 

Water  vapor  in  the  transmission  path 

40 

Subroutine 

TRIPNT: 

ALFA 

Ledsham-Pike  a  factor 

59 

ALPHA 

ALPHIR 

Angles  used  in  triple-point 

134 

ALPIF. 

computation 

107 

ALT 

Height  of  burst  above  ground  (kft) 

40 

ALTSRG 

Ratio  of  ALT  to  slant  range  from 
burst  to  a  point  on  the  ground 

89 

CAPD 

Function?  for  Newton-Raphson 

84 

CAPQ 

iteration  to  determine  ALFA 

85 

CONT 

Variable  for  program  flow 

125 

DDELP 

Derivative  in  Newcon-P.aphson 
iteration 

64 

DELPD 

Desired  overpressure  (Range 
solution  only) 

124 

DELPG 

Overpressure  on  ground 

62 

DELPR 

Overpressure  received,  used  in 
iteration  process 

31 

DPDX 

Derivatives  in  N-R  iteration 

75 

DRBAR 

76 

FR 

Yield  amplification  factor  (fused 
shock  region) 

46 

INCOMP 

Variable  recording  error  in  input 
data  (Range  solution  only) 

95 
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Subroutine 

TRIPNT:  (cont.) 

Name 

Definite  on 

Line 

PHIR 

Angle  used  in  triple-point 
calculations 

136 

PPP 

Exponent  used  in  functional  curve 

68 

SHB 

Scaled  height  of  burst 

42 

SMF 

Function  used  in  Newton-Raphson 
iteration 

83 

XI 

Incremental  value  of  shock  strength 

54 

XITER 

Variable  for  program  flow 

31 

XK 

48 

XXX 

Variables  for  program  flow 

127 

KKK 

128 
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Appendix  C 

Details  of  Data  Entry 


Input  data  cards  are  to  be  prepared  by  the  following 
instructions . 


Column 

Format 

Symbol 

Input 

Comments 

CARD  1 

1-12 

E12.5 

VIS 

- 

The  visibility  in  miles 
at  sea  level 

13-24 

E12.5 

PO 

- 

The  water  vapor  pressure 
in  mmHg  at  sea  level 

25-36 

E12.5 

RHO 

- 

The  albedo  factor  (a 

decimal  fraction) 

37-48 

E12.5 

ZH 

- 

Altitude  of  haze  layer 
(ft.  above  sea  level) 

CARD  2 

1-3 

13 

II 

Integer,  number  of 
vehicles  to  be  entered 
in  the  vehicle  array 
(right  justified) 

4-6 

13 

CYCLE 

Integer,  number  of 
weapon  bursts  (right 
justified) 

CARD  3 

1 

11 

KTEMP 

0 

Option  for  use  of 
standard  atmosphere 
(remainder  of  Card  3  is 
blank  for  this  value  of 
KTEMP) 

KTEMP 

1 

Option  for  nonstandard 
atmosphere  (complete 

Card  3  as  shown  below) 

2-7 

F6.0 

TEMPS 

- 

Desired  temperature  at 

sea  level  in  degrees 
Rankine 
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Column 

Format 

Symbol 

Input 

Comments 

CARD  3 

8-13 

F6.0 

TM(1) 

— 

Temperature  at  altitude 
of  receiver  (degrees 
Rankine) 

14-19 

F6.0 

TM(2) 

- 

Temperature  at  height  of 
ground  (degrees  Rankine) 

20-25 

F6.0 

TM  (3) 

- 

Temperature  at  height  of 
burst  (degrees  Rankine) 

CARD  4 

(Multiple 

card) 

1-12 

E 1 2 . 5 

onev  rn 

- 

X  coordinate  position 
in  feet  for  vehicle  I 

13-24 

E12.5 

POSY  (I) 

- 

Y  coordinate  position 
in  feet  for  vehicle  I 

25-36 

E12.5 

POSZ(I) 

- 

Z  coordinate  position 

in  feet  for  vehicle  I 

(Coordinate  points  for 
each  vehicle  are  rela¬ 
tive  to  an  arbitrarily 
established  origin  of 
coordinates  at  sea 
level . ) 


37-48 

E12.5 

VX(I) 

X  component  of  velocity 
in  ft/sec  for  vehicle  I 

49-60 

E12.5 

VY  (I ) 

- 

Y  component  of  velocity 
in  ft/sec  for  vehicle  I 

61-72 

E12.5 

V2(I) 

- 

Z  component  of  velocity 
in  ft/sec  for  vehicle  I 

(There  will  be  one  card 
number  4  for  each 
vehicle  entered  in  the 
vehicle  array.) 
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Column 

CARD  5 

Format 

Symbol 

Input 

Comments 

1-3' 

13 

Num(l) 

Any  integer  between  1 
and  100  of  the  readers 
choice  is  to  be  entered 
(right  justified) 

4-6 

13 

Num(2) 

- 

see  above 

7-9 

13 

Num(3) 

- 

•»  i» 

10-12 

13 

Num (4) 

- 

ft  ft 

13-15 

13 

Num(5) 

- 

tf  ft 

16-18 

13 

Num (6) 

- 

fl  ft 

19-21 

13 

Num(7) 

•• 

tf  tf 

22-24 

13 

Num (8) 

- 

ft  tl 

25-27 

13 

Num  (9) 

- 

ft  ft 

28-30 

13 

Num(10) 

- 

91  ft 

CARD  6 

1-12 

EJ2.5 

VGAMA 

- 

Vulnerability  level  for 
gamma  rays  (total  gamma 
f luence) 

13-24 

E12.5 

VXRAY 

- 

Vulnerability  level  for 
x-rays  in  cal/cm2 

25-36 

E12.5 

VTHERM 

** 

Vulnerability  level 
for  thermal  energy  in 
cal/cm2 

37-48 

E12.S 

VBLAST1 

- 

Vulnerability  level  for 
overpressure  in  psi 

49-60 

E12.5 

VNEUT 

* 

Vulnerability  level  for 
neutrons  (total  neutron 
f luence) 
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Column 

Format 

Symbol 

Input 

Comments 

CARD  6 

61-72 

E12.5 

VQ 

- 

Vulnerability  level  for 
dynamic  pressure  in  psi 

(Entry  of  a  zero  will 
cause  the  associated 
effect  to  be  deleted 
from  the  analysis.) 

Note: 

The  appropriate  sequence  of 
repeated  for  each  burst  that 
program. 

cards  7,  8,  9,  10  must  be 
is  to  be  treated  in  the 

CARD  7 

1-12 

E12.5 

TIMEX 

- 

Incremental  time  in  sec 
at  which  a  burst  occurs. 

(Time  begins  at  time 
equal  zero  and  TIMEX 
may  have  any  value 
including  zero.  Bursts 
separated  by  a  zero  time 
increment  are  treated  as 
separate  bursts.) 

13-24 

E12.5 

W 

- 

Weapon  yield  in  KT 

25-36 

E12.5 

HG 

- 

Height  of  ground  in  ft. 
above  sea  level 

37-38 

12 

0PT1 

0 

Option  for  user  to  enter 
the  coordinate  position 
of  burst  (right  justified) 

OPT1 

1 

Option  for  automatic 
placement  of  the  burst 
(right  justified) 

39-40 

12 

IODAT 

0 

Option  for  printout  of 
vehicle  losses  only 
(right  justified) 

IODAT 

1 

Option  for  printout  of 

all  effects  parameters 
(right  justified) 
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volumn 

Format 

Symbol 

Input 

Comments 

CARD  7 

41-42 

12 

OPT2 

0 

Crtion  to  allow  change 
in  velocity  vector 
(allowed  only  if 

OPT]  =  0) 

OPT2 

1 

Option  to  retain 
original  velocity  vectors 

CARD  8 

(required 

only  if  OPTi 

=  0) 

1-12 

E12.5 

BURSTX 

- 

X  coordinate  position 
in  feet  for  burst 

13-24 

E12.5 

BURSTY 

- 

Y  coordinate  position 
in  feet  for  burst 

25-36 

E12.5 

BURSTZ 

- 

Z  coordinate  position 
in  feet  for  burst 

(Coordinates  relative  to 
the  common  sea  level 
origin  of  coordinates 
described  above.) 

CARD  9 

(required 

only  if  OPTI 

=  0  and 

0PT2  =  0)  (multiple  card) 

1-12 

E12.5 

VX(I) 

New  X  component  of 
velocity  for  vehicle  I 
(ft/sec) 

3  i  2  4 

El  2 . 5 

VY  ( I ) 

- 

Y  component  of  velocity 
for  vehicle  I  (f  /sec) 

25-36 

E12.5 

VZ(I) 

- 

Z  component  of  velocity 
for  vehicle  I  (ft/sec) 

(There  will  be  one 
card  9  entered  for  each 
vehicle  in  the  original 
vehicle  array.) 

CARD  10 

(required 

if  OPTI  =  1) 

1-12 

E12.5 

spe 

- 

Spherical  probable  error 

radius  in  feet 
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Appendix  D 


CSSANE  Prograir.  Source  Listing 


PROGRAM  GUIDE  ( INCUT , OUTPUT) 

£  **  **«.**«*  *****  *****  **  **  ***  *******  *****  *  *  ***  *  **  ***  **  ******  ************* 
C  FLOW  OF  ANALYSIS  IS  CONTROLLED  BY  FROGRAM  GUIDE. 

C  SUBROUTINES  CALLED  BY  GUI^R  ARE  TARGET , G A “ A , XR AY , THERM, PL ASTi , AND  NEUT 

COMMON/CAPRI/OISFR  (lJr.  )  ,PCSZ  (ILL)  ,KTEMP,  T  EMDS  ,  TM  (  3  )  ,  VEH  IC  ( 1  j  r  )  ,  1 1 
i  ,V9LAST1, VG,ICDAT 

COMMON/GRNO/HZVHG,H°,W 
COMMON/POS/0, HURSTX, BURST Y,9U°STZ 

COMMON/ INFO/ POSY  (11(1)  ,  pOS Y  ( 13 2  )  , VX ( 1 0 C ) , VY (lu ^)  , V7 ( 19 0) 

COMMON/BIT/  P0,VIP,RH0,7H 
DIMENSION  NUM  (10) 

INTEGER  COUNT, CYCLE, OPTi,CPT2 
C  ENTER  VISIT.  LEVEL,  H20  VAPOR  PRES.,  ALEECO,  HA7E  LEVEL 
C  ENTER  NUMBER  OF  VEHICLES  AND  NUM3-R  OF  pUcSTS 
C  ENTER  TEMPERATURE  VALUES  IF  NCN-ST ANDARD  ATMOSPHERE  DESIRED 
C  (ENTER  ZEPOES  FOR  STC.  ATMOS.) 

READ  22,VIS,P0,RH0,ZH 
READ  2l,II,PYrLE 
READ  4U,KTEMP, TEMPS, TM 
COUNT--  3. 

DO  lCu  1  =  1, IT 
VEHIC \ I) =1 

190  RFAO  22,POSX(I) ,FCSY(I),PCSZ(I) ,VX(I) ,VY (I) , V 7 ( I ) 

READ  21, (NUM(I)  ,1  =  1, lr  ) 

PRINT  25 
PRINT  26 
DO  1C1  1=1,11 

101  PRINT  27,I,P0SX(T)fP0SY(I)  ,POSZ(I)  ,  V  X  ( I )  ,VY(I),VZ(I) 

PRINT  62 
PpINT  84, CYCLE 
PRINT  86, VIS 
PRINT  3 8 , RHO 
PRINT  qp,PO 

If  (KTEMP.EO.O)  GO  TO  102 
PRINT  92 
PRINT  96 
PRINT  97 

PRINT  94, TEMPS, TM(l),TM(2),TM(3) 

102  PPINT  93 

ENTER  VULNERABILITY  LEVELS  FOR  TYPE  VEHICLE 
E_NTER  ZERO  FOR  EFFECTS  TO  BE  CELETED 

READ  22,VGAMA,VXRAY, VT HE RM ,VBL ASTI, VN PUT, VO 
STATEMENT  104  BEGINS  LOC°  THAT  CYCLES  ONCE  FOR  EACH  BURST 

104  READ  24,TIM5X,W , HG, OPT l,IOQA7,OP72 
R  =  G 

COUNT  =CO'JNT  +  i « 
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NEXT  in  LTN^S  UPDATE  VEHICLE  POSITIONS  TO  TIME  OF  BURST 
TIMEX  IS  INCREMENTAL  TIME 

i j 5  PRINT  10, COUNT 

PRINT  26 
AA  =  0. 

00  109  1=1,11 
IF(VEHICCI)) 1.7,109 
1'7  POSX(I)  =  POSXm+VX(I)*TIMFX 
POSY  ( I )  =P0SY  ( I)  +  V  Y  ( I )  TIMEX 
P0SZ(I)=F0SZ(I) fV7(I)^ TIMEX 
AA=A A+i. 

PRINT  27,1 , POSX ( I ) , PQSY(I) ,POSZ(I) ,VX(I) ,VY(I),VZ(I) 

IF  (POSZ(I).GE.HG)  GO  TO  139 
VFHIC(I)=. 

PPINT  72,1 
1C  9  CONTINUE 

IF(AA.EO.O)  GO  TO  311 
llu  RFB=113* (W**.A) 

IF  ( OPT 1 )  3u  5, 3  J  0 

BRANCH  3r'  G  FOR  US^R  PLACEMENT  OF  «URST 

3‘C  READ  22, BURSTX, BURSTY, 9URSTZ 
IF  (OPT2. EG. 1)  GO  TO  3"3 
00  3C2  1=1,11 

3"2  RPAO  22,VX(I> ,VY(I> ,V7(I> 

3o-  3  PRINT  73 

PRINT  78,BURSTX,3URSTY,3URSTZ 
PRINT  79, HG 
PRINT  83, TIMEX 
GO  TO  333 

BRANCH  305,  SUBROUTINE  TARGET  PLACES  BURST  ACCORDTNG  TO  A  SPHERICAL  LY 
NORMAL  DISTRIBUTION 

DIRECT  ATTACK  ON  FIRST  OCCUPIED  SLOT  IN  VEHICLE  ARRAY 

3  35  READ  22,  SDE 
DO  310  1=1,11 
IF (VFHIC (T) ) 30  6 , 31T 
786  TARGX  =  POSX  (T) 

TARGY-=POSY  (I) 

TARG7=POS7 < I > 

GO  TO  315 

310  CONTINUE 

311  PPINT  35 
GO  TO  9j  3 

C  IRC  USED  TO  PRECYCLE  THE  RANDOM  NUMBER  GENERATOR  IN  TARGET 
C 

715  IRC=NUM ( COUNT) 

If  (SPE.LT.1GC)  GO  TO  322 
PRINT  72 
PRINT  77,1 
PPINT  79, HG 
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PPINT  8", TIMEX 

CALL  TARGET (SPE , T ARGX , T APGY, T ARG7, IRC ) 

GO  TO  323 

C  FOR  SPE  LFSS  THAN  100,  PLACE  EUPST  ON  VEHICLE  FOS1I7ION 

322  BURSTX=TARGX 
BURST Y  =  T  ARGY 
BURSTZ=T  ARGZ 
PRINT  76 

323  PRINT  74, =>,PUPSTX, BURSTY, PURST7 
33c  PRINT  28 

PRINT  29, COUNT, W 
PPINT  28 
H8=BUPSTZ 

C  FIND  SLANT  RANGE  TO  EACH  VEHICLE  RFMAZNIN^  AND  CHECK  FOR  FIRFDfiLL 
C  INTERSECTION  (NEXT  10  LINES) 

00  ?5c  1=1, IT 
IF  €  VEHIC  C  T  > )  748,350 
340  OISPX= (POSX (I)-°URSTX) **2 
OISPY= (POSY(I)-^URSTY) **2 
OISPZ=  (POSZ(I)-®U°FT7>  **2 
DISPR (I) =SCRT (OISPX+DISPY+DISP7) 

IF(OISPR(I) .uT.RFP)  GO  TO  350 
VEHIC  ( I )  =  .' 

PRINT  3  0  ,  I ,  RFB 
350  CONTINUE 

COMPUTE  EFFECTS  LEVELS  AND  DETERMINE  LOSSES,  (NEXT  10  LINES) 

112  IF(VGAMA) 115,120 

115  CALL  G AMA(II,PURSTX, BURSTY, 3URSTZ,PCSX,P OS Y,PCCZ,VGAMA, VEHIC) 

120  IF(VXRAY)125,130 
125  CALL  XRAY 
130  IF  (VTHERM ) 135,140 
175  CALL  THEDM  ( VTHERM , DF ^ ) 

140  IF(VPLASTl.r_fJ,?  ,  ,  ANO .  VO .  EO  ,  0  .  )  GC  TC  170 
145  CALL  RLAST1 
17 C  IF(VNEUT) 175,19 0 
175  CALL  ME'JT 

LOSSES  EACH  EFFECT  FRINTED  IN  SUBROUTINE 
190  CONTINUE 

IF  (CYCLE -CCUNT)9QC, 20  0,104 

DETERMINE  AND  PRINT  ALL  VEHICLES  LOST  IK  ENTIRE  MISSION  PHASE 

20  0  PRINT  60 

DO  2tf  1=1,11 
IF  (VCHIC(I)  >9'10,2C5,2G6 
205  PPINT  62,1 
276  CONTINUE 
9 j  0  STOP 

10  FORMAT  (1H1,UX,* . VEHICLE  POSITIONS  AT  TIME  OF  3UPST  NUMBER*, 

113,*  — - *,//) 
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21  FORMAT  ( 1° 13) 

22  FORMAT (7E12. 5) 

24  FORMAT  <3E12.C,  312) 

25  F0RMATC1H1,25X,39H - VEHICLE  INITIAL  CONDITIONS - //) 

26  FORMAT  OH  VCH ICLE  ,  i 2X  ,  15HF CS ITI ON  (  X  ,  Y  ,  Z )  ,  *  ( FT  ) *  ,  liX  ,  2 7H  VcL0C I TY  ( X  , 
1Y,7)  (FT/SEC)  ) 

27  FORMAT(3X,I3,6X,TF10.2,3X,3Fir.2) 

28  FORMAT  (/,1X,  4 7  ( 1H*)  , / , IX , 47 ( 1H* ) ) 

29  FORMAT  (/7X,*  LCS^E^  CROM  NUMFEc*,  13, 2X  ,*  (*,  F7.  1,  *KT)  *  ) 

30  FORMAT  (//, 1 3X ,* VEHICLE  NUMBER* , 13 » *  IS  INSIDE  THE  FIFrQALL  RAOIUS 
1*,F5..,*  FEET,*,/1CX,*  ANO  HAS  BEEN  REMOVED  FROM  FURTHER  CCNSIOERA 
2TION. * , //) 

32  FORMAT  (//,i:X,* - VFHICLE  NUM3Ep*jI3,*  HAS  INTERCEPTED  *,/, 

114X,*GRGUND  LEVEL  ANO  HAS  FEEM  REMOVED.*) 

35  FORMAT  </,lX,*ALL  VEHICLE  SLOTS  ARE  EMPTY,  °RCGRAM  IS  TERMINATED.* 
1) 

40  FORMAT  (Il,4F6.n> 

60  FORMAT  (//,4'H  V EHICL-S  LCST  OVER  ENTIRE  MISSION  PHASE) 

62  FORMAT ( 12X , I  3) 

70  FORMAT  (/////, IX, 134(1 H- ),/, IX, 1 34 <1H*m  ,/, IX, 13 4(iH-) ,/,* - WE A 

I® ON  BURST  POSITIONED  BY  USER - *) 

72  FORMAT  (/////, IX , 134 ( i H- ),/, IX , 134  < 1H-) , / , IX , 13  4  ( 1H- )  ,/,  * - WEA 

1PON  BURST  POSITIONED  dY  SUQOOUTINE  TARGET - *) 

73  FORMAT  (1QX,*TUE  VEHICLE  IS  VEHICLE  ,'/*,I7) 

74  FORMAT'  3X,*TURST  DFT  ONATION  AT*  ,  c  1 0  .  2  ,  *  FT  FROM  TARGET  VEHICLE 

1  ,/,  8X , *NUCLE  AR  DETONATION  COORDINATES  A  RE*  ,  3F1 j • 1 , /) 

76  FORMAT (  3X,*--THE  SPE  IS  LESS  THAN  1  j  FEET,  RURST  HAS  PE^n  PLACED 
1  AT  THE  TARGET  COORDINATES.--*) 

78  FORM  \T  (  8X  , * NUCL E  AF  DETONATION  COORDINATES  AR E* , FI 3 . 1 ,* , * , Flu . 1 , * , 
1*,P10.1) 

79  FORMAT  (3X,*GR0UNB  HEIGHT  A°OVE  SEA  LE VEL  =  * , 4 X , F 9 . 1 , *  FT*) 

80  FORMAT  (3X  ,*  INCREMENTAL  TIME  OF  BUR ST  =  * , 8X , F3 . 1 , *  SEC*) 

52  FORMAT  (// /P X  ,* INPUT  °  AR  AM  FTE pS--* ) 

94  FORMAT  (10X  ,*N’UM3ER  OF  UUPSTS  ENTERED  I?*,9X,I2) 

96  FORMAT  ( 1 QX , * VIS  I B ILIT Y  AT  SEA  LFVEv  I . )  IS*,2X,F3.2) 

38  FORMAT  (1JX  ,*ALBEDC  (GROUNC  REFLECTANCE)  I3*,4X,F8.2) 

90  FORMAT ( 1  OX ,*H?0  VAPOR  PPESSURE (SEA  LcV£L,MMHG)*,F3.2) 

92  FORMAT  (13X , *  NON S T ANTA PC  A7MOSPHrRE  ENTERED*) 

96  FORMAT  (IjX  ,*TEMBEpATJRES  IN  C“GPEFS  RONKINE  — *) 

97  FORMAT  </,loX,*AT  SEA  LE VE L*  ,  3X  ,  * AT  R c CEI VER* ,3X , *  AT  GROUND* , 3X  ,  *A 
IT  BURST  * ) 

94  FORMAT  (10X,4(2X,  F13.1)) 

98  FORMAT  ( IT X ,  *ST  ANDAP'O  ATMOSPHERF  SPECIFIED*) 

END 
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SUBROUTINE  TARGET  < SPE ,X , Y ,7 , IPC  ) 

C  SUBROUTINE  TARGET  POSITIONS  THE  BURST  ACCORDING  TC  A  NORMAL  DIS7RI9U- 
C  TION  OF  ERROR 
C  calling  ROUTINE  IS  GUIOP 

COMMON/PQS/R, BURSTX, BURSTY, 8UPSTZ 

sumx=;. 

SUMY=  G  • 

SUMZ=n. 

C  XYZ  IS  VEHICLE  COORDINATE  POINT  ON  WHICH  EURST  IS  TARGETED 
SIGMA=SPE/1.5?3 
00  10  1=1, IPO 
10  Xi=RANF{OUH) 

DO  20  1=1, 2F 
X2=RANF  (DUN) 

IF (I«  GT» 12 )  GO  TO  14 
SUMX=SUMX+X2 
GO  TO  20 

14  IF ( T « GT, 24 )  GO  TO  16 
SUMY=SUMY+X2 
GO  TO  2? 

16  SUMZ=SUM7+X2 
2C  CONTINUE 

RX=( (SUMX-6.) ♦SIGMA) **2 
RY=( (SUMY-S. >*SIG*A)**2 
RZ={ (SUM7-6.) *SIGMA)**2 
R=RX+RY+R7 
RX1=SPRT ( °-R Y-R  Z ) 

RY1=S0RT (R-RX-R7 ) 

RZl=SnRT (R-RX-RY) 

R=SORT (R) 

8URSTX=X+RX1 

BURSTY=Y+RY1 

BURST7=Z+RZ1 

END 
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SUBROUTINE  BLAST i 

Q  **:*9******.e**********$.******.***¥***-**¥*******'*****-$-#******IMMHfS»S-J!!-J 

C  SUB°OUTINE  8LAST1  CCHFUTES  THE  BLAST  EFFECTS  PARAMETERS  AND  MAKES 
C  SURVIVAL  TESTS 


C  CALLING  POUTINE  IS  GUIDE 

C  BLAST!  CALLS  SUBROUTINES  TRIPNT,  MOTION,  SETUP, MAC'JEE ,  AND  ATMOS 

COMMON/C  ARRI/OISPS  (10.),  PCS?  ( 1-j  T  ) ,  KTEMF ,  TEMP$,TM(3)  ,VEFT0(1G  )  ,1 
1  ,V8LAST1, VOtICOAT 

COMMON/GRNO/HZ,HG,HP,w 
COMMON/BIN /POSX 1 , POSY  1 

COMMCN/TRI/PZ,PG, cp,  PZR,PGR,P8R,  SR ,FR ,SRE, ALP1ER , ALP1E, 

1  R90,PBPZ,R,PBPW, INCOME 

COMMON/SENSE/CFS 

COMMCN/COLC/CFii<2) ,Cr22(2) ,CF33(2)  ,CF44(2),CF55(2),CFR(6>,CFi(7) 
1  CF2(7) ,CC7  (7) 

COMMON/TAB/  TABU  (69)  ,  TABID  (69)  ,  TAB2I  (62 )  ,  T  A9  2C  ( 62 )  ,  T  AB3 1  ( 69 )  , 

1  TAB  30 (69) ,  TAS4K18) , TAB  40 (18) , TAB  FT (26) ,T ABED  (26)  , 

2  TAB6I  (69) ,TAS6H(S9) 

COMHON/CON/UL2(7) ,UL3(8),UL4(8),UL5(5), 

1  C4(7) ,CC(7) ,C6 (?) ,C?(S) , C 8 ( 3 ) ,C9<8) ,Cl3  (S),Cil  (8)  , 

2  Ci?.(8)  ,P2(8)  ,P3(3)  ,A1E  (41)  ,ACP(13)  ,CF(13) 
COMMON/PNT/RA,STC,PFI?,SHE,$T,XKK, ALTSRG,HT,XKKX 
CONMOM/GVP/DELD,QELPO,OELPR,CONT,NIT,XIT£R 
C0MM0N/T1LKIJP/L1,LF,NA  (6 )  , XL  ( 1C  * )  , NNEX 

DIMENSION  F(?),H<3),PR<?>  , TEMP ( 3) , RFO (3)  ,SS(7) 

EOUIVALENCE  (P,PZ) , (H,HZ> , (PR,P2R) , (RH0Z,RH0(1) ) , (RH0G,°K0(2>  )  , 

1  (RHOB,PhO(3) ) , (SBZ,SS(1) ) , (SS0,SS (?) ) , (SSB,SS(3) ) 

DATA  CFii/. 599829 71^,. 5999J 26553 3/ 

OATA  CF22/-. 81312 116 J,-. 3 3 4545670° 2/ 

DATA  CF 33/. 06 3796 96 84, .23 153383 : 84/ 

DATA  CF44/.  1735 12  593,-. ’■'42235;  7  7752/ 

OATA  CF55/o  137932961,.  3  8  26  3  3  228  6  7/ 

OATA  CFS/1.85/ 

DATA  CFR/. 5, .3, 1.2,1. 6,3.6.13./ 

DATA  CF1/0 .,-,6444,-.  1525,-. 05,-. 3  J7,  .  00  01111  ,C./ 

OATA  CF°/. 6499, 1.2943, .4725, . 2  . 5 3 , . C 4 24 , - . J 0 3  466 , 0 ,  / 

DATA  CF3/. 007C75, -. 15^024, .195,. 354, .5C408,. 577 04, .591/ 

R9Q  =  9r.  .  0/5  7.296 
DO  10  J=l, 18 

TAB4I <J)=T AB4T( J) /57,°96 
10  TflB4D(J) =TAB4D(J) /57. 296 

NN=?n 

DO  2C  J=1 , 69 
N=NN-J 

TA«3T(J)=TAB1D(N) 

20  TAB3D( J) =TA01I(N) 

HOR  =  W 

PRINT  1555 
N  =  i 

30  IFCVEHIC(N))  34,32 
32  N=N+1 

IF(N.GT.IT)  GC  TO  1130 
GO  TO  30 


78 


GNE/PH/72-3 


34  SR=OISPR(N) 

H7=P0SZ(N) 

PSL=2il6. 217/144. 

NIT  =  P 
XITPP=0. 

XKKX-r.. 

CONT=''  • 

DELPD=9. 

INCOMP  =  '< 

SPE=G . 0 

ALP1E=G.9 

ST0=9 . 0 

RA=0 • P 

ST-C.O 

HT=G,C 

TSA=O.G 

KCASE-0 

OELP=C 

SRSV=SP 

120  IF  (H3.LT.25J09.)  GO  TO  130 

CALL  SETUP  (ACF,l,2,l?,j,o,0,d,0> 

CALL  MACUPE  (CP,H3,j, :  ,j,0,0,LrR,CFF) 

W=CFF*WOR 
130  KCASP=KCASE+1 

21u  IF  (HB«LT,250o00«)  GO  TO  22" 

PRINT  1290 
PRINT  1290,  H0 

GO  TO  113' 

220  ALTM=A8S  (HZ-rP) 

IF  (SRSV-ALTM)  240,25?,2EO 
240  PRINT  14  j  ?  ,  HB, HZ, SR 

GO  TO  lu6j 
25C  HZi=l . 

HZ§AV=H7 

INTC=V 

FLAG=r. 

SFP=0. 

OT=.  1 
TT= .  0 

SR=SR/1C  00 . 

GO  TO  270 

260  IF  (HZ. EQ «H71 )  GO  TO  4^0 
HZ=HZ1 

270  OO  290  J  =  1 v 3 

CALL  ATMOS  (H ( J ) , TEWP ( Ji  , DEN , RHO ( J)  , TR , PR { J) , SS ( J) , VC ,KER) 
IF  (KER.NE.l)  GO  t e  3.0 
290  P(J)=PR(J) *P^L 
GO  TO  351 

300  PRINT  127.,  J , KER 

GO  10  Id 6C 

350  IF  (KTEMP.EQ.Q)  GO  TO  43" 

PSL=14.696*TEM°S/51P.67 
360  00  370  J-1,3 

P(J)=P(J)VTEMP(J) /YM<J> 

PR ( J>  =P ( J> /FSL 
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RHO ( J ) =RHO ( J) *TN ( J) /TFMP ( J) 

37T  SS(J)=49.32*TM( J> **.5 

470  pnpz=PZ/P° 

SSZR=SSZ/illo,4437 
PPRW= (PPR/W) **. 327733 
PZRH= (P7R/W) **. 33733T 
443  IF  (SR. MR. 3..)  SO  TO  49j 
PRINT  116° 

GO  TO  1C  60 

490  CALL  TRIPNT  (KCASEV 

IF  (INCOMo.EQ.l)  GO  TO  1060 
550  O0RWFR= (PRR/ (W*FR) ) ** . 333333 
P7RWFR=(PZR/(l-i*FR)>  **.  3333  37 
RBAR=SR*PPRWFR 

CALL  SETUP  (TAT11, 1,2, 69,0, o5C  ,uT) 

call  MACURE  (TA°10,P9AP,0,0,u,:,0,LER,S')ELP) 

IF  ( H7 . NE . H° )  GC  TO  5o0 
DELF=SOELP*PRR 
ALFA  =  i.  j 
GO  TO  570 

56C  CALL  SETUP  < TA32I , i , 2, F? , 0 .  1 , - , o  ,  0) 

CALL  MACUPE  (TA320,P9AR,C ,0,3, C,0,LER, ALFA) 
0ELP=S0ELP*Ppp*F3F7**ALcA 
57°  CPSTL0=0ELP/P7 

SPV  =  SS2*(1.0<-6.0*rPSIL  0/7. '■')**<.? 

Q=2.5*0ELP**2. 0/ (7«:*nZ+nELP) 

CALL  SETUP  (TA96I,l,2,69,o,C,C,0,r) 

CALL  M AC Upc  (TA360,R3AP,C  ,C,0,C,'JjLFR,TSACAP) 
TSA=TPACAP/(SSZP*F7PWCD) 

CALL  MOTION(N,SFV,SP,HZ1,OT,SFR,FLAG,SEPO,TT,P7SAV,SEPOT,INTC> 
IF  (INTC.LC. 1 )  GC  TO  c7i 
VFHIC (N) =0 
PRINT  156^, N 
GO  TO  lu’C 
571  SP=FR/10r. 

IF(FLAG.LT,3.)  GO  TO  r'72 
IF (SEPOT .L  T  .  .  • )  GO  TO  574 
PDINT  149  ;,N 
GO  TO  1 j  70 

57?  IF (APS (SEP C) .GT.2CP .)  GO  TO  2Fn 
574  OO  58  ?  J=1 ,7 

IF  (PBAR.LT.UL2 ( J) )  GO  TO  590 
580  CONTINUE 

TOPZ=.252609+(1. 3/11.21) *ALOG(PPAP) 

GO  to  6^0 

590  TOPZ=C4 ( J) *P9AR**2+C5 (J> *RBAR+C6(J) 

600  00  61?  J  =  1 , 8 

IF  (RPAR.LT.UL3( J) )  GO  TO  620 
610  CONTINUE 
J  =  8 

620  TMV=C7 ( J)*RPAP**P2( J) +C8 ( J ) *Rp AP+C9 ( J ) 

OO  63c  J  =  1 , 8 

TF  (RRAP.LT.UL4  ( J) )  GO  TO  64  1 
630  CONTINUE 
J  =  8 
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640  TPRZ=Ci: ( J) »pqA9¥*p3( J) +C11( J) *PBAR+C12( J) 

P000P.=  T0P7/(SSZR*CZFW"D 

PMV=5. 0*EDSILO*SS7/ (7. C* (1,0*6. 1 *EPSILC/7 . -) * *. 5) 

POMV=TMV»POOOP 

P0D-RH07.*  (  7 .  C  +6 «  ^  *EPSILO)  /  (7  .  +  FPSILC) 

POOO=TO°7*COOCP 
SP=SP*1‘ Jr . 

IF(V8LAST1.E0.3.)  GO  TO  660 
IF(VRLASTl.GT.OELP)  GO  TO  660 
V  EH  I C  ( N )  =  *  . 

PRINT  150- ,N,OELP,VqLASTl 
GO  TO  6*0 

660  IF(VQ.EQ.  .)  GO  TO  60  ■ 

IF  (VQ.GT.Q)  GO  TO  63" 

VFHIC ( N> =" , 

PRINT  155-,N,n,vQ 
68G  IF  { ICO AT . EO . 3 )  GO  TO  Iu60 
730  PRINT  143", ‘5 

PRINT  144' ,  P  O  S  X 1 , D  O  S  Y 1 , HZ  1 

PPINT  1450,H7,SR,TSA,SFV,G,DELP,PMV ,FCO,PCOO 
PPINT  146T , FOOOP , PDMV 
IF  (Hn  ,LT , 25u C  3  « )  GO  TO  740 
GC  TO  1060 
74  0  PPINT  132** 

1  j  6°  CONTINUE 
1070  N=N+1 

IF(N.LE.II)  GO  TO  30 
1130  CONTINUE 

1160  FORMAT  (1H  ,////, 1CX, *  OHSL ANT  RANGE  FCUNO  EQUAL  TO  ZERO  IN  OVERFRE 
1SSURE  SOLUTION  PROCEEDING  TO  NEXT  CASE) 

1270  FORMAT  (1H3,3HATMCS  ?.R,I2,2H  =  ,I3> 

123o  FOPMAT  ( 1H J , /// , 4 8X , 37H9LA?T  EFFECT^  ARE  ESSENTIALLY  RECUCpD , // ,4E 
1X,24HT0  ZERO  AT  THIS  ALTITUDE) 

1290  FORMAT  ( 1H 3 , 47X , 5HHB  =  ,1PE12.5) 

132"  FORMAT  (1H0,22X,39HYIELQ  CORRECTION  FACTOR  IS  EQUAL  TC  ONE) 

1400  FORMAT (/////45X, 4 ; <1H*) ,///45X,U3H*  AM  OVERPRESSURE  SOLUTION  C l 
1MNOT  *,/,45X,40F*  RE  OBTAINED  WITH  THE  GIVEN  *,/,<*5X,< 

20H*  INPUT  GEOMETRY  * , // ,6 OX , 4HH E  =,E12.6,/,< 

3GX,48HZ  =,E12.5,/,6tX,4HSR  = , El2 . 5 , // ,46 X, 40H *  THE  PROGRAM  WILI 

4  PROCEED  WITH  *,/,45X,4CH*  THE  NEXT  CASE 

5  *»///45X,4" ( l  H  * ) ) 

1430  FORMAT (//,5X, ♦SHOCK  FRONT  PARAMETERS  AT  INTERCEPT  OF  VEHICLE  NUMB! 
1R*» 13) 

1440  FORMAT  (5X ,* VEHICLE  POSITION  AT  SHOCK  INTERCEPT (X,Y»7)*,?X,3F!0.j 

1//) 

1450  FORMAT (1H j ,DX, 50HRECEIVE=  HETGHT  AT  SHOCK  INTERCEPT  (FT) 

1  ,F1C ,4,//,lCX ,5  .HSLANT  RANGE  AT  SHOCK  INTERCEPT  (FT) 

2  ,FlO.4,//,10X,p0HTlME  CP  SWCCK  ARRIVAL  (SEC) 

3  ,F10.4,//,i  )X,FpH'?HOrK  FRCNT  VELOCITY  (ET/SEC) 

4  ,Fl'**4,//,lJX,EvjPPEAK  DYNAMIC  FRESSUPE  (FSI) 

5  ,Flw  .4,//,l"X,5-’HPEAK  OVE ppRESSURE  (PSI) 

6  ,Fi0«4,//,lCX,5CHPEAK  MATEpIAL  VELOCITY  (FT/SED 

7  »F10*4,//,1vX,5jH(-EAK  OVERDENSITY  (SL UGS/ET**  3) 

8  ,F10«4,//,10X,5CH POSITIVE  CtRATION  OF  OVERDENSITY  (SEC) 

9  ,Fir,4) 
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1460  FORMAT  (1HO,3X,5QFPOSITIVE  DURATION  OVERPRESSURE  (SEC) 

1  ,Flv.4,//,l;X,5GhP0SITIVE  DURATION  MATERIAL  VELOCITY  (SEC) 

2  ,F1G.4,/) 

1490  FORMAT ( l^X  ,* - VEHICLE  NU«°ER* , IT , *  IS  OUTPUNHING  SHCCK  FRONT,*) 

15'.C  FORMAT  (10X,* - VEHICLE  NUMBco* ,  13  ,*  HAS  BEEN  SUBJECTED  TO  AN  OVER 

1PRESSURE  OF* ,/14X ,F1I .4,*  AND  WAS  BEEN  REMOVED  FROM  T  w£  PROGRAM,*, 

2/,10X,* - THE  OVERPRESSURE  VULNERABILITY  LEVEL  IS*,F10.4) 

1550  FORMAT  (IjX,* - VEHICLE  NUMSE0*  ,  IS ,  *  HAS  BEEN  SUBJECTED  TO  A  OYNAv 

lie  PRESSUPE  OF*,/, FI’.. 4,*  AND  HAS  BEEN  RE^OVEC  F^OM  THE  PROGRAM  ,*, 

2/  ,10X ,  * - thp  DYNAMIC  VULNERABILITY  LEVEL  IS*,C10.4) 

1555  FORMAT  (///, IX ,* RESULT S  FROM  BLAST  EFFFCTS  COMPUTATIONS*) 

156o  FORMAT  (// ,i:X,  — -VEHICLE  NUM3FR*  ,  P  ,*  HAS  INTERCEPTED  *,/, 

114X,  *GROUMC  LEVEL  AND  HAS  eEEN  REwOVEr.»  *) 

END 


o  o  o  o  o  o 
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SUBROUTINE  THEpM  (VTHERM,RpB) 

Q  V-*  W  WWWV-  y-ww  W  ■&&**.**■  +  **  9  +  if.  9.  WW  w  WW  ¥+*****■**■***■  V*  * 

C  SUBROUTINE  THPRM  COMPUTES  THERMAL  ENERGY  LEVELS  FCP  EACH  VEHICLE  AND 
C  MAKES  A  SURVIVIVALBEITY  D- TERMINATION. 

C 

0  THERM  IS  CALLER  °Y  ROUTINE  PLAST1 

C  SUBROUTINES  CALLEO  BY  THERM  ARE  TRANS, ATMCS, SETUP ,MACURE 
C  *****************  ***♦»»■*****************#¥************■*****  ***#♦***■**< 
COMMON/BIT /  P0,VIS,PH0,2H 
COMMON/T AR  TH/  TO M 1 ( 12 > , T pT 1 ( 12 ) , TBT 2  ( 1 2) , 

1AAA1  (22)  ,  AAA?  (12  0)  (39)  ,B32  (3S"')  ,TEV1(8)  ,TPV2(8) 

COMMON/CAPRI /niSP3  (U  .  ) , PCSZ  (1  uS  ) , KTEMP,  TEMP S, TM ( 3)  ,  VEH  IC  ( 10 C  )  ,  I! 
1  , VBL AST 1»  VG, ICHAT 

COMMON/GRNO/HZ, HG ,H°, N 
COMMOM/POS/R,PURSfX  ,  BURST  YjBUR'JTZ 

COMM ON /I NFC /POSY ( 1L D , FOSY (10 : ) , VX { 10 u ) , VY (10  0) ,VZ CIO") 

CHK=G 

WMT=VI/ 10  j  C 
PRINT  813 
HOB=Hn 

3  CALL  ATMOS(HOB,tEmF,OEN,RHO,TR,PR,SS,VC,KER) 

IFiKER.NE.l)  CO  TO  5 

IF(KTEMP.EQ.n)  GO  TO  7 

OEN=OFN* (TN(3)/TENP)* ( 518. 57 /TEMPS) 

GO  TO  7 

5  PRINT  fl-j  j  ,  KFR 
GO  TO  3 j  0 
7  CONTINUE 

TMIN={6J  .*  <V»MT**.4>  )  /100  0. 

TMAX2  =  .S j* (WUT**. 42)* ( DEN*  *, 42 ) 

C  DETERMINE  THE  THERMAL  YIELD  (NEXT  13  LINES) 

H=  H5*3.C48E-4 

TEFF=E.'.P  ( (-3.57R7123E-1-8. 8Q 4057E-3*H +7 , 13 ' 80  1E-4*H*H-1, 254  80  "E-5 
1H*H*H+6.423205E-8*H*H*H*H) *2,3^.2585893) 

I F ( (i,5*Rr B) .LT.HB)  GO  TO  15 
C  GROUND  RU°ST  CONDITION 
WTH=«32*W*1.E12 
PI=2,*3. 1416 
CHK=i 
GO  TO  20 

AIR  PUR^T  CONDITION 

15  WTH=<TEFF*W**.94)*i.E12 
PI=4«*3» 1416 
CHK=  2 
20  N  =  1 


78  BEGINS  LOOP  TO  TEST  EACH  VEHICLE 

78  IF(VEHIC(N)) 85,80 
80  N=N+1 

I F ( N , G T ♦ 1 1 )  GO  TO  300 
GO  TO  78 
85  QMHp=G 


ooo  o  cr  ooo 
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o°=n 

C  VEHICLE  FOSITTCN  ANO  SLANT  RANGE  TO  CURST  AT  TIME  OF  2N0  *AX. 
PX=P0SX(N) +VX  (N) *TMAX2 
PY=P0SY ( N) +  VY (N) *TMAX2 
P2=P0SZ(N) +VZ(N) *TMAX? 
l)ELX=A33  (nURSTX-PX) 

OELZ=BUPSTZ-PZ 
DX=(PURSTX-PX)**2 
0Y=(PURSTY-PY)**2 
0Z=(6URSTZ-°Z)**2 
DZZ= (PZ-HG )  **2 
SRR=S0RT (0X+0Y+0ZZ) 

SP=SOPf (OX+DY+OZ) 

SPMI=S‘R/5?8t 
SRCM=  v'SR*l2*2. 54)  **2 
IF<HB.GT.5 jOJo)  GO  TO  1)1 

DETERMINE  THE  REFLECTED  THERMAL  ENERGY  (NEXT  15  LINES) 

IF(HR.GT.3G:  )C)  GO  TO  41 
IF (CHK.EQ. 2)  GO  TO  25 
SPN=SD/RFf? 

PFE=AC3S( <PZ-HG)/SP) 

CALL  SETUD  (AAA1,1,3,12,1C,0,C,C,0  ) 

CALL  MACURE  (AAA2,SRN,PFE,0, 0 , ? , 3 ,L ER  ,GAM) 

GO  TO  30 

25  XRO=CELX/(HP-HG) 

HRO= (PZ-HG ) / (HP-HG) 

CALL  SETUP  (3Bl,l,3,lP,2C,J,0tC,0) 

CALL  MACUPE  ( 092 , XPO . «pO , D , w , G , C ,L ER , GA M ) 

30  TYPE=1 

CALL  TRANS  (TYFE , FZ , S R , SRR ,CHK , TR) 

REFLECTFP  THERMAL  FLUENCE  IN  CAL/CM2  FOR  AIR  OR  GRND.  PUpST 
SRRCM=  (SRR*12,*2.54)  **?. 

QR= (WTH*TR*RHO*GAN)/(PI*SRRCM) 

40  TYPE  =2 

CALL  TRANS  ( TYP E , FZ ,S R , PRR  ,CHK, TO) 

QOFF=  (WTH*TP) / (PI’SRCM) 

QDFFUN=WTH/(PI*SRCM) 

DIRECT  FLUENCE  ( CAL/CM2)  NORMALLY  INCITENT  ON  HORIZONTAL  RECEIVER 
QPN=OOFF* ( A35  (PZ-HP) /RR) 

IF(PZ»GT.13G'J3j)  PRINT  865 
IF(OELZ.GT.:>  GO  t0 
QMHR  =  QR4-QON 

IF  (VTHERM. GT.OMHR)  GO  TO  50 
VFHIC(N) =: 

PRINT  860, N,0MHR, VTHERM 
GO  TO  150 

50  IF(VTHERM.GT.OOFF)  GO  To  150 
VF'HTC  ( N)  =  0 

PRINT  86J,N,00FF,VTHEpM 
GO  TO  150 

COMPUTATIONS  FOR  HIGH  ALTITUDE  PURST  (HCB>50303  FT) 
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10 0  TF(PZ.LT ♦e>3C  30  GO  TO  120 
QnFF=WTH/(PI*SRCM) 

QDFFUM=QDFF 

QON=QDFR*( A0S<°Z-HP)/SR) 

IF(VTHERM.GT.COFF)  GO  TO  150 
VEHIC(N) -0 

PRINT  86Q,N,00FF, VTHERH 
GO  TO  150 

l?.u  CALL  SrrU°  (T«V1,1,2,3  »3» 3*0 , : tQ  ) 

CALL  MATURE  (T9 V2  ,H03 , c , «  ,  G , G , C , LER , FUV) 

STHET=  iB-P7)/SR 
SSTH=S7HET**2 

R=9987./STHET*(i.-EX?(-4.55E-5* (HB- 2 8 £  0 0  , ) ) ) 

IF  (PZ«GT«28)i3)  GO  TO  130 

R=R+31746. /STHET*  (EXP (-3.1 5E-5*PZ> 414) 

130  QDFF^WTH/ (SRCM* (H9-P7) **2> * (FUV* EXP (-1 . 3lE-4*R>  +  (1-FU V) *EXP(-1.3lf 
1-5*P) ) *SSTH 
QOFFUN=WTri/(PI*SRCM) 

QON=QDFF* (H8-P7)/SR 
IF(VTHERM. GT.noFF)  GO  TO  150 
VEHIC(N) =0 

POINT  863,N,0DFF, VTHERM 
150  IF ( IOOAT . £0 . i )  GO  TO  25^ 

Pr.'lNT  570  , N 
PRINT  875 ,  OOFFUN 

PRINT  88G,HB,^RMI,THAV?,00FF,nDN 
IF  ',HR»  GE  #3C0‘wU)  GO  TO  250 
PRINT  890 , OR , OMHR 
250  N=N+ 1 

''  IF(N.LE.IT)  GO  TO  78 
300  CONTINUE 

800  FORMAT  (1,X,* - AT  HOB  SUBROUTINE  ERROR  NUMBER*, IT) 

010  FORMAT  (IHJ,1X,*RESULTS  FROM  THr  RMA  L  EFFECTS  COMPUTATIONS*) 

860  FORMAT  ( 1H1 ,9X, * - VEHICLE  NUMBER*, I:,*  HAS  BEEN  SUBJECTED  TT  THi 

1RMAL  FLUENCE  OF* , /, 14* ,F1Q . 4, * (C AL/CM** 2 )  AND  HAS  BEEN  REMOVE C  FR’ 
2M  THE  PROGRAM*, /,1CX,*  — -THE  THERMAL  VULNERABILITY  LEVEL  IS*,F10 
34,*  (CAL/CM**?) *) 

865  FORMAT  (14X,MU3ER  WARNING-RECEIVER  ABOVE  I'D  KFT--TH ERM 6L  VALUES 
1COMPUTEO  THIS  CASt  HAVE  BEEN  SUB I£CTED  TO  EXCESSIVE  ATTENUAT I CN . ) 

2) 

870  FORMAT  ( //  ,  RX , *  THERMAL  EFFECTS  PARAMETERS  FOR  VEHICLE  NUMBER*, 13 
1//) 

830  FORMAT  ( 1H j , 9X, 50HHFIGHT  CF  BURST  (FT) 

1  ,F10 .4,//, i'X, 5 'HSLANT  RANGE  AT  TIME  OF  PEAK  RADIANT  POWER  (MI, 

2  ,F13.4,//,  nx,5CHTIME  TO  PEAK  R  A  C  I  ANT  POWER  (SEC) 

,  3  ,FlC.4v//,l'!X,ErM9Ir<FCT  RPEE  FJELO  FLUENCE  (CAL/CM2) 

4  ,F13.4v//,13X,5jH3IPECT  NORMALLY  INCIDENT  FLUENCE (CAL/CM?) 

1  5  »F1..  .4,/,  10X,  *  (NORMAL  TC  A  PLANE  HORIZONTAL  RECEIVER)*) 

890  FORMAT  ( 1W0, 9X, 50HREFLECTEO  FLUENCE  (CAL/CM**2) 

1  ,Flr  .4,//,10X,5'1HPE"LECTED  +  DIRECT  NORMAL  BLUENCE  (CAL/CM2> 

2  ,F1Q.4,/,10X,*  (CCMPUTEC  ONLY  IF  BURST  '-ELOW  RECEIVER'*) 

875  FORMAT  ( 1H 0 , 9X , 5 1 HUNAT TENL AT EO  FLUENCc  (CAL/CM2) 

1  ,  F 1  v  «  4 ) 

END 
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SUBROUTINE  tripnT  (KCASP) 

SUBROUTINE  TRIPNT  CALCULATES  LIMITING  ANGLE  FCR  REGULAR  REFLECTION 
AND  PREDICTS  WHETHER  RECEIVER  IS  IN  OR  OUT  OF  THE  FUScO  SHOCK 
REGION 

ROUTINE  REOUIP£MENTS-> 

NUMEROUS  PARAU£TE3S  FORM  MAIN  ROUTINE  THROUGH  COMMON 
CALLS  SUBROUTINES  SFTUP  ANO  MACURE 


5 

6 


CALLING  SFQUENCE 
WHERE- 

KCASE=1  FOR  OVERPRE^SLRE  SOLUTION 

2  FOR  TRIPLE  POINT  PATH  SOLUTION 

3  FOR  RANGE  SOLUTION 


'  *  *  *  4, 


CALL  TRIPNT (KCASE) 

COMMON/TRI/P7,PG,OP,  PZf-> ,  PGR ,  PPR  ,  SR,FR  ,SRE,  ALP1ER  ,  ALP1E  , 

1  P9C,F5°Z,R,PRRh, INCOME 

COMMON/GRND/HZ,HG,ho,  j 

COMNON/SENSF/CFS 


COMMON/COL  O/CFll  (2)  ,CF22  C  2 )  ,  CF33  ( 2 )  , C F44  ( 2)  , C  F55  <  2 )  ,  CFP  ( 6 )  ,  CF 1  { 7 )  ,j 

1  CF2 ( 7 ) , CF 7  (7)  j 

COMMON/TAR/  T A*U I (c9) , TAPIO ( 69 ) , T A9 2 1  (62 ) , T A8 2D ( 62) , T AD3I  ( 69 )  ,  j 

L  TA33C(69) ,TAP4I(1S) ,TA940(18) ,TAB5I (25) ,TA9ED(26)  ,  j 

2  TA36I (69) ,TABF0(69)  j 

C0MM0N/C0N/UL2(7)  ,UL3  (  8)  ,  UL4  (  8 )  , UL5  (5 ) , 

1  C4<7) ,C5(7),Ce<7) ,C7(8)  ,C?(8),C9(8) ,C1T (8),C11(P)  ,  1 

2  C12 (8) f P2 ( 3)  ,F3(  8)  , A1E( 41'  ,AOF(13) ,CF(13)  ] 

COMMON/PNT/RA,STO,PHIR,SH«,ST,XKK, ALTSRG,HT,XKKX  j 

COMMON/O VP/OELP , GEL  PC, PELPR, CONT, MIT, XI TEP 


OIMENSION  ALPHA (41) , I1 EOl (5,3) , ID (28) 

DATA (1HF01 (J) ,J=6,10) /5CH  TRIFLE  FOINT  PATH  SOLUTION 

1  / 

DATA  (IO(J), J=1,2«)/2uSR,2HHZ,2HHG,2HHB,1HW,3HTSA,2HFR,thScv, 

1  4HDCLF,3HPMV, 3HDC9,EHSDELP,4HRcAR,1HR,4H ALFA,3HSP£,5HALP1E,  ji 

2  3HST0,2HRA,2HST,CHHT,5HPC0OP,4HPCMV,4HP0Cn,lHQ,4HRH0Z,3HSS7,  \ 

3  3HCFF  /  J] 

PCPB=PG/P°  H 

ALT=(H3-HG)/1pJ3#T  1] 

R=i45.3*W**>.  ,4/n:o.')  I? 


SHB=ALT*PnRW 
SHBB=ALT/W**.3333 
IF  (SH3B.GT.2.5)  GO  TO  9" 
IF  (R.LT.ALT)  GO  TO  i 
FR=1  •  6 


t 


IF  (SP.EQ.T.)  GO  to  9 
XK=APS(HZ-HB)/(SR*1CT  .0) 

GO  TO  6 
XK=0 

IF((APS(XK-1.)) »LEo.Tj2)  XK=1. 
IF  (XK.LS.1.0)  GO  TO  180 
GO  TO  15 j 
XI  =  G  #  C 


RBAP=ALT’P8pW 
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CALL  SFTUP  (TAB2I,i,2,62,L,C,0,L,3> 

CALL  MACULE  (TA920,R3AR,C,G,d,,.,a,LER,ALFA> 

00  8 G  J"-l,41 
ALFA=r • 

XI=XI+0. 025 

IF  (XI.GT.l.Q)  XI  =  1 * 3 

n£LPG=PG/XT-PG 

2  P  S0ELP=DELPG/ (oqR*PGP9**ALFA) 

ODFLP=-DELPG*ALQG(pq°B)/ (FBR*PGPB**ALFA) 

K=1 

IF  (SDELp-CFS)  40,3r,?r 
3^  K=2 

4C  PPP=CF11(K)+CF22 (K) *AL0G1': (SQELP) +CF37 (K) * ( ALCG13 (SOFLP) **2) +  CF44 ( 

IK)* ( A  L0G1 .  (SOELP) **3)+CF55 ( K ) * ( A L0G1 0  (SOELP) * *4 ) 

PPAR-1 3 • **  FPp 
AA=ALOGlQ( 2, 71828) /SOELP 
0P=2.* AL0G1P (SOELP) *AA 
CC=3.*(ALCG1 j(SCSLP) **2) * AA 
QP=4.*  (AL0G1° (S0ELP)**7)*AA 

DPDX=CF22 (K)*AA+CF33(K )*9B+CF44 (K) *CG+CF55 (K) *DQ 
DRBAR=R3AR*AL0G(13. )*CpOX 
CF4=0 « 

00  50  IT=1,& 

IF  (RRAR-CFR(II) )  60,50,50 
50  CONTINUE 
21=7 

CF4  =  - . j  38 

6"  SHF=  (CF1  (II)  *R9AR+CF2  (II)  )*R3AP  +  CF3  (II) +CF-**ALOG13  <R**AR> 
CAPD=2,*CC1(II)  *R3Ap+CF2(  II)  +  CF4*ALCGK  (2. 71 8 28) /ROAR 
CAP0=CAP0*0oBAR*00FLP 
ALFO=ALFA 

ALFA= (3MF-ALFO*CAPO) / (l.-CAPH) 

IF  (ABS(  ALFA-ALFC) -.0  1)  73,7%23 
7  0  ALTSRG  =  4LT/ (ROAR*  (W/P^R) **. 337533) 

IF  (ALT3RG.GT.1.  5*)  GO  TO  85 
ALPHA (J) =ACOS (ALTSRG) *57.295 
IF  (ALPHA(J) .GT.AIE(J) )  GC  TO  120 
9T  CONTINUE 

9*  IF  (KCASE.NE.2)  GC  TO  li.u 

INCOHP=l 

PRINT  2J0, (IHF91 (J,KCASE)  ,J=1,5) 

PPINT  213 
PRINT  220, H8 
RFTURN 
1J0  FRrl.O 

IF  (XITER.GT.1..0R.NIT  »G7 « 1)  GO  TO  115 
110  XK=AB*{H7-HB>/(SR*lr>J0,0) 

IF  (  (ABS(XK-Ib))  .  LE » .  .Zc)  X«  =  i. 

IF  (XK.Lt.l.-)  GC  TO  130 
GO  TO  150 

1 2 j  IF  ( J . NL » 1 )  GC  TO  130 

Ai  P1F=A1E(  1) 

GO  TO  140 

133  ALP1E=  (  (AiE  (  J)  -  A1E  (J""i) )  *  (AL°HA  ( J»i )  -  Alt  ( J  - 1 )  )/(AlE'J)-Alr  ( J  - 1) - 
iPHA  ( J)  +ALnHA  (J-l)  )  f-Al"  (J)  > 


wooo 
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It*:  )  1.  PIER  =  M_PlE/57.296 

Sf<!E- ALT/COS  (  ALPIER) 

:r  (K'CASE.  EQ .  2)  GG  TO  180 
ALPHI=R90 

:f  ch7.eo.hp)  go  to  i?c 

<K=APS  (H7-HD)  /(SRM  "»')')•'») 

:r  c  (Ass<x<-i.) )  .LE..-  J2)  xxr=i. 

:r  < x k . l  .  l • ; >  r<Q  to  i6q 

.L *> I.  PRINT  240, YK 

PRINT  251, I 0 ( " ) ,W,ID(1) ,0ELP 
PPINT  253, IP (2) ,H2,in (3) ,rfi 
PRINT  250,10(4) , HP, 10(1) , SR 
PRINT  260 
OEIPO-OIL0 

i*  out  =  roNT>  i  * 

' RIM T  270 

<kk=xk 

<  K  K  x  =  y  K 
RETURN 
ic,'  «KK=XK 

*KsASlN(XK> 

IF  (M7.LT. HP)  4LPHI=ALPHI-XK 
IF  ( H7 . GT . Ha )  AL°HI=AL°HI+XX 
17 '.i  t  LPHTP=ALPHI-*LP1PP 

'  ALL  SETUP  (T*n4l,l,2»id,C»C»i.  ,4,0) 

!  ALL  MACU°E  (  T  A  24 0 ,  AL  3HI R  ,  y  „ ,  r  ,  0 , 0  , L CP ,  PHIR) 

='A=sPF*cns  ( At  pier-phi  3 )/ccs< At fwi-phir> 

A  =  A  P  S  ( R  A ) 

F'Pal.  C 

]'F  (XKKX.GT.  >  .)  GO  TO  l<*r 

::f  (pa.gt.sp)  go  to  ior 

CALL  SETUP  (T 4051,1, 2, 26, C,0, 0,0,0) 

CALL  MACURE  ( TA 1 5  0,  SHD  ,  -i  ,  4  ,  -  ,  0  ,  0  »LeR  ,  FP  > 

::f  ( shb.lt .i .54)  go  to  ion 

::F  (  (S9-9A)  ,LT.  J.l)  GO  TO  19  J 
FR  =  2«  33-  3  • 125*P3AR 
IV!  :<KK=XK 
190  RETURN 


220 
2  4u 
25" 
260 
27  J 


FORMAT  (lHi, 45X,5A1C ,///) 

FORMAT  (1H0, /// , 4 0X , 35HINPUT  PARAMETERS  ARE  NCT  COMPA TI BL E , // , 4 3Y 
134HFOR  THE  TRIQLE  POINT  PATH  SOLUTION) 

FORMAT  ( 1H  0 , 47X , 5HHB  =  ,lPr12.5) 

c0RMAT  (38H  ***  ARG  OF  ASIN  (X)  OUT  CF  RANGE.  X=, E16 . 8 , //// ) 
r0RM AT  (1H0, 22X,2(4X,A6,1PE12.5) ) 
rGRM A  T  (IX,/////) 

FORMAT  (IPX, 66HTHE  TNPUT  GIVEN  IS  INCCMP4TA9LE  WITH  A  POSSIBLE  PH' 
1SICAL  nONniTlON,//,lJX,39HTWO  ALTERNATE  SETS  CF  OUTPUT  ARE  GIVEN* 
?/,14X,61Hl-RECETVER  OIRECTLY  AEOVE  CP  PELOW  THE  BURST  DEPENDING  0 
3  THF,/,16X,57HINTTIAL  ORIENTATION  OF  RECEIVER  WITH  RESPECT  TO  THE 
4BURST  ,//,l4Y ,64H2-TH£  ALTITUDE  AT  WHICH  THE  DESIRED  GUST  OR  OVERF 
5ESSURE  OCCUR RS,//,34X,1H*,/,34X,1H*,/,34X,1H* ,/ , 34X , 1H* , / , 31X , 7H* 
(,  +  ***», /,32X,  5H*  ****./,  33X,3H***./.  3  4X.1H*) 

END 


88 
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SUBROUTINE  TRANS  (TYP- ,DZ s SR,SRR  ,CHK , TT) 

C0MM0N/TA3TH/  T3W 1 (12 ) , TBT1 ( 12 ) , T«T 2  ( 12)  , 
iAAAi  (22)  ,AAA2(12J)  ,P91  (39)  ,33?  (38')  jTEVl(3)  >TPV2\3) 
COKHON/BIT/  P0,VIS,RHP,7H 
COMMON  /GRND/  HZ,HG,HP,W 
ZS=HP 
HL=HG 

El=EXP(-4. 57E-P4?H) 

E2=EXP(-4.S7E-54FZ) 

E3=FXP(-4. 57E-54ZS) 

E4=EXP (-4.57E-54HL) 

IF  (TYPE. EG. 1)  GO  TO  330 

C  COMPUTATIONS  FCR  DIRECT  TRANSMISSION  (NEXT  23  LINES) 

C  NON-COALTITUOc  CASE 

CTHETA=(PZ-7S’,/SR 

IF  (APS(CTHETA) »LE..3"1)  GO  TO  30C 
IF(HB.LE.ZH)  GO  TO  53 
IF  (PZ.LE.ZH)  GO  TO  25 
TH=i. 

Tp=EXP(~. j85/CTHETA*(ET-E2)  ) 

GO  TO  100 

25  TH=EXD(-16.44SR/ < VIS* ( ZS-PZ) ) * ( E2-E1 ) ) 

TR=EXP  ( - « ..  8754SR/(ZS-DZ)4(E1-E3) ) 

GO  TO  13  0 

50  TF(PZ.LE.7H)  GO  TO  75 

TH=FXP 1-16.4/ (VIS4CTHCTA) 4 (E3-E1)) 

TR=EXP(.  0  8  75/CTHETAMEi-E2)) 

GO  TO  luO 

75  TH=EXP(-16.4/(VIS4CTHETA)ME3-E2) ) 

TP=1  • 

1)0  XXW=2.3*P0/CTH£TA4 (1344 <-6.1E-54ZS> -10**  (-6.1E-54P7) ) 

GO  TO  103C 

c  coaltitijoe  case 

330  IF(ZS.LE.ZH)  GO  to  325 
TH=1.C 

TR=EXP  ( (-4.E-6)  4^R4FXP (  (-4.57E-5)4ZS)  ) 

GO  TO  400. 

325  TH=EXP( (-7.5E-4) 4SR/VIS*EXP( (-4.57E-5)4P7) ) 

TR  =  i  a 

40  0  XXW= (3.23E-4) 4P04SR4 (1044  <-6.1E-5*ZS)  ) 

GO  TO  1 J  8  0 

C  COMPUTATIONS  FOR  PEFLECTEO  TRANSMISSION  (NEXT  18  LINES) 

8 o 0  CTHETA= (DZ-HL) /SRR 

IF (APS <CTHPT  A ) .LE..o31)  GO  TO  950 
IF (HP • LE. TH)  GO  TO  85' 

IF(PZ.LE.ZH)  GO  TO  825 

TH- EXP (-16. 4 /VIS4 <E4-Pi) )4EXP(-16.4/ ( VIS4CTHE TA ) 4 ( E4- El ) ) 
TR=EXD(-.i,3754(El-E3) ) 4FXP (-. Q975/CTHE7A4 (E1-E2) ) 

GO  TO  9o3 

825  TH=EXP (-lo.4/VIS4 (E4-E1) ?4EXP(-1&.4/ (VIS4CTHETA)4CE4-c2)  ) 
TR=EXP<-. J8754(E1-E3)  ) 

GO  TO  9 J 0 

850  IF(PZ.IE.7H)  GO  TO  875 

TH=EXP  (-16,4/VIS4  (E4-t  7 )  )4EX°(-16.  4/  (VIS*CTHETA  )  4(F4-E1)  ) 
TR=EXP(-. v  875/C  THETA4  (El -E?)  ) 


o  o 


GNE/PH/72-3 


GO  TO  9JQ 

875  TH=EXP(-16.A/VIS*  3 >  )  *EXP  (-16.  4/  ( VIS*CTHE  TA  )  * ( E4-E2)  ) 

TR=1« 

03  0  XXU=2.3*PO*A9S(10**(-F.iE-5*HL)-i3**{-6.1E-5*ZS)  )  f2.3*PC/CTHETA* 
l(l»j*M-6.iE-5*HL)-lG**  (-6dF-5*PZ)  ) 

GO  TO  1380 
950  TT=0 
RETURN 

FINO  TRANSMISSIVITY  FOR  INFRARED  (NEXT  11  LINES) 

10 8 C  IF (CFK . EQ «  2)  GO  TO  ljqo 
C  GROUND  °URST  CONDITION  (ASSUME  3 0 0 C  BLACK  BODY) 

FV=  •  1 
FIR=.9 

CALL  SETUP  (TPWl,l,2,12,0,0,G,u,0) 

CALL  MACURE  ( TBT1 ,  XXW ,  0,  3 , 0 , 0  ,  C  ,  I.ER ,TW) 

GO  TO  1°9F 

C  AIR  BURST  (ASSUME  600jPn) 

1090  CALL  SETUP  (T8W 1 , 1, 2, 1 2, C , 0 , 0 , 0 , 0) 

CALL  MACURE  (TBT2 , XXW > 0 , 0 , 0 , C , 0 , LER , TW ) 

FV= • 50 
FIR= . 50 

C  FIND  TOTAL  TRANSMISSIVITY 
1095  T7=FV*TR*TH+FIR*TW*(.7+.3*TH) 

ENO 

SUBROUTINE  MOTION (N,SR V, SR ,HZi , DT , SFR , FL AG , SE RD ,TT , HZ SA V , SRPD T , INT 
1C) 

COMMON/DOS/R, BURS TX, BURSTY, BURSTZ 
COMMON/ CRN D/HZ,HG,Hn,W 

COMMON/I N^O/POS X  (  HO)  ,POSY  (1  J  n  )  ,  VX  ( 1  u  "  )  ,  VY  (1 0  0 )  ,  VZ  ( 1C  3 i) 

COMMON/ BIN/POSX 1 , POSY1 
TT=TT  +  OT 

POSXl=POSX(N)+VX(N)*TT 
POSY 1= POSY (N) +VY (N) *TT 
HZ1=HZSAV  +  V7  (N) *TT 
IF  (HZl.Gt.HG)  GO  TO  IPO 
INT  C= 1 

130  SPSAV=SR»1jP0. 

OX= (POSX1-BURSTX) **2 
DY=(POSYi-qURSTY)4*2 
DZ=(HZ1-EURSTZ) 

SR-SODT(QX+CY+C7) 

C  FINO  PATE  AT  WHICH  SHOCK  FRONT  IS  CLOSING  CN  VEHICLE  TN  TIME  DT* 
SEPD=SRSAV-SFR 
SFo=SGV*TT 
SFPDT=SR-RFp 

OSEP=ABS (SERB) -AqS(SEnOT) 

IF (DSEP. GT , 0 . )  GO  TO  12P 
FLAG=FLAG+1. 

120  IF(SEPOT.oT.j.)  GO  TO  13C 
OT=-DT 1/10 « 

RF  TURN 

130  VCLOS-ABSOSEP) /AqS(OT) 

OT-SEPDT/VCLOS 

DT1=DT 

pnp 
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C 

C 

C 

C 


BLOCK  DATA 
********** 

BLOCK  DATA  CONTAINS  TABULATE!  VALUES  LSEO  IN  T»E  MAIN  PROGRAM 
AMO  IN  SUBROUTINES  TRIPNT , BLAST  1  *  TRAN$,  AND  THERM 

COMUON/GRNC/MZ,HG 

COMMON/T AB TH/  TBW 1  ( 12) , TPT 1 ( 12)  ,T9T2(12)  , 

1AAAK22)  ,AAA2(12  3)  ,B91(3Q)  ,BB2(38G)  ,TFV1 ( 8) ,T?V2 (3) 
COMMON/TRI/PZ,°G,FP,  PZR,FGR,PE9,  SR, FR  ,CRE, ALP 1ER , ALP1" , 

1  R90,PBPZ,R,PnoW 

COMMON/T A3/  TABU  (£9)  ,  TABIO  (69)  , TAB 2 1  (6?)  ,  TAB 20  (62)  ,TAP-*K69)  , 

1  TAB3Q  (69)  ,TAP£*I  (18)  ,TABAO  (18)  ,  TAB  51  (26)  ,TAB5D  (26)  , 

2  TAB6I (69) , TAB6D(69) 

COMMON/CON/UL2 (7)  ,bL3  (8) ,UL4 («) ,UL5 (5> , 

1  C4(7) ,C5(7) ,C6(7) ,C7(3) ,C8(8) ,C9( 8) ,01) (R ),C11(8) , 

2  C12(8),P2(8) ,P3(3) ,A1E(41) , ACF (13 ) , CF ( 13) 
nATfi  (TftPi  Tf  n  .  i  =  i  . p,qw 


****< 


33?.11B,36.5'53,4J.A47,^4.7j1,5C.0J3, 

100, JO/ 

DATA  (TABIC(J), 

J  =  i ,  69  > / 

1172CC  • ,  8 0  0  j  » , 

4240  .,  249  C . ,  1660., 

90  0  .  , 

5  4c., 

36. j.  , 

246.  , 

2  186.,  144,, 

1.7,,  88.8,  73.5, 

62.7, 

54.5, 

46.6, 

36.4, 

3  29.3,  24.1, 

20.2,  17.2,  14.8, 

12.9, 

11.3, 

9  .9, 

8.  •», 

4  7.9,  7.0, 

6.2,  *.2,  4.45, 

7 . 83  , 

3.7 :, 

2.  85, 

2.43, 

5  2.10,  1.85, 

1.62,  1.43,  1.27, 

1.12, 

.°8, 

.97, 

.77, 

6  .68,  .e"1. 

•  53>  .  U7 ,  •  4  i ) 

.36, 

.32, 

.  28  , 

.252, 

7  .222,  .193, 

•174*  .15$,  ,139, 

.124, 

.  113, 

.10,  , 

•  091, 

8  .083,  .,  755, 

•  -.69)  .  v  63)  •  .,574) 

. '333/ 

DATA  TAB2I/ 

1  u.,.1,.2, 

2  . 3C ,  .35,  .  4 ,  , 

•  45>  .5i) )  #55)  .$1, 

.f5,  .7 

%  .75, 

.  8  ’’ ,  .8 

~ ,  .9  , 

3  .95,  1.0,  1.5, 

2  •  £  )  2.5%  3oU)  3*5) 

4 . C ,  4. 

5,  5.0, 

5.5,  6 . 

0,  c. 5, 

4  7.0,  7.5,  8,^, 

8«5)  9  •  J  $  9*5)  1 J  • ) 

15.,  20 

25., 

30.,  35 

•  ,  40  . , 

5  45.,  5  ,  55,, 

6  )  65w)  7 j . )  75# ) 

8  u » ,  85 

.  ,  9  0., 

95.,  10.' 

•  I  1 5  0  »  , 

620 0», 258. ,300., 

353.  )4  JO.  )45v.  •  )5d./ 

DATA  T A 9 29/ 

1  .020,.  0625,  .125, 

2. 1800, .215 J, .245 1 , ,29 7C, .3323, .3650 

,  .39io,. 

41^3, .4 

350 ,.457 

: ,  ,46pj 

63.3,  72.)/ 
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1 

»  •  J  , 

u  .  1  , 

0.?, 

0.3,  0.4,  1.5, 

U 

.  6 , 

o.7,  <j»8, 

0.9, 

2 

1.  J, 

1.1, 

1.2, 

1.3,  1.4,  1.5, 

1 

•  6 , 

1.7,  1.8, 

i.% 

3 

2.0, 

2.1, 

2.2, 

2.3,  2.4,  2.5, 

4 

1.53, 

1.62, 

1.75, 

2.C2,  2.79,  2.9% 

% 

6% 

4.3%  5.70, 

5.46* 

6 

5.56, 

5.37, 

4.46, 

3.61,  3.38,  2.45, 

2. 

12, 

1.86,  1.65, 

1.5% 

6 

1.36, 

1.  ’5, 

1.16, 

1.08,  1  .k  3  ,  1 » u* / 

DATA  (TA96KJ), 

J=1 ,69)  / 

1 

.05.,  . 

’  625, 

.  75, 

.0875,  .130,  .125, 

. 

15C ,  .1*5, 

.20  . , 

2 

.226, 

.25  , 

.275, 

.700,  .325,  .35 

0, 

. 

375,  .4  .% 

.45  , 

7 

.  6  C  c, 

*  55  % 

.  6C  C  , 

.650,  .70%  .75% 

. 

8*0,  .860, 

.90  -  , 

4  .95%  1.0%,  1.1C5,  1.221,  1.35",  1.492,  1.64g,  1.322,  2.C14, 

5  2.226,  2»46.f  2.718,  3.0C4,  3. 72%  3.669,  4.755,  4.482,  4.957, 

6  5.474,  6.u5Jf  6.686,  7.389,  8.166,  9.C25,  9 . c74, 1 1 . 3 23 , 1? . 1 8 2 , 

713.464,14. 86- ,16.445, 18.174,29. €86,22. 199, 24. 52 2,27.113, 29. 96 4, 

833.115,36.  59%4j.*47,  ♦4.7Cl,5J.0u'',lJr.0U/ 

OAT  A  (TAR6C J)  ,  J*i,69)  / 

1.  PC  042,. 75 10 6,.  3  3 17%.  *3.255,.  0:  763  ,  .  j  C  63  8 ,  .  0  0  99  4 , .  0 14  32  ,  .  r  19c  4 , 

2. C2555,.33235».<i4"32».04Qj3,.*  5844,.')*85  3».w7924,.k.90  55,.11481» 
3.141b 5,.  169 -5,.  19  66  %  .  22955 ,.  26173  ,  .296  01,  .  32  °?9,  .  364  45  ,  .4  JC41, 

4.  437 1C  ,.47444,  .555,  .646,  .749,  .864,  .  994,  1.  14,  1.30, 

5  1.48,  1.6%  1.94,  2.15,  2.43,  2.73,  3. .7,  3.45,  3.86, 

6  4.32,  4.37,  5.4c,  6.c2,  6.71,  7.47,  3.72,  9.25,  1C. 29, 

7  11.43,  12.69,  14.(9,  15.64,  17. ’4,  16.23,  21.32,  23.62,  26.  17, 

«  28.9%  72.11,  36.55,  -*9.26,  44.1",  9".  %/ 

DATA  CD».2  ( J)  ,  J- 1,7)/''.  ",.35,.  45, 0.6, 1.7, 3.9,1.  ./, 

1  (UL3(J),J=l,9)/.J?5,.6,l.‘,2.5,7.t,10..,17.:,J.C/, 

?  (UL4(J)  ,J  =  1, 8 >/.2%C. 5,.  85,1.12,1.5,4.2, 17,.%%:/, 

3  (UL5 (J) ,J= 1,5)/. 265,, 47, 1.68, 5. 8 67375,u.O/ 

DATA  (C4 ( J ) , J=1 , 7) / 

1-. 7  246 16, -.0%-.  3  2,  .416667,-. . 71 25 4 , - . 0 ? 7659, -.0  31797/ 

0 AT  A ( C5 ( J)  ,J  =  1,7)  ,(C6(J) , J  =  1 , 7 ) , (C7 ( J ) , J  =  1 , 8 ) , ( C3 ( J ) , J= 1 , 8 ) / 
1.087493,-.  11 C ,  .13%-. 72 J* 27, .277338,  .  80934, . c29527, 

2. 167 C66, .166  2, . J724,.  145  ,.o21749, ,173  847, .2724  44, 
3«16PS75,-.6C6r61, 1.44,. 097778,. uT8552»0»C,-.  •  00 3 7 8 7 , 7 6. 202 , 

41.8138  69,-3. 554546,-7. 16, -.5422 23, -.11 42^3, -r.c  3167,0 .0  ,".7/ 

DATA  (C  9  <  J )  ,  J  -  1 , 9  )  ,CC1.'  (J)  ,J=1,8)  ,  ( C 1 1  f  J )  ,  J  =  l,  6),(C12(J),J  =  1,8)/ 
12.0443, 3, 928788, 3. 22,1. 944445, 1.43273'1, 1,06465 °,1.  0  56  28  5,1.  C, 
21.153846, -6. 2 97 2, -1.3 7 14 3, -3.05556, .6 1 8421 , . U5 ,667 ,- .C 01" 8276 , -3. 

3.  -.  0  0  3346, 7. 9?  t,  54, 3. 0  3  571, 6. 0194  4,-1.8  0976,-.  :  654667*  .0  224  069 , 0  . 0  , 
4. 0146514,- 1.6 379 71,-. 345, -1.37 8 89, 2. 35 11 9, 1.1 14 3,. 854207,1.0/ 

HAT A  (P2( J) ,J  =  i,8) ,  (P3(J) ,J  =  1,8)  / 

1  5* 2 • , 1 . , 3. , -3 .  ,  7*2., -2./ 

OATA  (A1E(J) ,J=1,41)/ 

1  40.  C,3y. 6, 79. 4, 3 9. 2,  S9. 7,39. %  39. 0,39.  0,39.1  ,39. 3 , 39 .5 , 39 . 6 , 39 . 9 , 
?  40. 0,4 0.4, 4  -.7, 4 1.0, 4 1.3, 41. 6, 42, 0,4". 4, 42. 7, 4 3.7 ,43.6,44. 0,44.5, 
*  3  45. 3, 46. 4?. %  4 9. 7, 49. 5, 5C.  3, 52.  3, 5 4. 0,56.0  ,5  8.2, 6u. 5, 67. 5, 67. 5, 

4  74.0,90.0/ 

0 AT  A (CF ( J)  >J  =  1,13)/1. 0,. 98, .96, ,917  ,.93, .66, .47 ,.343,  .265, 

1  .211, .  175  , .  143,0./ 

DATA (ACF(J) , J= 1,13) /25 00 0  ., 50000 . , 625 JG. ,750 j " . , 875 0 0  . ,  1 0 0  %  0  .  , 

1  1125"  0.  ,  I25.no.,  1375  C  . ,  150  JOi  .,1625  CO.,  175  j  0-.,  25  «0  JO./ 

DATA  (TBWt (J) , J=l,12)  / 

10. ,17  ,,20.  ,3j»,4Q.,59»,60»,  8*2  ,  ,100.  ,120.,  140.  ,150./ 

DATA  (T3T1 (J) ,J  =  1 ,12) / 

1. 89,  .  72, .  65,.  62,.  58,  *  r>  6  ,  «  54,  .51,  .48,  .45,  .42,  .41/ 
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DATA  (T9T2 (J)  ,  J=l,12) / 

!  .97, . 89,. 74, .7j,. 67, .64, .62, .58, .54,. 51, .48, .4 7/ 

DATA  (AAAI (J) ,J=1,2?) / 

1  2.0C,uu3E+0j  ,  3. OOGOOr+CO  ,  5.f  3GCCE*0'',  7«5j0jGE  +  00,  1 . 0  0  o  ■:  0  E  +  r»  1 , 

2  1.5CC0 CE* jI,  2.9JCC3E  +  ul,  2 . 5  C  u  ?  .  E  ♦  C  1 ,  3.  -  ’J  0  CL  E  +  °  1 ,  4.C00‘,CE  +  'j1, 

7  5.900 Out^Jly  1. 3  ;G )3r*"2, 


4  (\,  .17,. 35, .52,. 7C,. 37, 1. 95, 1.22, 1.4 :, 1.57/ 
HATA  (AAA?  (J)  , J=1,6C)  / 


1 

9. 588006-0  2, 

i.  3 :e?or-o 1 , 

1. 1318CE-31, 

1 . 29  j  7C  E-0 1 , 

1. 437 if E-3 1 , 

2 

1.5688QF-C1 , 

1.63P?3F-j1, 

1.5°633F-G1» 

1.384  C. E-31 , 

1.15j3CF-01» 

3 

1.98.3  UE- J1 , 

1.9957 OE- 01, 

1 • 9837“ £- a  1 , 

1 • 955  2C  E-3 1 , 

1. 9049:6-01, 

4 

1 . 8  C 120E- 31 1 

1.6  -833r-ul , 

1  •  3 115CE- 0 1 , 

9.597  L  E - j 2  , 

7.330-.  Jt-no, 

6 

3.07139E-01, 

3.0  j440r“jl, 

2.9163CE-01, 

2 . 559  ft  E-0 1 , 

2.24743F-01, 

6 

1 . 877  4  CE- 31 , 

1. 45613E-.1, 

l.?142'*E-Cl, 

6.25o::e-o?, 

4.302^5-0?, 

7 

?.68833r-Jl , 

3.574,]F-C1, 

3.2843  . F- 0 1 , 

2.8833  E-'l, 

2.4064CE-01, 

8 

1.882 COE- 01, 

l.  748l  Jr-^4.» 

8.4873.E  —  Jc, 

4.591  -i.E-32, 

2.800n6E-02, 

4.CC93JE-31 , 

3.9726 Oi-Cl, 

3.5  30  2  )F- 31, 

3.053  i C E - J 1 , 

2.4974CE-31, 

1. 87773 E- 01 , 

1.2°.  52JF-V.1, 

7,622005-02, 

3.7560. E- 32, 

2.1G3 jCE-02, 

: 

4, 3362 JE- 31 , 

4 . 1 7°40E-3 1 , 

3.7834JE-C1, 

3.21971 E-31, 

?.5587r.E-ul, 

7 

1  .87.3u--31, 

1.2196 J"-01, 

6.72C3CF-02, 

2.918CCE-02, 

1.4 03  COE- 0  2/ 

DATA  ( A  A  A  2 ( J  >  ,Js61,12  )  / 

i  4.50133E-31,  4.336 .15-91,  3.9l?5.E-01,  3. 3376l E-Jl ,  2.596206-31, 
?  1.8637-c-Ol,  1  •  1 94 8  '*r-u  1 ,  6.2740 9 F-02*  2.493\E-J2,  1, lOT 0vE-?2, 

3  4,60>.30F-'J1,  4.43C°3r-Cl,  3.99063E-C1,  3. 35  7  V.  E-3  l ,  2.6197:6-01, 

4  1. 859906-31,  1.163636-01,  5.9980-F-3?,  2.242.1 E-02,  5.7CCrCE-03, 

5  4. 66723 £-31 ,  4 . 49453E-0 1 ,  4  .  C 43 7C  E- u  1 ,  3 •  393  C  ■„  F- 0 1 ,  2. 63341*  c-9  i  , 

6  1 »  8572  j£-  31  *  1.14°43£-ni,  5.8l20*E-32,  2.*73Cc E-32 ,  7.  ?■  03  :  3E-0  7, 

7  4.751-1c-j1  ,  4.57513.-01,  4.1C96CE-01,  3.4336C F-o 1 ,  2.65210E-01, 

8  1.8534wE-'Jl »  1.131236-^1,  5  .c7833E-02,  1.  3613.E-02,  5.3  04  CE-03, 

9  4. 80.93  £-01,  4.62c2J“-,,l,  4.1493-E-31,  3 . 462  9.  E- 3  i  ,  2.6632LE-ul, 
J  1.851232-01,  1  *  1 2C30E-0 1 ,  5 . 4  ?6u  3P- G  2  ,  1.732:CE-32,  4 .  J  0  0  0  C  p -C  3 , 
S  4.9C  01CE-01 ,  4.7i93uf-91,  4.22860E-01,  3.515HE-31,  2.6949LE-G1, 
?  1.8456GE-J1,  1 • J  98  3  3c-w 1 ,  5.1543.E-02,  1.476-CE-02,  2.03^06E-u?/ 

OATA  (  3 *3 1  (  J)  ,  J=l,39>  / 

1  0.  ,  1.0>CC0E*oG,  2.r,:C0j£40G,  3.  00 0 . C E*3 9  ,  4. 0  0  j  j  .  E  +  0  u  , 

2  5.0  0(.vJ3p>OG,  6.  9rOCdr*-ll,  8 . 9  •  0- QE  ♦  0  ; ,  1 . - 0 j CC E*0 1 ,  1 .  2  2  '£  •*  C  E  ♦  0 1 , 

3  1.4UO  Ju->91,  1.6CCoOF  +  ?l*  1 . 8 >' 0 3 CE ♦  0 1 ,  2 .  C  3 0  Cl £  +  9 1 ,  2.  20 4 p  CE ♦  C  i , 

4  2. 6'J  .  OGE  +  01',  3.0  jCHOFtCl,  3.5.33CE  +  31,  4.000  C3E  +  01, 

5  0.  ,  2.5GC . 3c-C 1,  5.j'JO0E-O1,  7.5’GCrF-Jl,  1. 0 00 OC EO '  , 

6  1 . 5  0  j  j  1 E  ♦  0  0  ,  2. 0-C  G3F>CC  ,  3. 0  CGOOE  ♦  0  C  ,  4.3jCCwE*3C,  6.03C  j  C  E  <■  0  2 , 

7  e.OwCOOE+Ot,  1.3GC0jE401,  1.2CG3lE*01,  1.4'3  0':E+31,  1.80C  CCF  +  01, 
9  2.2C0  3oE+0l,  2.6CCCOF  ♦  l»  1  »  3.CG030E+01,  3.5n0  JLE  +  Ol ,  4.  G 3 G - GF *9 1/ 

OATA  ( 992 ( J) , J=l,95> / 

1  i.CC.GOEO:,  4.943rJ"-01,  1.6430CE-31,  4.29OC:E-02,  1.  T3r 0 C E-0 2, 

2  E^I^CjE-T?,  1. 994036-.  1,  4.534  )36-01,  6 . 67  0  0 .  E-3 1 »  9.72GCGE-01, 

3  1.1 67  GOF+Ow  ,  1.31C“JE*00,  1 . 41 G  3  u  E  *■  0  C  ,  1.495,.E+30,  1. 594C  CE*0  :  , 

4  1.6693  J  E  +  J  0  ,  1.722..  JE  +  jC,  1 . 76  40 .  E  ♦  0  .  ,  1.31"..E  +  00  ,  1. 835JC E*Q ?, 

5  7.  C7C  30E-31  ,  5.3753JE-C1,  4.1963CE-01,  3.  334 ;.C  E- j  1 ,  2.6930GE-01, 

6  2.6090 JE-01 ,  T.11CCJE-01,  5 . L 34  )  uE- 0 1 ,  6. 861  ,r E-91 ,  9. 79G3 C E-0 1 , 

7  1 . 16.  0  OE  +  ou  ,  1.30Gr0r40.,  1 . 400  OOE*  0  ,  1.48^  '-E  +  OC,  1.592.CE*'.", 

8  1.66830E+0C,  1.72CCOr  +  OC,  1 . 7620  0E  + 0  r  ,  i .806 CC E+0 0 ,  1.  8 32 j C F ♦ o * , 

9  4.47.  3GE-31 ,  4.314035-31,  4.2680  3E-C1,  4.3.2L:E-0i,  4. 372”' OE-C 1, 

*  4.4803CE-31,  4.724.0C-01,  4.571rCE-01,  7.206!CE-31,  9.633006-01, 
5  1.14O03E+  00  ,  1.2815 9 E+QG,  1.78200E*CC ,  1.46?.rE+33,  1.5820uE  +  ut, 

t  1.66t'3uE  +  3C  ,  1.712C0E+00,  1.7563uE*0C,  1. 799  E  +  00  ,  1 .  9  23?  CE  +  0  0  , 
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c 

3.16 -QCE-31, 

3 . 226  C  05-0 1 , 

3.4353:5-01, 

3.724GCE-J1, 

4.U040CF-01, 

X 

4.5670  iE-Jl, 

5. 35-01, 

6. 171Q05-01, 

7.310 .TE-Ol  , 

9.  425' lE-p 1 , 

r 

1  «1?.313E+0l  y 

1.258C0E+00, 

1 • 35C  0" E  +  C 1 , 

1.44 . or  5  +  0 3  , 

1.566 '65  +  0  0, 

<*■ 

1.6460  J5  +  iO, 

1. 7  ;2C  3r  +  Qr  , 

1.7493CE+ j J , 

1.7a*.  CcE  +  JO, 

1.8190 C5+g  .  , 

* 

2.47  005-11, 

2.49C-  3l-:i, 

2 . 7* G  3  3E- C 1 , 

2.96Q:CE-j1, 

3. 290005-01, 

* 

?.967oOc- jl. 

4.  6?5C‘3e-Cl, 

5. 94S30E-C1, 

7.  ]37CrE-gl, 

9. C820CE-01, 

♦ 

1  .  086005+03, 

l.l^lwJF+ro, 

1. 3260CE  +  0C, 

1.41.CCE+3J, 

1.542005+CC/ 

niT  a  (BH2(J),J 

=  a  6 , 1 9  3 )  / 

1 

1.64  iqif+g:  , 

1.68€jOE  +  ij, 

1.737305+00, 

1.78. "CE+g7 , 

1.  R10o05  +  C", 

2 

1. 962Q0E-11, 

2.  J  cnc-gi. 

2.17731^-31, 

2.36.'  O'E-Ol, 

2.672"0E-31, 

3 

2.  266j  15-11, 

3.957  05-ul, 

5.2630  .5-31, 

6.35  .  j'.F-^I  , 

8.60**0PE-C1, 

4 

1.0380j5+Cl, 

l.lHCOEKC, 

i.2°2u:>e+c: , 

1.381  : C E ♦ o  0 , 

1.529 005+00, 

5 

1 » 6 1  v.  0  u  5 +  Jg  , 

1 . 6  721 0  5  +  CC  , 

1.722 3  . E  +  O:, 

1. 77C  C..  E  +  00  , 

1.795«0E+ jO, 

6 

1 .64h  J  jE- J.r, 

1 .71LCJE-tl, 

1.9C60 jE-01, 

2.CJC  r 5- 3 1 , 

2.20100E-01, 

7 

2.7  37  Oat- 01 

3  ,  3440Jr-"l, 

4.65-g£-31, 

5.722Q.E-J1, 

8.C0?'0E-0i, 

3 

9.8110 jE- Jl, 

\ .  127CJC  +  'J0, 

1.2410CE  +  0  ’  , 

1.741..C5  +  0  j, 

1.  486.-0E  +  09, 

9 

1. 57911E+9G, 

j.  6  49C  j 5  ♦  0'j  » 

1.7-6C05+0C, 

1. 753  ■*  E+jO, 

1.794,  CE  +  C", 

* 

1.24nu'JE-Cl, 

i.'c  30  3  05-01, 

1  •  3  4  C  3  0  e  -  g  1 , 

1.46  3  T-E-01, 

1. 599DC5-C 1 , 

1  .96  0 jOt- Jl, 

2 . 429C35--1  , 

3.4460 jE-01, 

4.55.juE-j1, 

6.6740CE-D1, 

t 

8.527-15- J! , 

1  .  o  It*  J-I+OQ  , 

1.133u1E+0l, 

1.246  ,  *•  E ♦  0  3  , 

1.  4o9C:E+0w, 

♦ 

S 

1 . 5  ?  C  0  u  E  ♦  3  j  , 

1.60CCJ£»00* 

1.6630  E+CC, 

1.721*'c5  +  30, 

1.76C-.0E  +  3C  , 

<+. 

9.9C00  05-02, 

1  .  ICC  ■  0  L  * '  0  1 » 

1.c65jCE-G1, 

i.i5:plE-ji, 

1.217  iCF-Ol, 

1.485jSE-J1, 

1  .  3  2Cgl--0l, 

2.62'.  0"E-G1* 

3.536':5-gl, 

5.457  :  9  E - 0 1 , 

* 

7. ?69u3E- 01  , 

3. 341C  J r - 0 1 , 

1.  j22C.£  +  0C, 

1.143  /  E+30, 

1.  321CGF+U**, 

C* 

1.4511JU+3:, 

1.647*.  95  +  Ct, 

i.ft2i;:E+o*, 

1 » 6  9  5  3 .  E  ♦  0  ^  , 

1.735v0E+0j, 

<** 

> 

8.?c::  JCc- J2, 

3. 5  JC  ;  J5-G2, 

8.9.00  jE  —  uii, 

9.i:  j  'C  5-  3  2  , 

l.OOC  jgE-Ol, 

* 

1.17.  j3l-.i1, 

1.43fC3E-Cl, 

2. CSC JC5-01, 

2.78  .  .wE-Jl, 

H.41j-i.E-U1/ 

OATA  <U02<J),„ 

=191,295)7 

1 

6. 1C. 0 1 E  - 1 1 , 

7.*9C0Or-Jl, 

9.  :Pr CCE-01, 

1.028fE  +  0'3  , 

1. 2290CE  +  0  . , 

2 

1.37. JGE+1G, 

1.478:nc*9C, 

1.5620'jE  +  Or, 

i.649  :;c+rfo , 

1.781  IE  +  -J 

r 

7. 12-C0E-G2, 

7. 2  3G-1E-U2, 

7 . e61 C 1 F- 0  c  , 

8.30" *Cc-C2, 

8.  460  005-02, 

4 

9.32  0jQc.-02, 

1 .166o15-31, 

1.6470  3 5 - C 1 , 

2. 249  1" E- J 1 , 

3.6640 g£-Cl, 

5 

5.1733 05-01, 

6. 554COr-91, 

8.CCC3CE-01, 

9.244CCE-31, 

l.lll-CF+Oj, 

6 

1.288  JOE*  V.  , 

1  .4*10  35 ♦0  0, 

1.5C3CCE  +  U  C, 

1.6j;Cic>Qlf 

1,6630«E,03, 

7 

6.24'-  03E-32, 

6.3CC- JE-02, 

6.5ow’315-C2, 

6.9*  0:l'5-:2, 

7,20  ..uE-02, 

9 

8.4Cg0O£-O2, 

9.9  ’C03E-C2, 

1.35CPQE-01, 

1 . 8  ?  g  .  ■  E  —  g  1  , 

3. 0  79  OC  E-01 , 

9 

4. 39 : 0  0  5-31, 

5. 7  1C  0  0  5-01, 

’..6  JC  OE-Ol, 

8.30-  'wE-01, 

1 .  0  4  C  0  C  F  •?  J  C  , 

<? 

1.2r31uE+30, 

1 .329lJ5  +  00, 

1.432005+Qr, 

1.543.-  5  +  3  i  , 

1. 617-05+0", 

* 

5. 55  CO- E- 02, 

5.6:CC3C-G2, 

5.67609E-02, 

6.G0CCCE-32, 

6.43CQOE-02, 

<* 

T> 

7.2’w  j  o«=:- 0  2 , 

9.33c:ot-r.2, 

1.153-CE-J1, 

1.556CC E-01, 

2»63u9CE-Cl, 

? 

3.745C IE-01 , 

4.98C03r-Pl, 

6.2*00CE-0i, 

7.383  * 05-0 1 , 

9.46730E-J1, 

9 

1. 12C OiE+OC , 

1.252G3e>Gg, 

1.371G-E+Q0, 

1.481  3P+30, 

1.  5  67  j  C  E  +  0  0 , 

.« 

4,98  0  OE- g  2 , 

5  .  3  10  Cl £-0  2, 

5  •  ».5  J-l  E-  0  2  » 

5.20CCL5-C2, 

5.65 fj  05-02, 

a 

6.9CC  Q.3E-02, 

7.25u03E-j2, 

9.9CCQ05-0Z, 

1.310  QGE-31 , 

2. 19090E-01, 

* 

■’.21C00E-31, 

4.7?cnE-jl, 

5.50GO0E-01* 

6.56-3  ] t E - -3 1 , 

3.61P00E-J1, 

<  ? 

i.G4CQOc+ OC, 

l.moOF  +  Cl, 

l.TUCgE  +  0w, 

1.437 j.E+OO, 

1.5160OF+ j‘, 

% 

4.54000E-12, 

4,540005-02, 

4.54jOgE-02, 

4.  930  CPE-02 , 

5. C30  OCE-02/ 

OATA  ( R  9  2  ( J )  »,  J  =  2  8  6  »  3  9  ')/ 

1 

5.65.035-02, 

6.460306-02, 

8.62-Q0E-Q2, 

1.148:^5-01, 

1.895JCE-01, 

2 

2. 7990 OE- 01 , 

3.9  12C35-01, 

4.88Q3CE-31, 

5.898 :CE-11, 

7.  8  74f 05-^1, 

3 

9,6140 JE-01, 

i.ii8CJ'430, 

1.234  G-jE  +  SC, 

1.359  .15+30, 

1.463  iuE  +  G 

4 

3.84  00E-32, 

3.900005-02, 

4,0 1  uO.E-Oc, 

4.i:o'.E-02, 

4.30r-0E-J2, 

5 

4.6CD00E- 32, 

5. 6uC00£-02, 

6.7CC03E-02, 

8 . 9i>g  Ci  E-J  2  , 

1.421’?  QE-Oi, 

6 

2. 18C JoE-31, 

e.ggcGD^-oi, 

3,980005-01, 

4. 820  g" E-01 , 

6.520 JoF-01, 

7 

8.263005-11, 

9.97C j0c-01, 

1.10910E+GC, 

1. 238  Cl  E  +  03 , 

1.352:CE+C j, 

94 


oooooooooooooooooooooo 
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8 

3.330JDE-02, 

3.301037-02, 

3.30 JQ0E-02, 

3.400  :CE-02, 

3.57" OGE-02, 

9 

3.9300GE-Q2, 

4.37CCQ--D2, 

5.56C‘8jE-Q2, 

7. 170  0t.E-02, 

1*152. OCE-91, 

* 

1.71400c- JI» 

2.3  723  Jr- ) 1 , 

3,  :SG00E-G1, 

3.891 3lE-01  , 

5.4*9G.E-C1, 

7.042008-01, 

8»52CU~-',1, 

9.ri5o:e-oi, 

1 . 123  E  +  J  o , 

1.24100E+0-, 

<t 

2.86C00E-G2, 

2.8 jf C0r-)2, 

2.6CC0CE-G2, 

2.80  0  TE-J2, 

2.90- 0CE-C2, 

* 

7.30.0)8-02, 

3. 6;/ COO  r- 02, 

4ol.o:^E-02, 

5.8jr’  GCE-02, 

9 .  C  G  0  i  0  E  -  0  2 , 

■*» 

1 .32000E-G1, 

1.9  C  0  D  F-3 1 , 

2.52C0',E-01, 

3.18-j  CC E-Dip 

4. 57C0QE-U1, 

<5 

6.03. JOE- 31, 

7.43CQJ--C1, 

8.480?  GE-0 1 , 

9.96G  . . E-Gl , 

i.  113  0E  +  0  ., 

■b 

2 . 5  C . 0 3E-02, 

2 . 5  u  C  C*  0  E  -  J  2  , 

2. 50 J  Q0E-02, 

2.5:jc;£-j2> 

2.610:0E-P2, 

5 

2. 82ujoE—  D2, 

3,  08CCJ--02, 

3.7FCG:E-0?, 

4.700  CL E-j2, 

7.24)  ,  0  E-G  2, 

S 

1  .C640-1E-91, 

1.477008-01, 

i.92~3GE-Gl, 

2.4/9GCE-01, 

3.630:0E-01, 

? 

4.8550 JE-ul, 

6.G9GCu~-01, 

7. 293 JuE-01 , 

8. 7’.  Ol/uE— 01, 

1.  GOlGoE  +  O-./ 

OATA  (T8VKJ) 

,  J  =  .1 , 8  )  / 

10.,4f>'.'1C,,8CD.O.vi2C*0  .0 
DATA  (TBV2  (J) ,J  =  1,8)  / 
1»C8,»11,«20,.45,.51,«49 
END 


«  ,14G3GCi.  ,160060. ,  20 0  j  DC  *  ,240300./ 
,  .  4  2  »  .  3  2  / 


SUBROUTINE  ATMOS  ( 2 , TM , SIG^A  ,c>ho ,THET A , DELTA ,  CA  ,  A mu , K) 
CALLING  SEQUENCE 


CALL  ATMOS  ( 7. ,  TM ,  S IGMA  -  RHO  ,  THET  A  ,  BELT  A  ,  C  A  ,  AMU ,  K) 


Z  =  GEOMETRIC  ALTITUDE  CRT) 

TM  =  MOLECULAR  SCALE  TEMPERATURE  (DEGREES  RANKIN) 

SIGMA  =  RATIO  OF  DENSITY  TO  THAT  AT  SEA  LFVCL 
RHO  =  DENSITY  L 3-SEC** 2-FT** ( -4 )  CR  SLUGS-FT**3 
THETA  =  RATIO  OF  TFMP- RATURC  TO  THAT  AT  SEA  LEVEL 
OFLTA  -  RATIO  OF  PRESSURE  TO  THAT  AT  SEA  LEVEL 
CA  =  SPEED  OF  SOUN n  C^I/SER) 

AMU  =  VISCOSITY  COEFFICIENT  (L8-SEC-FT**2) 

K  =  1  NORMAL, 

=  2  ALTITUDE  GREATER  THAN  30  00  3(3  .  FT., 

=  3  ALTITUOF  NEGATIVE, 

=  4  FLOATING  POINT  OVERFLOW, 

=  5  ALTITUDE  GREATER  THAN  3 ?  C  0  0  C •  FT.  AND  FLOATING  FOINT  OVERFL. 


ATMOS  IS  CALLED  DY  ROUTINES  9LAST1  ANC  THCRM 


■  i  T- 


DIMENSION  HCRIM  ?  (11) ,  TMB  ( li) ,SIGMA9 (11), AIM (11) 
DATA(HPRIMeCX) ,1-1,11) / 


1 

2 
3 


■j  .  , 

173634. 510 , 
524934.380, 


36.'  39 . 239  > 
259186.3 ‘30, 
557742,790 , 


82020.997, 
25c275, 590  , 
656.167,980/ 


15^199.489 , 
344468* 190, 
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OATA(TMB(I) ,1=1,11)/ 

1  518.688,  3«9.938,  389.983,  508.738, 

?  508.788,  298.133,  293.188,  406.188  , 

3  2386.188,  2566.138,  2836.188/ 

DATA (SIGtfAB ( I) ,1=1,11) / 

1  l.CrOoOTOE  9J,  2.9706953E-:.l,  3.2665751£-t,2,  1. 211787GE-03, 

2  5.  8677311E-C  4  ,  5 .36  77  311E- °4 ,  1 . 76 385 95E-0  6 9.3921519^-0  3, 

3  7.7658593F-13 ,  5.S324877E-10 »  2.5?2b77lE"lG/ 

DATA (ALM(I ) , 1=1,11)/ 

1  -0,93356616,  3.  ,  0.00164592,  C,  , 

2  -0.. 9246883,  0.  ,  Q.0U219466,  0. 31.97280, 

3  Ci  •  0  J54864  J  ,  J.C227432C,  C.  33192024/ 

DATA  Q  /  0.318744176  /  ,  RE  /  2.9855531E  97  /, 

1  S  /  198.72  /  ,  PZ  /  2116.2  /, 

2  AHUZ  /  3.7372993E-C7/  ,RHOZ  /  0-,3f'23769  /, 

3  TMZ  /  518.688  / 

K=i 

IF  (Z)  IQ,  3u ,23 
10  K=3 

GO  TO  113 

29  IF  (Z.GT.73C2  10 . 

33  HPRIM=  (RE/ (RE4--".  .<  •  . 

00  4Q  M=l, 11 

IF  (HPRIM-HP=>i?-’3(K)  )  50,60,43 
40  CONTINUE 

M  =  12 

5  0  H-.-M-l 

69  IF  CALM  tM) )  7G,83,7u 
7  9  TM=TM3(H)  +ALN(M) * (HPRT M-HP^I M  0 ( M) ) 

SIGMA  =  FXP(  (l.r  +  (Q/ALM(M))  )  MALOG(TM°  (H )  /  TN  )  )  )  VSIGMA9  ( M ) 

GO  TO  90 
83  TM=TMBM) 

SIGM A  =  SIGM A3 ( M) *EXP<- <0* (HPRIH-HFRI6E  O' ) ) ) /TM P ( M ) ) 

90  RHO=RHOZ*S IGVA 

THET  A=TM/T  MZ 
oelta=sigma*theta 

(3 A  =  4 9 ,  j2177*S03T  (  TH) 

/iMU=AMU7*$QRT  (THETA**7)*  (  (  TMZ  +  S )  /  <  T  m  +  S ) ) 

110  i-^TUPN 
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SU9P0UTINF  MAGU-  (?,  XA1,XA2,XA3,XA4,XA5,XAF,  IF,7R) 

SUHOOUTINE  M  ACUR  2  EXECUTE’S  AN  N  iJlM^NSION  AL  TAOLE  LOOK  UP 
WITH  PXTPAFOLA  riCN  IF  0FS Io<r0 

CALLING  Sc  CUENCt- 


CALL  MAGURr(Z,XAl ,XA2, \A3, XA4, XA5,XA6,IF,7R) 

WHERE 

IE  =  ERROR  COOE 

0  INTERPOLATION  SUCCESSFUL 

-1  OP  1  ARGUMENT  EXCEED  LIMITS  OF  THE  T  A9L  c 

2  tnoefenoent  variables  not  in  ascencing  or der 


2< 

1)  = 

F(X1,Y1,Z1) 

Z<1?)= 

F(X3,Y1,Z1) 

Z( 

?)  - 

F(Xi,Yl,72) 

7(14)= 

F(X3,Y1,Z2) 

7  ( 

3)  = 

FvXl,Y2,7i) 

Z  <15 )  = 

F(X3, Y2, 71) 

?( 

4)  = 

F(X1, Y2,Z2) 

2(16)= 

C(X3,  Y2, 72) 

Z  ( 

5)  = 

F(Xi, Y3,?l) 

2  (  1 T  )  = 

F(X3,Y?f 71) 

Z  ( 

F )  = 

F  (XI , Y3 , 72) 

Z ( 19) = 

F ( X3» Y 3 »  7? ) 

?( 

->)- 

F(X2,Y1,?1) 

7 ( IP) = 

F  (  X4,  Y 1 , 71 ) 

2  ( 

8)  = 

F (X2, Y1 ,Z2> 

Z (2  0 )  = 

F(X4,Y1, 7?) 

1  ( 

P)  = 

F  <  X? , Y2 »  7 1 ) 

Z (21 ) = 

F(X4,Y2,71) 

7(1 

,C>  = 

F (X2, Y2,?2) 

Z  ( 2?)  = 

F( X4, Y?, 72) 

7(11)  = 

F<X2,Y3,71) 

2(23)= 

F ( X4 , Y  7 , Zi ) 

2  (12)  a 

F(X2,Y3,22) 

Z ( 24 ) = 

F(X4,Y3,Z2) 

**,**»»»*<•  »*,»»»»**, **«,**»*9«»»»«*****»»»*****»*»»**»W**V****,*  + 

nOMMON/TPLKUO/H.LFjMa  <F)  ,X  (lcO)  ,  NNEX 

OTMENSION  Id),  X  A  ( 6)  ,  N  S  ( 5 )  ,  KJ(32),  FATI0<5>,  KGP0UP(5), 

UTOTtS) 

IE  =  0 

NEXTRrNNEX 
V  A ( 1 )  =  X  A 1 
XA(2)=XA2 
XA (3)=XA3 
XA  (4)=XA4 
XA(5)=XAS 
XA (6 ) -XA6 
00  10c  1  =  1, LF 
L2=L1+N4 (I )-? 

.  FOUND=0. 

00  SO  J  =  L1  , L 2 

5  TF  ( X  ( J)  .GT.X(J-t))  GO  TO  10 
TE=2 
RETURN 

13  IF  (FOUNO.  NE.C-.1  GC  TO  50 
IF  CXA(T)  —  X( J-l) )  20,50,50 
21  IF  (J.GT.Ll)  GO  to  43 

IF  (NEXTR. EQ. r )  bO  TO  30 

F0UNP=1. 

NS ( I ) =L1-1 
GO  TO  5 J 
30  I E  =  - 1 

RETURN 

43  FOUNO= 1 • 
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ns ( r ) = j-p 

CONTINUE 

IF  (FOUND)  9  j  ,  6  J , PC 
f  r  IF  (XA  (I)  -  X  (l?)  )  -*0,8  .  ,70 

7,  IF  (Nt-XTR.  NF.-)  GC  TO  80 

IE  =  1 
PfTUPN 

8  0  NF  (I)  =i_2-l 

DO  Ll=L2+2 

10'  CONTINUE 

C  IN'  NS  ( I )  IS  THE  SUPSCPIPT  IN  THE  ARRAY  X  SUCH  THAT 

C  X ( NS ( T ) )  IS  LESS  THAN  THE  ITH  ARGUMENT 

nn  Ho  1  =  1  ,lf 
K=NF  (T) 

RATIO (I)  =  ( XA (I) -X  (X) ) /  (X (K+l) -X  (X) ) 

C  IN  PATIO(l)  IS  THE  RATIO  OF  X  ARG,  R AT  10 ( 2 > =R AT  10  OF  Y  ETC. 

110  CONTINUE 

NGDOUP ( 1 )  =  NS  (1) 

NSUM=MA(1) 

DO  12.  1  =  2, LF 
NGROUP(I)=NS(I)-NSUM 
NSUM=NSUM+NA ( T) 

12'  CONTINUE 

C  IN  NGROUP(I)  IS  THF  SU3SCRIPT  OF  THF  TTH  VARI APLE  SUCH 

C  THAT  THE  TABLE  VALUF  IS  LFSS  THAN  THE  CORRESPONDING  ARGUMENT 

C  THIS  IS  IN  T coMS  OF  MIS  VARIABLE  ONLY 

C  FOR  A  FUNCTION  OF  OEG^EE  NO  WE  NEE02 *  * ( NO-1 )  VALUER 

C  FpOM  THE  Z  ARPAY 

ITCT  <LF)=1 
00  13..  1  =  2, LF 
J=LF-I*i 

ITOT<J)  =  ITCT  <JH)*NA(JU> 
in  CONTINUE 

C  IN  ITOT(J)  IS  THE  NUMBER  OF  LOCATIONS  IN  THE  Z  ARRAY  NEEOEO  TC  CM 

C  THE  JTH  SUBSC°IPT 

KF=2**LF 
M  W  =  -  2 

00  170  1  =  1  ,  KF, 2 
IFIRSTrl 
M W  =  M W  +  2 
00  16U  J  =  1  , L F 
MM=2** (J-l) 

IF  (ANO(MM ,HW) , EQ.C.)  GO  TO  14C 

IMON=NGROUP(J)fl 

GO  TO  1ST 

140  IMON=NGROUP ( J) 

ISO  IFIRST  =  IFIRST*(IM0N-1>  *ITOT( J) 

160  CONTINUE 

ISEC=IFIRST>ITOT(l) 

WJ(I)=Z(IFIRST) 

WJ(I41)=Z( ISEC) 

170  CONTINUE 

00  180  1  =  1, LF 

KF=KF/2 

OO  180  J=i  , KF 
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180  WJ(J)=WJ(2*J-1) +(WJ(2*J)-HJ(2*J-i) )  *RATIO(I) 

7R=W J ( 1) 

RFTUPN 

END 


k 


w,  * 


c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  SETUP  <  X  ,NEXTP  ,NO  ,NAi ,NA2  ,  N A3  ,  NA4,  KAUNAS) 
SUBROUTINE  SETUP  SETS  UP  ARRAYS  FOR  TABLE  LOOK  UR 
CALLING  SEQUENCE- 

CALL  SETUP  (X,NEXTR,NO,NAi,NA2,NA3,NA4  ,NA5,NAG) 

WHERE 

X  =  TABLE  OF  INDEPENDENT  VARIABLES 


c 

c 

NEXTR 

— 

0 

i 

NO  EXTRAPOLATION 

EXTRAPOLATION  IS  DESIRED 

c 

ND 

s 

NUMBER 

OF  DIMENSIONS  (WHPN  7=FfX,Y), 

r' 

u 

CD 

z 

c 

NA1 

NO. 

OF 

VALUES  FOR  FIRST 

INCEPFNOENT 

VARIABLE 

c 

HA  2 

“ 

NO. 

OF 

VALUES  POR  SECOND 

independent 

variaple 

c 

NA3 

NO. 

OF 

VALUES  FOR  THIRD 

TNCEPENDENT 

VARIABLE 

c 

NA4 

NO. 

OF 

VALUES  FOR  FOURTH 

INDEPENDENT 

VARIABLE 

c 

NAS 

r 

NO. 

OF 

VALUES  FOR  FIFTH 

INDEPENDENT 

VARIABLE 

c 

p 

NA6 

- 

NO. 

OF 

VALUES  FOR  SIXTH 

TNCEPENDENT 

variable 

u 

c 

**************** 

C0MM0N/T3LKU°/Li 

************************************************ 

,LF,NA (6) , XL (ICO) ,NNEX 

DIMENSION  XU) 
llTOT(S) 

DO  10  1=1, NA1 
10  XKI)sX(I) 
K=NA1+1 
L=NA1+NA2 
DO  12  I=K, L 
12  XL(I)=X(I) 
NNEX=MEXTR 
NA  (1) =NA1 
NA  (2) =NA2 
N  A  (  3)  =N A 3 
NA  ( 4 ) =N A 4 
NA (5) =NA5 
NA ( 6) =NA6 
Ll=2 
LF=NO-l 
RETURN 
END 


X  A  (6)  ,  NS  (  5)  ,  WJ  (32)  ,  RAT  10  ( 5)  ,  NGROUP(S), 
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Appendix  E 


A  Sample  Problem 


Input  Data 


Visibility 

at  sea  level 

10  mi. 

Krtter  vapor 

pressure 

at  sea  level 

8  maHg 

Albedo 

.20 

Haze  layer 

above  sea 

level 

10,000  ft 

Vehicles 

x,y,2 

Position 

x,y,z  Velocity 

(ft) 

(ft/sec) 

1 

X 

3000 

Vx  6000 

y 

0 

Vje  100 

z 

2000 

V  0 

Z 

2 

2000 

600 

100 

0 

1500 

0 

3 

1000 

400 

100 

0 

1500 

-500 

4 

1000 

800 

0 

0 

20000 

0 

5 

1500 

1000 

1000 

0 

8000 

0 

Vulnerabilities 


Thermal 
Overpressure 
Dynamic  pressure 


10  cal/cm' 
psi 
r»  psi 


Height  of  ground  above 

sea  level  1000  ft 


Bursts 

1 

2 


x,y,z  Position  (ft)  Yield  kT 

x  11,000  600 

y  0 

z  35,000 

Targeted  10 

(spe  ■  800  ft) 


Time  (sec) 

0.0 


15.0 


100 
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OLtjfit  Listing 


VCHICLF  INITIAL  CONCITICNS 


VEHICLE 

PCSITION(X,Y 

,7) (FT) 

VELOCITY (X, Y,  2) 

(FT/SFC) 

1 

o  .-ail  J  C 

-  j .  r .. 

2  J  *  3 .  j  T 

60G  »C  j 

ic:.  u 

-C.  ;3 

2 

p  /*  r 

Cv  X  •  •  » 

ir:.cu 

is:-...,** 

6Cr.r J 

-c. .  0 

**  L  •  j  " 

T 

V 

1  j  C  <-  •  0  . 

1  -  U  •  l  c 

1  So  j  •  3  ^ 

ACC. 03 

-r.'.u 

-5  0  a • o . 

4 

l :  *j  ; . )  r. 

•  r  •  U  w 

2,30.  .  j  j 

8  0C.C  ’ 

-0.rl 

*  O  •  ‘J  J 

5 

1 F  i  C  •  •:  V 

1  j  0  u  .  v 

8 ;i j. ju 

uof  .o ) 

-l  •  Uli 

“  C  •  v  0 

NPlJT  PA  ?APrT Ppc-- 

NU-TER  OF  a,,RS t$  EflTEPE'*  I* 

2 

VISIBILITY  AT  SEA  LEVEL (MI.)  IS 

10.03 

ALBc')0  (GROMNC  RrFLECTANCp)  IF 

.2*“ 

H20  VAPOR  PRESSURE (FF A  LEVEL, MMHG) 
STAN'JAC*C  A T NOSPMEpF  SPECINIFH 

8.»'1 

BURST  NtWDCR  1 . 

vehicle 

POSITIONS, v 

,■»)  (FT) 

VELOCTTY(X,Y,7) 

(FT/SEC): 

1 

T  1  »  L  «  J  " 

.  .UV 

2  ,  0  .  A  j 

63C.-0 

1J0.CC 

r 

-  .  .  w 

2 

2  u  V  '  .  it. 

1  W  •  ••  V 

15. ".33 

5  c  W  .  0  a 

-o. 'C 

_  «  i 

~  •  *.  U 

3 

1  J  o  l  .  0  r 

lw-T.cr 

15^  i. Jl 

400.0 J 

•0  ♦  .u 

-5  j  C.  jO 

4 

.  i  j  -  c .  •) : 

3  •  5  % 

20  iui.JJ 

Sui .00 

-c.lg 

-  w .  j  i 

5 

153t.DC 

500  .00 

10  u  G  •  C  0 

■t.  Co 

_  o  -.! 

•  •  •  J  -  i 

WEAPON  PURST  =OSITIOM5D  RY  USp’--~ - 
NUCLEAR  INTONATION  COORPINATCS  APE 
GPOUNr  HEIGHT  AaOVr  SEA  L“VtL= 
INCFEMcNT A L  i*IRr  OF  nUPST= 


11.00.0, 
ICJw.O  FT 

-;.o  sec 


•••Of  350  0  j»  3 


LOSSES  FROM  BURST  NUMBER  1  (  630. EKT) 
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RESULTS  from  thermal  EFFECTS  COMPUTATIONS 

THERMAL  EFFECTS  PARAMETERS  FOR  VEHICLE  NUM7ER  1 


UNATTENUATEO  FL'JENCE  (CAL/CM2)  12.5574 

HEIGHT  OF  SUP^T  (FT)  350- 0. Cl 3r 

SLANT  RANGE  AT  TIMP  OP  PEAk  RADIANT  POWrR  (MI.)  6.4193 

TIME  TD  PEAK  RADIANT  POW£p  *SEC)  .44h4 

DIRECT  FREE  FIFLO  FLUEHCE  CCAL/CM2)  8.J913 

DIRECT  N C°M ALLY  INCIDENT  FLl,cNCE  (CAL/CM2)  7.8778 

(NORMAL  TO  A  PLANE  H0RI7  0NT AL  RECEIVE-"') 

THERMAL  EFFECTS  PARAMETERS  FOR  VEHICLE  NUM">Ep  2 

UNATTENUATED  cLUdNCE  (CAL/CM2)  12.7364 

HEIGHT  OF  PURST  (FT)  35U.C.QL0J 


SLANT  RANGE  AT  TIKE  OF  PEAK  PADIANT  POWE^  (MT.)  6.5563 

TIME  TO  PEAK  RADIANT  POWER  (SEC)  .4444 

DIRECT  FFEE  FIELD  FLUpNCF  (CAL/CM2)  7.5485 

DIRECT  NORMALLY  INCIDENT  FL'JENCE  (CAL/CM2)  7.3044 

(NORMAL  TO  A  PLANE  HORIZONTAL  RECEIVER) 

THERMAL  EpFECTS  FARAMETEFS  PQR  VEHICLE  NUMBER  3 


ONATTFNUATEC  FLUENCE  ( CAL/CM2)  11.6936 
HEIGHT  CF  DURST  (FT)  350-0. OHuO 
Scant  RANGE  AT  TIKE  OF  PEAK  RADIANT  POWER  (MI.)  6.6522 
TIME  TO  PEAK  PADIANT  POWE°  (SEC)  .4444 
DIRECT  FFcE  FIELD  FLUFNCE  (CAL/CM2)  7.2296 
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CIRfCT  NORMALLY  INCIDENT  FLUEMCE (CAL/CM2)  6.9411 

(NO0  iAL  70  A  °LANL  HORIZONTAL  RECEIVE3) 

f.’HICLE  NUMR-P  4  HAS  BE-N  SUBJECTED  TO  THERMAL  PLUENCE  OF 

4~.7344(CAL/CH?)  AM J  MAS  3rL N  3EMOVED  F°OM  THE  PROGRAM 

- THE  THERMAL  VULNERABILITY  LEVEL  I3  1 j.  JGL  u  (CAL/CM2) 

THERMAL  EFFECTS  PARAMETERS  FOR  VEHICLE  NUMBER  4 


UNA TT EMU A TEC  FLUENCE  (CAL /CM 2)  45.3626 

HEIGHT  OF  BURST  (FT)  350CCo GCuj 

SLANT  R A Nl> E  AT  TIME  OF  PEAK  RADIANT  POWER  (MI.)  3.3775 

TIME  TO  r/ftK  RADIANT  POWER  (SFC)  .4444 

DIRECT  FREE  FIELD  FLUENCE  (CAL/CM2)  43.7344 

DIRECT  NORMALLY  INCIDENT  FLUrNCE  (CAL/CM2)  26.8237 

(NORMAL  TO  A  PLAN0  H0°I70NTAL  RECEIVE3) 


— —  —  —  V EH ICLE  NUMRE3  *  HAS  BEEN  SUBJECTED  TO  THFRMAL  FLUENCE  Oc 
15. 14i:  (CAL/CN?)  AND  HAS  SEEN  REMOVED  FROM  TEE  PROGRAM 
- - JWE  jhc?mAL  VULNERABILITY  LEVrL  IS  Iu.JjTO  ( C A L / 0 M 2 ) 


THERMAL  EFFECTS  PARAMETERS  FOR  VEHICLE  NUMBER  5 


UNATTENUATED  FLUENCE  (CAL/CM2)  17.766' 

HEIGHT  OF  n UR*' T  (FT)  35i.eC.30u0 

SLANT  RANoE  AT  TIME  OF  PEAK  PAOIANT  PGWER  (MI.)  5.3969 

TIME  TO  PEAK  RADIANT  POWER  (SEC)  .4444 

DIRECT  F°EF  FI'LD  FLUENCE  (C AL/CM2)  15.1413 

DIRECT  NC3M ALLY  INCIDENT  FLUENCE (CAL/CM2)  14.3466 

(NORMAL  TO  A  °L ANF  HORIZONTAL  RECEIVE3) 
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RESULTS  FROM  3LAST  EFFECTS  COMPUTAT ION'S 


SHOCK  FROM!  p/^AMEI  6=3  AT  INTERCEPT  0e  VEHICLE  NUMRPp  i 
VEHICLE  POSITION  AT  SHOCK  INTERCEPT ( X , Y , 2)  21259. 


3043. 


RECEIVER  HEIGhT  AT  SHOCK  INTERCEPT  (FT) 


2C00.0JU-' 


SLANT  R A N&r  AT  SHOCK  INTERCEPT  (FT) 


34691. 5577 


TIME  OF  SHOCK  ARRIVAL  (SEC) 


15.239* 


SHOCK  FCCNT  VELOCITY  (FT/SEC) 


1139.9193 


PEAK  DYNAMIC  FRES?U°F  (PSI) 


.0214 


PEAK  OVERPRESSURE  (PSI) 


.9085 


PEAK  MATERIAL  VELOCITY  (FT/SEC) 


51.2153 


PEAK  OVERDPNSITY  (SLUGS/PT ** 3) 


•  0-J  23 


POSITIVE  DURATION  OF  OVERDENSI TY  (SFC) 


2.7949 


POSITIVE  DU°ATION  OVErDRESSURE  (SEC) 


2. 8494 


POS IT IV z  DURATION  MATERIAL  VELOCITY  (SEC) 


2.4096 


SHOCK  FRONT  PARAMETERS  AT  INTERCEPT  OF  VEHICLE  NUM°ER  2 
VEHICL'P  POSITION  AT  SHOCK  INTERCEPT  (  X  ,  Y ,  Z)  2G268. 


10  0. 


RECEIVER  HEIGHT  AT  SHOCK  INTERCE°T  (FT) 


1500.0000 


SLANT  RANGE  AT  SHOCK  INTERCEPT  (FT) 


34758.4690 


TIMF  OF  SHOCK  ARRIVAL  (SEC) 


15.1483 


SHOCK  FRGNT  VELOCITY  (^T/SEC) 


1141.5747 


PEAK  OYNAMIC  PRESSURE  (PSI) 


,0213 


PEAK  OVERPRESSURE  ( I ) 


.9156 


PEAK  MATERIAL  VELOCITY  (FT/SEC) 


50.7852 


PEAK  OVEFUENSITY  <SLUGS/FT**3) 


0G  24 
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•\r 


PCSITI Vc  DURATION  OF  OVERnENSITY  (SEC) 
POSTTI'/R  DURATION  OVF  PoPESSU^f  (SEC) 
POSITIVE  OU PATTON  MATERIAL  VELOCITY  (SEC) 


2.7743 
2. 829- 
3*3841 


•VEHICLE  NUM9£D  1  MAS  INTERCEDED 
GROUND  LEVEL  AND  MAS  nEEN  REMOVED. 


■VEHICLE  POSITIONS  AT  TIME  OF  3’JRST  NUM°£R  2- 


VEHICLE 

1 

2 


PCSITION(>',Y,7)  (FT) 

12  j  0 l ♦ 3  0  IEcj.vL  2  J  0  1.30 
1 1 )  P  i  .  u  3  luO.C:  15--..  33 


VELOCITY (X,Y,2) 
60 -.CO  1QC.-0 
6o-j»un  -C  •  .  u 


(FT/SEC) 

n  ** 

—  •  J  . 

_  -  rv  n 
•«*  •  +>  • 


*  * 
*  *4 


■WEAPON  DURST  POSITIONED  8Y  SUOROUTINF  TARGET - 

THc.  TARGFT  VEHICLE  1*  VEHICLE  It  1 
GROUND  HEIGHT  A^OVE  SEA  LE  Vl  L=  1-3  .  f.  RT 

INCPE MENTAL  T I Mr  OF  BURST=  15.3  SEC 

8URST  OPTCNATICN  AT  5'59.179FT  FROM  TARGET  VEHICLE 
NUCLEAR  OFTONATION.  COORDINATES  A°E  12555.8  1504.6 


2195.3 


LOSSES  F°CH  9URP-T  NUMBER  2  ( 


lu.CKT) 


RESULTS  FRQ*-'  THERMAL  EFFECTS  COMPUTATIONS 

---VEHICLE  NUMBER  1  HAS  °EcN  SU3JECTFD  TO  THFPMAL  FLUENCE  OF 
1147,9366  (C A L/CM2)  ANC  HAS  BEEN  REMOVED  FROM  the  PROGRAM 
- THE  THERMAL  VULNERABILITY  LEVEL  IS  lO.uUbu  (CAL/CH?) 


CHE/ PH/ 72 -3 


THERMAL  EFFECTS  F4CAMET f F.S  FOR  VEHICLE  NUM°£P  1 


UNATTEMUfTFn  FLUENCE  (CAL/CH2)  12..2.998'j 

HEIGHT  CF  9U°3T  (FT)  21S5.32F2 

SLANT  RA^GE  AT  TXhE  OF  PEAK  PAOIANT  POWER  (Ml., )  .0931 

TIME  TO  PEAK  °ADI6NT  POWER  (SEC)  .1266 

DIRECT  F°Et  FIELD  FLUrNCE  (CAL/CM2)  1147*9365 

DIRECT  NORMALLY  IMCTOFMT  FLUENCE (CAL/OM?)  422.7175 

(NORMAL  TO  A  CL ANF  HORIZONTAL  RECEIVE0) 

REFLECTEC  FLUcNCE  (CAL/C^* *2)  .1733 

REFLECTEE  ♦  OIF-CT  NORMAL  pL UENCE  (CAL /CM 2 >  G.OOL'u 

(CON°UT  cn  ONLY  Ic  PU°FT  MELON  RECEIVER) 


—— VEHICLE  NUMrEr  2  HAM  ~F£N  SU3JECTFD  TO  THERMAL  FLUENCE  OF 
59. 8434 (CAL/CM?)  AM 0  HAS  9 SEN  REMOVED  FROM  THE  PROGRAM 
- -THE  Tmcjmal  VULNc°A°ILITY  LEVEL  IF  Iu.ClC.  (CAL/CM2) 


THERMAL  EFFECTS  PARAMETERS  FOR  VEHICLE  NUM °ER  2 


UNATTENUATEO  FLUEMCE  (CAL/CM2)  69.5182 

HFIGHT  OF  DURST  (FT)  2195.3262 

SLANT  RANGE  AT  TT  ME  OF  PEAK  RADIANT  POWER  (MI.)  .4082 

TIME  TO  PEAK  RADIANT  POWER  (SCO)  ,1266 

QIFFCT  FFEE  FIELD  FLUENCE  (0 AL/CM 2>  59.8484 

OIPECT  NORMALLY  INC 10° NT  FLUENCE  (CAL/CM2)  19.3056 

(NOPMAL  TO  A  CLANE  HORIZONTAL  RECEIVER) 

REFLECTEC  FLUENCE  (CAL/CM**?)  .0664 

REFLECTED  +  DIRECT  NORMAL  FLUENCE  (CAL /C M2)  0.0000 

(COMPUTEC  ONLY  IF  3UPST  DEL GW  RECEIVER) 


‘RESULTS  FROM  3 LAST  EFFECTS  COMPUTATIONS 
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